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Abstract 


It is first demonstrated that in a forecast with the primitive equations, pure geostrophic or 
purely non-divergent initial wind fields will result in noticeable meteorological “noise”, especially 
in a baroclinic atmosphere. It :s then shown in a linear example that this noise is greatly suppressed 
if the divergence of the initial wind field is set equal to the divergence implied by the usual 
geostrophic theory of numerical weather prediction. The absence of noise in Charney’s 1955 
test computation with non-divergent “balanced’’ initial data is explained. 

Finally, the quasi-geostrophic expansion is used to define initial data for a non-linear forecast 
with the primitive equations. Two cases are considered, that where only the geopotential is 
known from observations, and that where only the non-divergent part of the wind is known 
from observation. It is shown that the initial tendencies computed from initial data of this kind 
preserve the assumed geostrophic character of the flow, and thereby make it possible to get more 


accurate forecasts with the primitive equations. 


I. Introduction 


It is well known that the “primitive equa- 
tions” (the Eulerian hydrodynamical equations 
modified by the assumption of hydrostatic 
balance) are very sensitive to the initial (t=o0) 
fields of wind and pressure; unless these fields 
are suitably adjusted to one another, the fore- 
cast may contain primarily noise in the form of 
high-frequency gravity-inertia waves, and the 
more important motions of lower frequency 
will thereby be hidden. HINKELMANN (1951) 
has shown that the amplitude of this unwanted 
noise can be made smaller than that of the 
desired low frequency motion by simply using 
the geostrophic values for the initial wind 
field. His theoretical treatment was for a linear 
system, however. It was later shown in a non- 
linear example by CHARNEY (1955) that better 
results could be obtained by relating the initial 
wind and pressure (geopotential) fields by the 
“balance equation”: 
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Here the initial wind v is given by the stream- 
function y, so that v is initially non-divergent. 
The symbols / and 9 are the Coriolis parameter 
and the geopotential, respectively, and v is the 
horizontal gradient operator. (1) is derived 
from the equations of motion by simply re- 
quiring not only that the horizontal divergence 
v-v vanish at t=o, but also that 9 (v- v)/at 
vanish at f=o. 

It is easy to demonstrate that, in a linear sys- 
tem, the relation (r) is not sufficient to elimi- 
nate gravity-inertia waves. In a linear system 
(1) reduces to 


v.Ivv=v?o, (2) 


which, if lis considered constant, is simply the 
geostrophic relation. We will now consider 
a particular simple linear system and show that 
geostrophic initial conditions, although setting 
an upper limit to the amplitude of the mete- 
orological noise, are not precise enough, and 
that the initial wind field must contain some 
divergence if the gravity-inertia waves are to 
be satisfactorily eliminated from the forecast. 
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2. Noise from geostrophic initial data 


As a simple linear system we take that studied 
by HINKELMANN (1951). À constant basic current 
U (from west to east along the x-axis), in a homo- 
geneous incompressible atmosphere of average 
height H, has small perturbations superimposed 
on it. The perturbations are assumed to have 
infinite extent from north to south (y), and the 
Coriolis parameter | is considered a constant. 
Three types of wave motions are found to be 
possible; a “meteorological” wave moving 
approximately with the speed c, & U, and two 
gravity-inertia waves moving approximately 
with the speed 3 ~U+ V gH+ [?/r?. (vis the 
wave-number =2x/L, where L is the wave 
length in the west-east direction.) 


In general all three waves will be present. 
Following Hinkelmann, we let the total ampli- 
tude in the geopotential be 1, and let Ay, Ay, As 
be the amplitude in the geopotential associated 
with each of the three different waves (i. e. 
A,+A,+A,=1). We now wish to know the 
values of A,, As, and A, when the initial wind 
field is geostrophic. This is easily determined 
by setting Hinkelmann’s parameters 6 and ¢ [in 
his equations (16) and (14b)| equal to zero. 
Substitution of the resulting values of F and G 
in his equation (15) gives the desired result: 


A, =1+ O(w?), 
I [2 =e ’ (3) 
Arm -Asu|ıtaen| + O (u?) 
Due (4) 


The non-dimensional term /?/y?.¢H in (3) is less 
than 1, even for very long waves. 


The use of strict geostrophic initial conditions 
for the primitive equations in this simple system 
will therefore result in some noise. This noise 
will have an amplitude (in @) equal approxi- 
mately to u times the amplitude of the low 
frequency wave. In the barotropic system we 
have used, wis indeed small (u ~ 1/19). However, 
in a baroclinic system, with internal gravity 
waves, the value of ¢ is effectively reduced, and 
u will not be so small. Geostrophic initial con- 
ditions will therefore no longer be even parti- 
ally satisfactory for predictions with the pri- 
mitive equations. 
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3. Elimination of the noise in a linear 
system 


We now ask for a better definition of the 
initial u- and v-fields in Hinkelmann’s linear 
barotropic model—better in the sense that A, 
and A, will be smaller than the values given by 
(3). From Hinkelmann’s analysis we readily 
find that the u- and v-fields (u, and v,) associ- 
ated with the meteorological wave (subscript 1) 
are related to the geopotential field (g,) of that 
wave by the following equations: 
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VgH + Rp? VgH + P/y? 


1= 


iv 


= [+0] 781 - 


If the (total) initial v- and v-fields are related to 
the (total) initial @-distribution by the exact 
forms of (5) and (6), only the meteorological 
wave will be present, and there will be no noise 
motion. If the approximate forms of (5) and 
(6) are used, an analysis similar to that in section 
2 will show that A, and A, are at most of order 
u? instead of u. 

In the preceding analysis in section 2 both the 
initial v- and v-fields were assumed to be given 
geostrophically by the initial @-field. As we 
saw, this did not eliminate the noise completely, 
but only reduced it to something like 10 % of 
the geostrophic wave motion. (5) and (6) on 
the other hand show that only v in this linear 
system should be geostrophic if noise is to be 
eliminated to order w?2. The initial u-field 
should not be geostrophic (ugo=0o since 
dp|dy=0 by assumption), but should be given 
instead by (5). 

Since the entire analysis has been made assu- 
ming that du/dy,dv/dy and dp/dy are zero, the 
horizontal divergence is given by 

v= ou de dv = ou Se 
ae Oy Sone na 
The approximate form of the relation (5) 
between u and @ is therefore equivalent to 
requiring that 
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Thus, the elimination of the noise motion 
requires a certain amount of divergence in the 
initial velocity field. Although this has been 
demonstrated here for a linear system only, it 
is obvious that the presence of non-linearities 
will not make this needed divergence zero, but 
on the contrary will only make the initial- 
value problem still more complicated. We may 
therefore conclude that the use of “‘balanced’’ 
initial data, as suggested by Charney, is not 
sufficient to eliminate noise motions in general. 

We will now show that the initial field of 
divergence given by (7), whose introduction is 
necessary to eliminate the unwanted gravity- 
inertial waves in this linearized system, is 
exactly that field of divergence implied by the 
usual quasi-geostrophic vorticity equation. The 
quasi-geostrophic equations for this simple mo- 
del are (we still neglect 6 =dl/dy): 
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ae only aN = 0; 
8 
= +gHv-v=o, (8) 
€ =dvlax = 1-192 p/dx?. 


These may be solved for v - v by eliminating 
Ip | dt. 


We get 
9°? Fy op 
2 u : =. IE PL 
[ sh | RE Tr 
or 
Gh aie Op 
ee 7 —c ) 


gH+P/y? 0x  \gH+ Pr? 9x 
This relation, by its agreement with (7), shows 
that the initial field of divergence, which is 
necessary to eliminate the gravity-inertia waves 
in this linear initial-value problem, is just that 
field of divergence which is implied by the 
usual quasi-geostrophic system of forecasting. 
Although the analysis has been confined to a 
linear system, we may infer that the same pro- 
cedure should be followed for a flow pattern in 
which non-linear terms are important.! 


1 HINKELMANN (1959) has recently performed a 
numerical experiment verifying this result. 
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4. The absence of noise in Charney’s test 
computation 


The analysis above has shown that the more- 
or-less complete elimination of noise for the 
primitive equations is possible only if the initial 
wind field contains a divergence field equal to 
the field of divergence implied by the quasi- 
geostrophic theory. This seems to be in disa- 
greement with Cuarney ’s result (1955) that the 
use of the balance equation (1) with a non- 
divergent initial wind field is necessary to eli- 
minate the noise. We will see that there is no 
disagreement because the idealized flow pattern 
used by him to test his hypothesis was a special 
case. 

Charney’s initial flow pattern was of the 
simple form 


TEN 
y=Asin sin 


where w is the stream function, and A, Land D 
are constants. This flow pattern existed in the 
upper layer of a two layer model atmosphere, 
with the lower layer being so deep that its 
motions could be ignored. 

The (non-linear) quasi-geostrophic equations 
for this model may be written: 


le a Iv-v=o, 


ea (11) 


(12) 
(13) 


(The Coriolis parameter | was also taken as 
constant by Charney.) Here g’ is a reduced 
value of gravity equal to (1-—¢’/e) g, where 0’ 
and o are the densities of the upper and lower 
layers. These equations are readily combined 
to give an equation for v : v: 


ee pale gad ee 1 


This equation gives the quasi-geostrophically 
implied value of v - v which, according to the 
results of section 2 and 3, should be included 
in the initial wind field if no noise is to appear. 

It now becomes clear why Charney’s neglect 
of this initial divergence field did not give rise 


ie +gHv-v=o, 


= v’y 
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to gravity-interia waves in his computation, for 
it is apparent from the simple form of (10) that 
the right side of (14) vanishes. This being the 
case vv also vanishes, and it was therefore 
unnecessary to include divergence in the initial 
wind field in his special case. On the other hand, 
the results of Charney’s test do show conclusi- 
vely that, in a non-linear flow pattern, the 
relation (r) is better than the geostrophic rela- 
tion as a constraint on the initial data. Evidently 
(1) gives a better definition of the initial vorticity 
than does the geostrophic formulae. 


5. Initial data for non-linear problems 


When the equations of hydrodynamics are 
not linearized, it is no longer possible to separate 
the different types of motion, and by a clever 
choice of initial conditions, to exclude forever 
any unwanted “noise” from the forecast. As a 
compromise, however, we can require this of 
the initial data—that when at t=0 we compute 
the local rates of change of dv/dt, Ap/dt, etc., 
the latter quantities should all have values very 
close to what they would have (initially) in a 
quasi-geostrophic forecast. This is a reasonable 
reqirement first of all because we know that the 
geostrophically computed values are at least 
approximately correct. Secondly, as is apparent 
from section 3, this requirement does eliminate 
the unwanted noise almost completely in a 
linear system. 


A technique for imposing this requirement 
comes very neatly out of the familiar theory 
of the quasi-geostrophic expansion (PHILIPPS, 
1939). So far the most useful application of 
this theory has been to derive the well-known 
geostrophic system of prognostic equations 
(ELIASSEN, 1948; IuDIN, 1957; KIBEL, 1957; 
Monin, 1958). Here we will emphasize its 
usefulness in deriving satisfactory initial data 
for the primitive equations. In the process we 
will see that the appropriate initial values of 
v + vand k- vxv— the importance of the 
former being apparent from section 3 and the 
importance of the latter from Charney’s 1955 
computation—will both appear as improve- 
ments on the use of geostrophic initial data. 

We first set the horizontal wind v equal to 
geostrophic wind v, plus the geostrophic devia- 
tion v’, 


’ 
v=v+V, 


(15) 
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whereupon the horizontal equation of motion 
can be written 


Pulp dv 
v= 7k x at (16) 
k is the vertical unit vector. The acceleration 
dv dt is now approximated geostrophically, 


1 ACTE) 


and this equation is used to evaluate v : v’ and 
C=k-vxv: 


; 21.06 
V-v = -5| Git ve vee. (18) 
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The Coriolis parameter ] has so far been con- 


sidered as a constant. As the final step we com- 
pute also vv, and E,, now letting / vary, 


(20) 
(21) 


and combine these expressions with (18) and 
(19) to get: 


-v,= —I-1v,- 
Vev,= SITE wl 


IC, = v2 —1-1vl- vo, 


=o. I 
If = v?p--vœo.vl- 


2 [2 Pp {29 ÿ 
PTox?9y? \dxdy] I 


If the right side of these equations can be 
evaluated, the field of v itself can be computed 
by integration of the v : v- and £-fields. 

Let us for the moment take the simple view 
that the initial distribution of @ is known, while 
that of v is only poorly specified by the wind 
observations. In this case we consider (22) and 
(23) as determining more accurate values of v 
than are given by the hypothetically poor 
wind observations. Evaluation of the right side 
of (23) gives no difficulty. However, (22), 
which is the usual geostrophic form of the 
vorticity equation, contains the term ¢,/dt= 
I"! v29p/ot. This term we treat by introducing 
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the geostrophic and hydrostatic form of the 
thermal energy equation for adiabatic motion, 


d (op 0 bs 
Herden cn 


and the continuity equation, 


v-v+owldp=o. (25) 
[¢ =0(p) is normally taken as a standard value of 
(Ap/p) x (2 In O/dp) where © is the potential 
temperature.] (22), (24), and (25), together 
with appropriate boundary conditions are, in 
fact, the quasi-geostrophic prognostic equations 
and can be solved for both 9p/dt and v - v. We 
are therefore able to determine v, since both & 
and v - v are now known. 

Let us now assume that the initial data for a 
forecast with the primitive equations is the 
field of v we have computed in this way 
together with the original field of y. We 
examine first the initial value of 0v/dt computed 
from the original equation of motion (16), 
realizing that v’ in (16) is now equal to the 
value given it by (17). 


ov ov : 
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The last step here is possible because, as shown 
by Cuarney’s scale theory (1948), wov/op is 
smaller than any of the other terms in (26) for 
geostrophic types of motion.! The difference 
between the two horizontal advection terms is 
also smaller (presumably) than the individual 
terms. Assuming finally that 0/at is of the same 
order of magnitude as v- v, we end up with 
only dv,/dt=!k x v(Ap/dt) on the right side. 
[[= constant here, as it was in (17).] We may 
then conclude that the initial values of dv/dt are 
such that 


Oe 1.102 9 
SE at 


(27) 


A EL (28) 


1 This is not true, though, for extremely long waves 
(Burcer 1958). For such waves, however, dv/dt = o. 
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We now investigate the value of 0(dp/dp) /at 
computed with the complete (i.e. non-geo- 
strophic) form of the thermal energy equation: 
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On the right side of this equation v is again the 
v determined from (22) and (23). w is identical 
with the w appearing in (24) and (25), as is also 
the case with v(2p/0p). o, however, is presum- 
ably the actual value of (dp/dp) (0 In @/dp) in- 
stead of a standard value as in (24). Applying a 
subscript g to the value of o(op/op)/dt com- 
puted geostrophically on the left side of (24), 
and eliminating the common value of w be- 
tween (29) and (24), we obtain 


Le) (EC) 
[Of Lo 


Here the final approximate equality is true only 
if (c/o) is not very different from 1. 

The importance of establishing the final 
approximate equalities of (26) and (30) is that 
they show that, although the fields of v and y 
are changing with time, they will remain in 
approximate geostrophic balance.‘ This, of course, 
is the implication of the approximation used 
between (16) and (17). 

Equation (23), which gives an improved value 
for the initial vorticity, is the same as the balance 
equation (1), except that the nonlinear terms 
in y and the term vy : vl have been “‘approxi- 
mated” by their geostrophic equivalents. As 
Mon has shown (1958), (23) and (r) are both 
of the same accuracy as far as the quasi-geostro- 
phic expansion is concerned. 

Let us now consider the opposite situation, 
namely that situation where the wind observa- 
tions are “better” than the observations of y. 
As CHARNEY has demonstrated in his scale ana- 
lysis (1948) ¢ is larger than v - v (again, except 
for extremely large scale motions) so that we 
may assume that ¢ is given satisfactorily by the 
wind observations, but that v - v is not. In this 


1 A more complete demonstration would require 
also an investigation of the boundary condition on the 
pressure at the surface of the earth. To save space 
it is not reproduced here, especially since there are 
several alternate ways of formulating this condition. 
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case we take the balance equation (r) to de- 
termine the unknown 9, 


Jy Jy Jy if 
an = Js ET 
V2q = V lop +2| 54 y2 (a aD) 


after y has first been determined by solving the 
Poisson equation 


veyp=C=k-v X Vobserved ° (31) 
v-v has still to be determined. We determine 
it again from the quasi-geostrophic system ex- 
cept that we use the non-divergent wind v,, = 
=k x vy in place of the v, in (22): 


Iv.v= - rm 4 29) (32) 


The counterpart of (24) is written 


9 [op Au oy x 
(2) neo (2)-m on 


by making the substitution y— 1%. The con- 
tinuity equation (25) is still retained. Both 
dy/dt and v - v can now be determined from 
(32), (33) and (25), just as 2y/0f and v : v could 
be determined from (22), (24) and (25). Finally, 
it is again possible to show that the approximate 
relations (26) and (30) are valid in this case also. 

In actual practice, both wind and geopotential 
observations are available, and it would seem 
advisable to somehow combine the two analy- 


sis approaches described by (22), (23) and (24) 
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on the one hand and by (r), (32) and (33) on 
the other hand. An extension of the type o 
procedure suggested by Sasakı (1958) may be 
very useful in this connection. 

The above methods of determining “good” 
initial data for the primitive equations are 
clearly based on the assumption that the wind is 
nearly geostrophic. The extent to which the 
method will prove useful in low latitudes, or in 
high latitudes when vv approaches / in magni- 
tude, can perhaps be answered only by experi- 
ence. 

In closing, it should be emphasized that the 
approximate equalities of (26) and (30) are 
not goals to be attained for their own sake. The 
whole idea of using the primitive equations in- 
stead of the usual geostrophic system is to obtain 
an increased accuracy, not merely to duplicate 
the results of the latter method. This increased 
accuracy in Jv/dt, for example, is to be obtained 
by including a term like w/v/dp in (26), a pro- 
cedure which is not possible with the geostrop- 
hic system. However, since wdv/dp is, by the 
scale theory smaller than (v- v)v or Ikxv’, 
it does no good to include wdv/dp unless the 
larger terms (v - v)v and /k x v’ are at the same 
time evaluated with sufficient accuracy them- 
selves. This may be considered as the fundamen- 
tal reason for being so careful with the initial 
data for the primitive equations. 
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Disturbances 


By ROBERT G. FLEAGLE, University of Washington! 2 


(Manuscript received October 27, 1959) 


Abstract 


Effects of heating on large scale disturbances are studied by assuming that heating is tuned 
to the disturbances, with rate of heating proportional to individual pressure change but displaced 
by an arbitrary phase angle. The linearized two-layer baroclinic model reveals that the heating 
which accompanies condensation in the middle troposphere may lead to marked increase in 
instability and to shift of maximum instability from five to three thousand kilometer wave- 
lengths, and this calls attention to a serious defect in the two-layer model which must influence 
the accuracy of numerical calculations using the non-linear equations. 


Introduction 


MILLER (1946) and MILER and Mantis 
(1947) have shown that the oceanic regions 
off the east coasts of North America and Asia 
are highly favored regions for cyclogenesis 
in winter. These are of course regions of large 
air-sea temperature difference so that the 
conclusion is almost inescapable that heat 
provided by the warm ocean plays an impor- 
tant role in cyclogenesis. For these reasons 
investigation of the distribution of heating 
within large scale disturbances and of the 
effects of this heating is needed. 

A general discussion of the distribution of 
heating and its effects is not possible because 
the physics of the essential transfer processes 
as well as the physics of clouds is incompletely 
understood. Nor is it possible accurately to 
measure the relevant quantities at many 
oints and so to map heating over volurnes 
of the size of large scale disturbances. It is 
possible, however, to sit in a comfortable 
chair and to make certain useful generali- 
zations which can serve as a starting point 


1 Contribution No. 49 Department of Meteorology 
and Climatology, University of Washington. 

2 Written while National Science Foundation Senior 
Post doctoral Fellow at the Imperial College of Science 
and Technology, 1958—1959. 
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in investigating the problem. One realizes 
that the major source of heating in the free 
atmosphere is the latent heat released in con- 
densation within clouds, one realizes that this 
heating is directly and rather simply related to 
vertical displacement, and one can make 
numerical estimates of heating to better than 
order of magnitude accuracy. One realizes, 
too, that there must be an imperfect correla- 
tion between meridional wind velocity and 
transfer of heat and water vapor at the ocean 
surface. Latent heat is distributed through a 
deep layer, perhaps the lower half of the 
atmosphere on the pressure scale, and reaches 
its maximum value at, say, 700 to 900 mb. 
Heat transferred to the atmosphere by contact 
with a warmer surface also is distributed 
through a deep layer but probably decreases 
in magnitude with height from the maximum 
value at the surface. 

Qualitative and crude though these generali- 
zations are, I should like to understand how 
heating of this sort influences the develop- 
ment of large scale disturbances. For this 
reason I have incorporated heating in a linear 
treatment of the large scale motion. The 
objective is to investigate the effects of various 
distributions of heating on the rate of growth 
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of unstable disturbances. Heating is assumed 
proportional to individual pressure change; 
this represents the physics of the condensation 
and evaporation processes with some accuracy, 
and it may be considered to represent in a 
purely formal way heating by conduction or 
by radiation so long as this heating is tuned to 
the developing disturbance. 


The Equation of Growth 


The most important theoretical properties 
of large scale disturbances in a baroclinic 
atmosphere are now well known for normal 
middle latitude conditions, at least for the 
linear case. Rate of growth increases with in- 
creasing dimension of the disturbance to an 
equivalent wavelength of about five thousand 
kilometres and decreases at greater wavelengths, 
and disturbances normally require a time 
period of the order of a day in which to double 
in amplitude. This is probably the most impor- 
tant result in theoretical meteorology; it 
reveals that large scale disturbances originate 
in initial impulses of large spatial dimensions 
and not in the inaccessible region of very 
small disturbances. The mathematical tour de 
force which affords this result is quite formi- 
dable (Kuo, 1952; GREEN, 1960) and there are 
aspects of the theory which have not yet 
been clarified. Nevertheless, it is possible to 
derive the results stated above by simplified 
approximate methods. It seems to be true, 
fortunately, that the approximate methods 
are most nearly valid for dimensions of the 
greatest growth rates, and these are the 
dimensions of greatest interest. 

In the following I have extended the develop- 
ment of the linearized two-layer baroclinic 
model given by ELIASSEN (1952) and Prius 
(1954) by including heating in the thermody- 
namic equation. The assumptions on which 
the two-layer baroclinic model rests have 
been discussed widely and they are summarized 
by PETTERSSEN (1956). The linearized vorticity 
equation may be written 


210) 


L (0) +v(8 70/0) +52 QUIOy ~f) =0 (1) 


where the symbols have the commonly ac- 
cepted meanings as follows: 
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x, y coordinates toward east and north respec- 
tively 

p pressure | 

u, v perturbation velocity components in x 
and y directions respectively 

Ë dviax —duloy 

U undisturbed velocity component in x 


direction . 

w individual rate of change of pressure with 
time : 

f Coriolis parameter 

B  ofloy 

L djat+Ualax 


The terms 22U/9y? and 9U/dy will be omitted 
in subsequent steps, but they may be rein- 
troduced at any point, if desired. The ther- 
modynamic energy equation may be written 


os ee 20 U2 aln@ | ain 
CpT dt 09 09% dy op 


(2) 


where @ and 6 represent respectively the undis- 
turbed and disturbed potential temperatures, 
T the undisturbed temperature, q the heat 
added per unit mass, and C, specific heat at 
constant pressure. 

In order to express the heating as proportional 
to the individual rate of pressure change (w) 
and to leave the phase relation (5) between 
heating and w arbitrary we may write 


LA, dg ORG 
oar meres he ig (3) 
re dq/d : 
where a= GT. a 6) dl The arbitrary quan- 


tities ‘a’ and ‘6’ are chosen to represent any 
one of several heating distributions. It follows 
from equation (3) that 


O-1L(8) +valnO]oy 


AlnO/dp (1 — ae?) 


(4) 


From the hydrostatic equation 6/0 = 
= (dz/dp)/9Z/Op where z and Z represent, 
respectively, the perturbation height and 
undisturbed height of the pressure surface. If 
this relation is substituted into (4) and if the 
derivative with respect to p is replaced by a 
finite difference, then the individual rate of 
pressure change at a surface 2, is expressible as 
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pares Ce 
RF (ps =p) Lis 20 + 
+ U À (2 22) | + 24 
ss a dy 
es IlnO r u 
| 5 (1 - ae’ ] 
3 
We make the following assignments of 


pressure surface: 


Zo : Zero pressure 
21,1250 mb 

22: 500 mb 

23: 750 mb 

24: 1000 mb 


If equation (5) were applied at z2, only the 
mean heating of the atmosphere but not its 
vertical distribution could be incorporated in 
the problem. Accordingly, the equation is 
applied at z, and 24. In order to avoid intro- 
ducing additional dependent variables, how- 
ever, we shall make the assumption that 
Zy—29 and 2,5—2, may be replaced in (5) and 
in the analogous equation for w, by 23—2, 
that is, that the temperature distributions at 2, 
z, and z, are sufficiently the same to be used 
interchangeably. Finally we introduce the 
geostrophic relations 


and for convenience assume 0U/dp independent 
of pressure and T, ~ T;~ T. Equation (5) then 
may be written 


gps[Ls(z3 - 21) + A U9z3/0x] (8) 
RTp; (9 InOJop)s(1 = ageis) 


O3 


Similarly 


_ gPı lLıl&s =) + AUD l0x] 0) 
Fous RTp;(o InOJop)ı(t - a, e's) 


These relations are to be used to eliminate w 
from the vorticity equation applied at p, and ps. 
Tellus XII (1960), 2 


We assume that dw/dp is directly propor- 
tional to and has the same phase relations as æ 
at each of the pressure surfaces. For simplicity 
the proportionality factors are set equal to 4/3p9; 
this permits an asymmetrical @ distribution 
in the general case but reduces to the much 
used parabolic profile in the case w, = wy. 
The finite difference representations for Jw/Op 
may be written for © = o at p, and py. 


(11) 


Finally, the geostrophic relations are introduced 
into the vorticity equation, so that at p, 


d d 
D [V22 + 223 -21)] + U, ER [V?2, + 


+4 (23 =2)] + (8 + A240) eo (12) 
and at ps 
d N) 
5, [V?23 — A(2s - 21)] + Us 3e [Vs - 
Ar? 24 073 = 
es — 24) ] + (B ~ 26840) 0 (13) 
where 
Aisha oy hon LLC a 
2 T4 3 p2RT . ( dp N aye ) (14) 
pe dd Vire (32 Br 
Az gep RT 5 op ne age ) (15) 


We consider only solutions in which the 
wavelength, speed, and rate of growth are 
independent of height. Therefore, let 


(16) 
(17) 


These may be readily generalized to include y 
dependence. Upon substituting into equations 
(12) and (13), a quadratic frequency equation 
follows: 


= 24 eixx-iot 


a = 2 eestor 
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The solution may be written 


DE et ll 


FE BETT et Pam LE es 
Er BAU 
(au (7-2). Une 
os )F 
Ta Ber (19) 
where 


For I = I's equation (19) takes the form 


N 
pe ee ia {4u*(x =") - 


AIDES 


(20) 


which is a slightly more general form of the 
frequency equation developed by Phillips 
(1954). Equality of I and J’; implies that 
À? = 3°, and this can occur in a variety of 
ways. For example, for a, = a3 and 6, = 0g, 
I’, equals 7’, if the static stability at z, is three 
times that at z,. Or, the same result may occur 
for equal static stabilities if 6, = 5, =o and 
I — a3 = 3(1 —4,); this implies that heating at 
z, exceeds heating at 24. 

For the case of no heating (a = 0) and for 
the case of heating in phase with w, J" is real; 
and the second term of equation (20) represents 
the growth of unstable baroclinic waves 
(positive root) as well as the reduction in 
amplitude of a damped wave (negative root). 
For the case of heating displaced in phase from 
w, I’ is complex; and equation (20) is diffi- 
cult to interpret except by specific examples. 
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Results 
Curve A in Figure 1 represents the case of no 
heating for the following values: 


B:1.6 10: sec, cm (lat. 43.deg.) 
MU 220% 10°. cm Sees, 


22: 1.46 x 10716 cm? (9 In O/dz at soo mb: 
16% 10° cm!) 


Maximum instability is reached at a wave- 
length of 5x 103 kilometers, and all waves 
shorter than about 3.6 x 10° and longer than 
about 7x 10% kilometers are stable. These 
results should be compared with the result of 
integrating the three dimensional linear equa- 
tions, as has been done by Kuo (1952), for 
example. The models are not identical, but 
this probably is of little importance. Curve K 
is taken from Kuo’s paper for JU/dz =3 x 10-8 
sec! and 9/n@/9z=1.5 x 10-7 cm-l. The 
sharp cut-off to instability on both the short 
and long wave sides shown by curve A is 
evidently incorrect, and the wavelength of 
greatest instability is about 10% kilometers 
too great in the two-layer model. Corre- 
sponding defects must be expected in the non- 
linear treatment of the two-layer model as well. 

The effect of heating in the two-layer model 
may be understood qualitatively from equa- 
tions (3), (14) and (15). Consider first the case 
which results from vertical motion in saturated 
air; condensation and heating occur with 
upward motion, evaporation and cooling with 
downward motion. For this case 6 vanishes. 

It is evident that as ‘a’ increases I” decreases 
and instability increases in magnitude and 
covers a wider range of wavelengths. For 
a = 0.7, which corresponds to heating at the 
saturated adiabatic rate, equation (20) yields 
curve B in Figure 1. Maximum value of the 
exponential factor is more than doubled, the 
wavelength of maximum instability shifts 
from 5x 10° kilometers to 3x 10% kilo- 
meters; and instability is extended to much 
shorter waves. This result may be achieved 
also by calculating static stability with respect 
to saturated air and setting ‘a’ equal to zero in 
equation (20). 

We may imagine the reverse distribution of 
heating in subsiding air, cooling in ascending 
air; although this situation is unrealistic for the 
normal middle latitude planetary disturbance. 
For this case 6 =z and J” is increased over the 
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adiabatic value. This results in decreased in- 
stability; for a = 0.7 maximum value of the 
exponential factor of only 10-8 sec-! occurs at 
about 6 x 103 kilometers. This is not shown in 
Figure 1. 


LS 


Factor (lo sec”) 


Exponential 


Wavelength (10° Kilometers) 


Fig. 1. Exponential factor as a function of wavelength 

calculated from equation (20) for no heating (curve A), 

for heating east of pressure troughs (curve B), and from 

equation (21) for heating in pressure troughs (curve C). 
Kuo’s (1952) result is shown as curve K. 

An arbitrary phase difference between heating 
and individual pressure change may be repre- 
sented by choosing an arbitrary finite value for 
6. For convenience we will choose 6 to be 2/2, 
representing maximum heating where in- 
dividual pressure change vanishes. For this 
case we may express I” by I, (1 +ib), and 
equation (20) becomes 

a Bb ja 
oO; = - 
CPS Poy EL ET) 


(+2 - ar) [aue(r-F2E=2)) - 


4 
À pète | 
ei = AU (21) 
For b= —0.7, which represents heating in the 


pressure troughs and cooling in the pressure 
ridges at the same rate referred to earlier, 
equation (21) yields curve C in Figure 1. 
This curve exhibits instability at short wave 
lengths qualitatively similar to the linearized 
three dimensional wave studied by Kuo. It 
also exhibits instability at long wavelengths, 
but here there is considerable difference from 
Kuo’s result. 

Heating in the pressure ridges with cooling 
in the troughs may be represented by setting 
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b equal to +0.7. This results in instability very 
much like that shown by curve C but displaced 
downward slightly and with no instability at 
long wavelengths. This is not shown in the 
figure. 

Heating distributions intermediate between 
those already discussed may be readily cal- 
culated for any arbitrary phase angle, 6. Their 
qualitative nature may be determined by 
superimposing two curves appropriately 
weighted, say curves B and C. 


Significance of Results 


In the following I should like to push 
interpretation somewhat beyond the valid 
limits of the preceding analysis; that is, to try 
to draw insights into the mechanism of the 
real atmosphere from the behavior of this 
greatly simplified and restricted model. There 
is no doubt that this is speculative and poten- 
tially dangerous, but often progress is made 
through initial insights which are largely 
intuitive. 

Let us imagine a series of moving sinusoidal 
disturbances and a distribution of heating 
which is fixed with respect to the earth. This 
condition might represent roughly winter- 
time heating and evaporation off the east 
coasts of North America and Asia. The results 
represented by curve B suggest that when 
the disturbance is situated so that the middle 
troposphere pressure trough is located to the 
west of the warm source (say, along the coast), 
waves should amplify rapidly. When the pres- 
sure trough is located above the warm source, 
rate of growth is then only slightly greater 
than with no heating; and when the northwest 
current lies above the warm source, rate of 
growth is much less and the disturbance may 
be neutrally stable. Finally, when the ridge is 
located above the warm source, rate of growth 
is similar to the rate with trough above the 
warm source. So one may imagine energy 
being fed into a train of moving waves in a 
series of pulses, each succeeding wave being 
amplified rapidly as the trough passes the 
coastline. 

Figure 1 indicates that dependence of in- 
stability on wavelength is strongly influenced 
by the distribution of heating. However, this 
may not be valid, for comparison of curves A 
and K shows that the two-layer model is 
seriously deficient in failing properly to describe 
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the dependence of instability on wavelength. 
Curve A indicates growth sharply limited to a 
narrow range of wavelengths, and it is easy to 
realize that this may result in serious errors in 
prediction. Comparison of curves A, C and K 
suggests that such errors in prediction arising 
from the short and long wave stability of the 
two-layer model might be reduced substan- 
tially by arbitrarily introducing heating in the 
troughs and cooling in the ridges. 

The two-layer model requires that maxi- 
mum @ occur at 500 mb for all wavelengths. 
This property violates the observed structure 
of baroclinic disturbances as well as the solution 
of the relevant differential equation (Kuo, 
1952); this we must recognize as payment for 
the simplicity of the method used here. If the 
restriction is removed (for a simple example, 
see FLEAGLE, 1957), short waves become un- 
stable though less unstable than waves of about 
5000 kilometers length. 

Long wave stability arises in equation (19) 
through the strong effect of the squared ß 
term at long wavelengths. One should not infer 
from this, however, that long waves are 
necessarily stable. The long wave instability 
revealed by curve C in Figure 1 reflects the 
appearance of a phase difference between the 
fields of vertical velocity and of individual 
temperature change at constant pressure. It is 
not difficult to demonstrate that if other 
phase differences were present, say between 
vertical and meridional velocity components, 
similar long wave instability may result. Phase 
differences of this sort are always observed in 
atmospheric disturbances and they must be 
expected in any but very much oversimplified 
theoretical treatments of baroclinic waves; 
there is, therefore, no convincing evidence for 
“stable wavelengths” in the real atmosphere, a 
point which I did not comprehend when the 
1957 paper was written. There is, however, a 
considerable difference in growth rate be- 
tween wavelengths of about 5000 kilometers 
and all others, and this is the most important 
point. 


Critical Observations 


(1) Replacement of the vertical derivative by 
finite differences violates the physics expressed 
by the governing differential equations; the 
maximum vertical velocity should not be 
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prescribed as occurring at an arbitrary surface. 
As pointed out earlier, this has serious conse- 
quences in introducing a false short wave 
cut-off to instability, but the wavelengths of 
greatest instability probably are not greatly 
affected. 

(2) Observations indicate a complicated 
spatial distribution of heating. The most 
important heating probably results from 
condensation in the middle troposphere, but 
this probably is not uniformly distributed be- 
tween 750 and 250 mb as assumed in the illus- 
trative calculations. Very intense heating may 
occur at and just above the earth’s surface. 
Heating is distributed horizontally by the 
geographical distribution of land, sea, ocean 
currents and surface characteristics, and not, as 
assumed here, entirely by the dimensions of the 
growing disturbances. 

(3) Linear theory becomes increasingly inac- 
curate as amplitude increases and is incapable of 
representing such important features as the 
development of fronts. 

These criticisms are sufficiently serious that 
one should be wary of drawing generalizations 
of wide applicability. I believe that the proper 
attitude toward the results discussed here is that 
they provide qualitative estimates of certain 
effects of heating by which one may judge 
what is and what is not to be expected from 
adiabatic theory. It also provides tentative 
relationships, such as the large heating effect 
when heating and vertical velocity are in phase 
or the effect of phase shift in introducing in- 
stability at all wavelengths, which it may be 
interesting to verify by independent means. 


The results indicate the following: 

(1) Heating of the magnitude which accom- 
panies condensation in the middle tropo- 
sphere may result in a significant change in 
growth rate of the most unstable waves. If 
such heating occurs in the ascending air east of 
the pressure trough, increase in rate of growth 
achieves its greatest value. If such heatin 
occurs in the descending air to the west of the 
trough, instability is strongly inhibited. 

(2) Heating which occurs in the pressure 
trough may result in only slight increase in 
instability at the wavelengths of greatest in- 
stability, but it has the effect of introducing into 
the two-layer model instability at all wave- 
lengths. 
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The Atmospheric Budget of Angular Momentum 


By G. B. TUCKER, Imperial College, London! 


(Manuscript received September 26, 1959) 


Abstract 


The mean meridional circulation derived from wind observations is combined with surface 
torque and mountain torque measurements to obtain a complete angular momentum budget 
for the northern Hemisphere between 15°N and 70°N for summer and winter. A previously 
proposed scheme of vertical flux of angular momentum involving downward flux everywhere 
except in the lower layers of the tropical easterlies is confirmed. The layer of zero vertical flux in 
the tropics is approximately at the level of maximum easterly wind. 


1. Introduction 


The mean meridional velocities at the surface 
over the oceans, observed by RIEHL and YEH 
(1950), have been used by PALMÉN and ALAKA 
(1952) to obtain an estimate of the various 
terms in the angular momentum budget of the 
atmosphere over a limited latitude range. The 
variation of mean meridional velocity with 
height had not been observed and therefore a 
hypothetical decrease to zero at 700 mb. with 
a compensating flow above was assumed. Re- 
cently however, the mean meridional circula- 
tion in summer and winter has been derived as 
a function of latitude and height from an analy- 
sis of observed winds in the area 160° W—O— 
40°E, 75°N—15°N (TUCKER, 1959). These 
observations have been combined with obser- 
vations of surface torque (PRIESTLEY IOSTa) and 
mountain torque (WHITE, 1949) to arrive at a 
fairly detailed and complete angular momentum 
budget for the Nothern hemisphere in winter 
and summer. 

The mean meridional (and zonal) velocities 
were previously derived (TUCKER, 1959) in an 
analysis based upon charts of the field of mean 
meridional motion in the Northern hemisphere 
drawn from as complete a network of surface 


1 Now with the Meteorological Office, London. 


and upper air radar and radio-wind stations as 
possible. These charts were constructed at vari- 
ous standard pressure levels from the surface 
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Fig. 1. Mean meridional velocity. Unit: m/sec, southerlies 
shaded. (After Tucker, 1959). 
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to 300 mb. Complete longitude cover was 
possible for the surface but not for the upper 
levels and the meridional velocities were correc- 
ted for the limited longitude field. The results 
are represented in Figure 1 which depicts the 
mean meridional velocity for summer (June, 
July and August) and winter (December, Janu- 
ary and February) as a function of latitude and 
height. A two-year period of observation was 
chosen for the computations, from October 
1949 to September 1951 inclusive. 


2. The budget of angular momentum 


(a) Theory 


Consider an element of volume dV and an 
element of the surface of this volume dS with 
do the projection of dS on the meridional plane. 
If 2 is the absolute angular momentum per unit 
mass, then following Widger (1949) the equa- 
tion of balance of angular momentum can be 
written as 


Sd SoQdV= [oQv,dS + [prdo+ [rr.dS (1) 
ot ÿ S G Ss 


where +, is the total zonal (eastward) frictional 
stress at the boundary, r the distance from the 
earth’s axis, v, the inward component of velo- 
city across the surface dS, and the integral is 
taken over a complete latitude belt extend- 
ing from the surface to the top of the atmos- 
phere. 

The term /prdo represents the torque due to 

oO 


pressure differences across mountain ranges 

where pdo does not vanish (WHITE, 1949). The 

term /rt.dS represents an exchange of angular 
5 


momentum due to frictional interaction with 
the earth’s surface. Both these terms are only 
effective at the earth’s surface. 
So@v„dS represents the convergence of an- 
3 


glar momentum advected across the boundary 
surface. If we consider a ring of the atmosphere 
with horizontal surfaces, H, bounded by la- 
titudes ®,, ®, and vertical surfaces, A, bounded 
by pressures pj, pa then writing Q=ur +wr?, 


peed de (Nero fornia ls 
Ss A ÿ i 


+ Here o/or’wdH |r: (2) 
H H 


pı 
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where u, v, w are components of motion in the 
x (east), y (north) and z (vertical) directions 
respectively. Squared brackets imply the differ- 
ence between the value of the expression inside 
the brackets at the two boundary surfaces 
indicated outside. 

We will consider the long term average of 
total angular momentum around a latitude belt 
to be constant, that is 


a) 
ER peodV=o (3) 


Neglecting the correlation between @ and u, v 
or w the equation of balance of angular 
momentum becomes 


J pr do at Ir, A ave 
o FI 
+ Dez +@ Jeter re 
A A ®, 


- peed ie Dieu Le Ry 
H H 


pı 


+ pee ae PT TT ee (iii) 
A 12: 

+ Heer eu ce (iv) 

0 (4) 


The bars in terms (ii) imply a mean value with 
regard to time and around a circle of latitude; 
these terms represent the effect of a mean 
meridional circulation. Terms involving  rep- 
resent the effect of the earth’s rotation (w- 
angular momentum); terms involving u rep- 
resent the effect of the motion of the atmos- 
phere relative to the earth (relative angular 
momentum). The prime implies a departure 
from the mean, and primed terms represent 
the effects of turbulent motion on all scales 
less than the mean meridional motion, this 
is called the eddy flux. 

The following sections will be devoted to the 
evaluation of the terms in equation 9 through- 
out the atmosphere using latitude intervals 
(®,—®,) of five degrees, and pressure inter- 
vals (p1—p2) of 200 mb. 


(b) Surface torques (terms (i), Equ. 4) 

In order to determine the sources and sinks 
of angular momentum at the bottom of the 
atmosphere the surface frictional torque must 
be combined with the torque due to pressure 
differences across mountains. 
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Table 1. Surface friction torque (S) and torque due to pressure difference across pour D: 
Also conversion table for 1 dyne cm? throughout a five degree latitude belt. Unit 10" & cm? sec”. 
Brackets indicate extrapolated values. 


woah. vba vl sets 2) 2 TS SE ES eee 


Latitude (°N) 70 65|65 60,60 55155— 50150 4545 40!40 35135 30[30—25|25—20|20—15 

SY) S000 (0) |(— 1-0)| — 3.0] — 6.0 9-0|—10.5 10.5 10.5| + 1.5) + 6.0| + 12.0 

WINTER ie ....[(+ 1.0)}| + 0.5 o | — 1.0| — 2.0]/—— 2.0 o |+ 1.5} + 2.0] + 2.0|(+ 1.0) 

S ....|(—0.5)} — 1.5 | — 1-5] — 3-0] — 3.0]— 1-5|.—1-5 o + 4.5) 7-51 2:5 

SUMMER { 94 ....[(+ 1.0)| + 0.5 | — 0.5] — 1-5] —2-0]— 2.0] —1.5]| — 0-5] + 0.5} + 2-0](+ 1.0) 
Conversion table 

I dyne = 2.1 3-0 4-1 5.3 6-5 Te]. 8.9] 10.3 ET.6) 4 712.3] 112-9 


(i) Surface frictional torque. The contribu- 
tion of the surface friction term f1rtxd Hsurface 
H 


within a latitude belt is 


TT a 
RATER 
0 di 


where c is the surface drag coefficient, R the 
earth’s radius and longitude. This term is not 
accurately known mainly because of uncertain- 
ties about the value of c. PRIESTLEY (1950) 
estimated a generation of 4.1 x 107° g cm? sec”? 
between o and 30°S using c=0.0013, while 
WIDGER (1949) using c=0.003 computed the 
larger value of 7.3 x 107° g cm? sec”? between 
10°N and 30°N for January 1946. The most 
reliable estimate is probably that of PRIESTLEY 
(1951a) using ocean surface wind data and 
assuming a value of 0.0013 for the drag coef- 
ficient. Priestley recognised that for various 
reasons (notably because the inclusion of land 
surfaces would require a higher value of drag 
coefficient) his value may be an underestimate 
by as much as so per cent. Priestley’s values are 
used in this work and are multiplied by the arbit- 
rary factor of 1.5. This is probably the largest 
source of error in the following analysis. The 
surface friction torque (S) for each five degrees 
latitude belt is given in Table 1; values have 
been rounded off to 0.5 x 1025 g cm? sec? 
A conversion table is includedin Table 1 relat- 
ing angular momentum flux for the five de- 
gree latitude belt to linear momentum flux 
per unit area. 

(ii) Mountain torque. The contribution of the 
term / prdo in equation 3 has been calculated by 

oO 


Wire (1949) who used normal surface pres- 
sure fields for each month. 


Using White’s data the torque due to moun- 
tains (M) has been determined for Winter and 
Summer and is included in Table 1; values have 
been rounded off to 0.5 x 1025 & cm? sec”?. 

A comparison of surface friction torque (S) 
and mountain torque reveals several interesting 
points: 

(a) S and M are generally of the same sign 
(i.e. in the mean, pressure is highest on the 
windward side of mountains). 
(b) The seasonal variation of S is larger than 
that of M; the two quantities are of the 
same order of magnitude in mid-latitudes in 
summer but elsewhere M is only about 1/5 
of S. 
(c) North of 65°N there appears to be an 
input of westerly angular momentum into the 
atmosphere due to the mountain effect (M 
positive). Widger (1949) has located a weak 
source of westerly angular momentum due to 
surface easterlies as far south as 60°N in Janu- 
ary 1946, but information is still scanty as to 
the magnitude, extent and persistence of 
polar easterlies. 
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Fig. 2. Sources and sinks of angular momentum at the 

earth’s surface, and total horizontal transport across lines 

of latitude necessary for balance. No transport is assumed 

to occur across 70°N. Arrows indicate the direction of 

transport of westerly angular momentum. Unit: 10% g 
Gis seca: 
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Table 2. Northward mass transport between 1000 mb and 700 mb. (unit 10% tons sec-!), 


Latitude °N | 60 


Palmen and Alaka. | | | | | 


O 


15 | 10 


| 


Fa 190 tes 67 


Inter ost — 2 E12 2178| 7 228 40| 246 8 49 92 IIo 
SUMMER). cance os —7| + 8 4+ 22| +15) + 30 Of 77 46 I 


Table 3. Upward mass transport through the 700 mb. surface (unit 10° tons sec-1). 


: 6 65 60 55 50 45 fe) (6) 2 
Latitude °N 4 35 3 5 
Se eg | 50 | ze] co le, | le 
Palmén and Alaka | | | | | 26 42 pl 7 | 56 2207 
NAN stern. | Se SS ey o O0 03 Où CEA 51 26 9 
SNS see Lees tO 8 I 27 33 Ch || 260 


(iii) Total surface torque. Values of S and M 
(Table 1) are combined in Figure 2 to give the 
total surface torque together with the horizon- 
tal transfer of angular momentum necessary for 
balance. In order to derive this horizontal trans- 
fer it was assumed that the source of angular 
momentum north of 70°N is negligible and 
therefore there is no horizontal transfer across 
this latitude. 


(c) The contribution of the mean meridional 
circulation (terms [ii], Equ 4) 


The mass transport across a given vertical 
surface on a complete latitude circle is given by 


M,=— fvdp 
pi 
aor ee (s) 


where 7 is the mean with regard to pressure 
of v. M¢ has been calculated between 1000 mb. 
and 700 mb. for each 5 degree latitude circle, 
these values given in Table 2 are compared 
with mass transports used by PALMEN and ALAKA 
(1952). These authors used the January surface 
meridional components computed by RıEHL 
and VEH (1950) and assumed them constant to 
goo mb. decreasing to zero at 700 mb. PALMEN 
and ALAKA used a 1010 mb. surface whereas 
this analysis is based on a 1000 mb. surface. 

The mass transport through a given horizon- 
tal pressure surface between two complete cir- 
cles of latitude is given by 


M, = 27R2?(sin 6, - sin $1) ew (6) 
Tellus XII (1960), 2 
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where w is the mean with regard to latitude of 
w. Macomb has been calculated for 5 degree 
latitude belts, these values are given in Table 3 
and compared with the values of Palmén and 
Alaka. 

There is reasonable agreement between the 
values of Palmén and Alaka for January 1946 
and the winter values except that their maxi- 
mum descent occurs between 15°N and 20°N 
whereas this analysis shows it to occur some 
10° further north. 

The flux of angular momentum due to the 
mean meridional circulation involves the com- 
putation of the four terms in (ii) of equation 4. 
Consider the block A (p,,®,), B (p;,®,), C 
(p1 Da), D (pa ®;) in Figure 3. Suffixes «,B,y,6 


latitude latitude 
2, Ai 
D | A 
as pressure 
Ful K ) . 
——— Le  — 
Fo) Fo) 
FT  — <<. 
FW) Fu 
DL Eee 1er = pressure 
C 5) ge) er: 2 
Fig. 3. Angular momentum flux across the boundaries of 
the block. 
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Fig. 4 and Fig. 5. Flux of relative angular momentum and © -angular momentum due to the mean meridional 
circulation. Arrows indicate the direction of flux. Relative angular momentum is shown above horizontal arrows 
and on the right of vertical arrows. w-angular momentum is shown below horizontal arrows and on the left of 


vertical arrows. Figures inside circles denote net convergence (shaded) or divergence within each block. Unit: 
95. Sicm® sec3, 
1075 g cm? sec 
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refer to the surfaces AB, BC, CD and DA 
respectively, the bar — denotes a mean value 
over the surface considered. 

The northward fluxes of relative angular 
momentum across AB, CD are F, (u) and F,(n) 
respectively, where 


F.,(u) = Mg, tz R cos di 
F, (u) =Mg, 1, R cos dy 
and 
F. (u) - F,(u) ES oruv sa fe in equation 4 
pi 


The northward fluxes of w-angular momen- 
tum across AB, CD are F;(w) and F, (a) 
respectively, where 

F.,(@) = Mg, oR? cos? ¢, 
F,(®) = My, oR? cos? dg 
Fa (w) - F,(@) = & no in equation 4 
A oy 
The upward fluxes of relative angular mo- 


mentum across BC, DA are Fg(u) and Fs (u) 
respectively, where 


Fg (u) = My, tip R cos (ES) 
Fs(u) = My, 45 R cos (a #) 


and 


Fg(u) — Fs (u) SIRES EEN in equation 4 


The upward fluxes of w -angular momentum 
across BC, DA are Fs (w) and F5 (w) respecti- 
vely, where 


139 
Fs(w) = M,,®R? cos? (er) 


Fs(w) = M,, oR? cos? (ei =) 
and 
Fg (@) — Fa (®) = onen Va in equation 4 


These terms have been computed for Winter 
and Summer and are presented in diagramatic 
form in Figures 4 and 5. The encircled figures 
give the convergence (shaded) or divergence of 
angular momentum in each block due to the 
mean meridional circulation. Since the total 
mass transfer across a latitude is zero the w- 
angular momentum terms cancel out when 
summed through the depth of the atmosphere. 
This is not true of relative angular momentum 
transfer which contributes to the total horizon- 
tal flux, the significance of this “circulation 
transport” has been a source of disagreement 
amongst meteorologists. The contribution of 
this relative angular momentum to the total 
horizontal transport has been calculated from 
Figures 4 and 5 and is given in Table 4 to a 
nominal accuracy of 0.1 x 102° g cm? sec-2, the 
ratio (in per cent) with the total horizontal 
transfer necessary for balance is also given. In- 
cluded in this table are similar values for January 
computed by Palmén and Alaka. 

Tucker (1959) described a method of ad- 
justing values of v to correct for taking a 200 
degree latitude arc instead of the complete circle. 
No similar correction can be applied to the 
values of u. In winter therefore the very strong 
zonal velocities near Japan are not sampled. To 
attempt to appreciate this the geostrophic zonal 


Table 4. Contribution of the mean meridional circulation (Fy) to the total horizontal flux of an- 
gular momentum (F,,) across a latitude circle. 


Latitude | | 65 | 60 | 55 | 50 | 45 | 40 | 35 | 30 | 25 | 20 | 15 
Py = CLP 0.1 0.5 0.8 II 2-4 3.1 1.0 2.31 + 5.5] + 6.5 
100 F,,/F,|— 20 20 20 8 5 Io 9 2 + 6 +18] + 36 
rte. Fy Using Petterssen’s 2.3 2.9 4.8 2.9| + 5:51 13.114 16.4 
100 Fy/Fr values of u II 13 14 27 + 14| + 42] + or 
1 

Fu Palmen and Alaka | | | | Bil, AIT 14 

100 FF Io 19| + 30 

Loge nee OR u En eS AE Se SS N ee a ee 
Ie, fo) + O.I O.I 0.8 0.8 o + 1.3] + 2.0] + 1.7} + 0.8] + 0.4 
Summer + +1750 IB ATO O + 20 4 II 7 fo) = of ie |) Bo hy 18| + 40 
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velocities for Winter computed for complete 
latitude circles by PETTERSSEN (1950) were com- 
bined with the values of v to assess the contri- 
bution of the mean meridional circulation to 
the horizontal flux of angular momentum. This 
composite model was used only for this cal- 
culation and the results are encorporated in 
Table 4. 

The results agree qualitively with those of 
Palmén and Alaka — the contribution of the 
mean meridional circulation is generally small 
but significant. In low latitudes the contribution 
is quite large and the composite model shows 
it to be predominantly important. In middle 
latitudes the indirect circulation and the increase 
of westerlies with height cause the contribution 
to be southward so that the northward hori- 
zontal eddy flux must be greater than required 
by surface stress measurements. 


(d) The horizontal eddy flux 


To complete the angular momentum budget 
the horizontal and vertical eddy fluxes must be 
included (terms (iii) and (iv) in equation 4). Not 
withstanding much work by PRIESTLEY (1949), 
STARR and WHITE (1952 a, b; 1954 a, b) and 
WIDGER (1949), the horizontal eddy fluxes are 
known only very approximately. Therefore the 
meridional circulation transport, Fu, was sub- 
stracted from the total transport, Fr, and the 
residue was attributed to the horizontal eddy 
flux. A percentage of this cross-latitude eddy 
transport Fr—Fm, was then allocated to each 
layer on the basis of the computations of STARR 
(1951) and Starr and Wuire (1952 a, b). At 
the time the computations for the present study 
were carried out Starr and White’s values were 
available only for 13°N and 31°N, consequently 
it was assumed that the relative importance of 
each level at 31°N applied to all higher latitudes. 
Subsequent work (Starr and WHITE 1954 b) 
shows this to be a reasonable assumption. The 
percentage of total horizontal eddy transfer allo- 
cated to each layer are given in Table 5, and 
the values are entered in the complete angular 
momentum budgets (Figures 6 and 7). 


(e) The complete angular momentum budget 


All the terms in equation 4 have now been 
computed except the vertical eddy flux (iv), 
this can be determined from the necessity for 
balance. The complete budgets for winter and 
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Table 5. Percentage of total horizontal eddy 
transfer allocated to each layer. 


Latitude °N T5 | 20)| 25 18011855 
Layer (mb.) 

o— 200 ..| 10 | 20 | 25 | 30 | 30 

200— 400 ..| 60 | 50 | 50 | 45 | 45 

Winter 4400— 600 ..| 20 | 20 | 20 | 20 | 20 

0600 7,800. 210.793 270 5 5 5 

800—I 000 ..| © fe) fe) fe) o 

6 — 200). | 01 1752| 20 1252025 

200— 400 ..| 70 | 55 | 55 | 50 | 50 

Summer? 400— 600 ..| 15 | 15 | 15 | 15 | 15 

600— 800 ..| Io | 10 | Io 5 5 

800—1 000 ..| 5 5 fo) 5 5 


summer are given in Figures 6 and 7, the en- 
circled figures in each block are the convergence 
or divergence due to the existence of a mean 
meridional circulation (the encircled values in 
Figures 4 and 5) and the vertical eddy fluxes 
have been computed to give balance. The ar- 
rows across the boundaries of each block and 
the associated figures indicate the direction 
and magnitude of the eddy fluxes. The fluxes 
through the earth’s surface were discussed in 
section 2 (b). 


3. The vertical eddy flux of angular mo- 
mentum 


Since it was realised that the hypothetical 
upward flux of angular momentum in the 
tropics presented by WIDGER (1949) must be in 
error (PALMÉN 1951), much attention has been 
focussed on the vertical transfer in the general 
circulation of the atmosphere. Recently SHEP- 
PARD (1953, 1954) has proposed a scheme where- 
by the upward transfer of westerly momentum 
near the surface in the Trade Wind region 
decreases rapidly with height and reverses to 
become a downward flux over the remainder 
of the troposphere — at least up to the level of 
maximum westerlies. He has incorporated a 
direct low-latitude circulation in his scheme to 
balance the divergence of vertical and me- 
ridional eddy flux at the upper (low latitude) 
levels and convergence at lower levels. In the 
region of surface westerlies Sheppard suggests 
that the downward flux is maintained with 
height up to the level of maximum westerlies 
where there is a convergence of horizontal eddy 
flux. He further suggests that no substantial 
meridional circulation is required for balance, 
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Fig. 6 and Fig. 7. Eddy flux of angular momentum. The surface (1 000 mb.) stress is taken from Fig. 2. Figures 

inside circles are the convergence and divergence of angular momentum due to the mean Meridional circulation 

(Figures 4 and 5), vertical eddy fluxes have been computed to give balance. Unit: 10° g cm? sec-?. A conversion 
table relating angular momentum to linear momentum is given beneath each figure. 


the convergence of horizontal eddy flux being 


approximately balanced by a divergence of 


vertical eddy flux. 


We are now in a position to compare the 
Tellus XII (1960), 2 


derived scheme of vertical momentum transfer 
with the empirical pattern proposed by Shep- 
pard for a uniform earth with no monsoon or 


topographic effects. 
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(a) Winter (Figure 6) 

North of 45°N the direction of the eddy flux 
is downwards throughout the atmosphere; this 
is strongest between 45°N and 50°N, exceeding 
4 x 1025 g cm? sec~? (0.6 dynes/cm?) throughout 
the column. The extension of the indirect circu- 
lation southwards to latitude 30°N (Figure 4A) 
counterbalances and exceeds the convergence of 
horizontal eddy flux causing an upward eddy 
flux in the upper levels. Also south of 25°N in 
the upper levels the convergence due to the 
direct circulation is insufficient to offset the 
divergence of horizontal eddy flux, again requir- 
ing an upward eddy flux. In the lower levels 
of low latitudes where the direct circulation is 
quite strong the resulting divergence of angular 
momentum far outweighs any convergence of 
horizontal eddy flux, consistent with a rapid 
changeover at about 800 mb. from an upward 
flux at the surface (in the easterlies) to a down- 


ward eddy flux above. 


Although there may be no a priori reason for 
demanding a downward eddy flux of westerly 
momentum in areas where westerlies increase 
or easterlies decrease with height in the mean 
(i. e. a positive coefficient of eddy viscosity), 
from physical considerations this seems more 
reasonable. The validity of Figure $ can there- 
fore be queried in the upper levels south of 
45°N. Since the upward eddy fluxes inferred 
here arise mainly from the existence of a strong 
indirect circulation it appears that the part of 
the Northern hemisphere not included in the 
observational analysis, in particular the strong 
winter jet in Southern Asia and near Japan, is 
strongly linked with the global meridional 
circulation. Thus it would seem that in winter 
a strong direct circulation may exist to the 
South of the Asian sub-tropical jet and this in 
conjunction with a strong horizontal tempera- 
ture gradient plays a major role in the mainte- 
nance of global westerlies. 
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(b) Summer (Figure 7) 


Apart from the upward flux at 600 mb., 
45—s0°N, which may be attributed to the 
accumulation of errors, the scheme of vertical 
angular momentum transfer is in good agree- 
ment with that proposed by Sheppard. The 
strongest downward transfer occurs at about 
30°N, reaching 14 x 1025 g cm? sec”? (1.3dynes/ 
cm?) at about soo mb. An estimate of the coef- 
ficient of eddy viscosity has been attempted, 
because of the crudity of the data this indicates 
only little more than an order of magnitude. A 
mean coefficient of eddy viscosity, K, was de- 
fined as 


where T,, is the mean shearing stress in the 
eastward direction, i.e. the vertical eddy flux of 
westerly momentum. The values computed for 
Summer are given in Table 6 and considering 
the crudity of the data they show a reasonable 
uniformity. Various estimates of the values 
of K for crosswind and downwind momentum 
transfer in the lowest few hundred metres of 
the atmosphere (RIEHL et al 1951, SHEPPARD 
and OMAR 1952, SHEPPARD 1954) place it at 
about 1x 105 cm?/sec with indications of an 
increase with height. An increase of K with 
height would, in the case of the Trades, 
provide downward fluxes aloft of the same 
order as upward fluxes near the surface where 
the velocity gradient with height (of opposite 
sign) is greater. Table 6 shows an increase of 
K from 800 mb. to 600 mb. and agrees fairly 
well with Riehl et al (1951) who found an 
eddy viscosity of 6 x 10° cm?/sec at about 700 
mb. in the E. Pacific above the trade inversion. 


(c) Level of zero vertical eddy flux in low 
latitudes 


Several studies of vertical eddy fluxes in the 
lowest layers of the Trades assume that the 


Table 6. The ‘mean‘ eddy viscocity in Summer. Unit 105 cm?/sec. 
ee ee ee ee En 
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turning point in the wind profile is the level at 
which 1,z (the shearing stress along the motion) 
is zero. This assumption may be tested on a 
climatological basis using the profiles of v = v( p) 
from Figure 1. 

The equation of mean zonal motion is ap- 
proximately 


The bar implies a mean in time and around a 


latitude circle, whence if we assume O 
der ek E 
ie A (8) 


Figures 6 and 7 show that in low latitudes a 
high value of 07,,/0y below 800 mb. is 
0.5 x 10% g cm? sec? (5° latitude)-1, or about 
1x10°° dynes cm? (100 m)-!, while a low 
Walton ob 7e es 30x10 gi cm? sect? 
(200 mb.)-!, or about 6x 10°? dynes cm? 
(100 m)-1. 97,,/dy can therefore be neglected 
in low latitudes at low levels and equation 7 
becomes 

de Bin 


(9) 


en 
Oz ge 


therefore 


(10) 


4 8 148 
Txz(p) = Tx2 (po) + = fv dp 
Po 


Txz(po). isthe surface torque. Values of 
the surface friction torque were taken from 
Table 1 and curves of v=v(p) from Figure 1 
were used to calculate the levels at which 
Txz(p)=0. Results are represented in Figure 8. 
Since the first level of observation above the 
surface was 850 mb. and the next 700 mb., the 
level of maximum easterlies cannot be accu- 
rately placed. However the line of zero zonal 
velocity has been drawn in Figure 8. RIEHL et 
al (1951) in a quasi-climatological study found 
the level of maximum north easterlies in the 
E. Pacific to be about 900 mb. at 20°N de- 
creasing northwards; SHEPPARD and OMAR 
(1952) in an analysis of winds over four tropical 
islands at about 20°N found a mean maximum 
of easterlies at about 0.5 km. for all wind speeds. 
These values agree fairly well with the level of 
zero Txz in Figure 8. 
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Fjg. 8. The level (dotted line) above which the upward 

eddy flux of westerly angulat momentum decreases to 

zero. The full line denotes the boundary between mean 
westerly and mean easterly flow. 


4. Conclusion 


The investigation of the angular momentum 
budget confirms the small contribution of the 
mean meridional circulation to the net hori- 
zontal flux of relative angular momentum 
except in low latitudes where it becomes compa- 
rable with the horizontal eddy flux. Perhaps 
the most important part played by the mean 
meridional circulation is in the horizontal and 
vertical flux of w-angular momentum. These 
large fluxes do not contribute to the net hori- 
zontal transfer but result in areas of convergence 
and divergence of advected angular momentum 
which balance corresponding areas of diver- 
gence and convergence of eddy flux in such a 
way that the direction of vertical eddy flux is 
downgradient. This confirms a scheme pro- 
posed by Suepparp and Omar (1952) and 
SHEPPARD (1954) namely that in low latitudes 
the divergence of angular momentum in the 
lower layers due to the mean meridional cir- 
culation is consistent with an upward flux in 
the surface easterlies which decreases rapidly to 
zero at the level of maximum easterlies and 
becomes downwards at higher levels. 
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Abstract 


The potential to kinetic energy conversion process in the lower stratosphere associated 
with the vertical exchange of warm and cold air is evaluated using adiabatically derived 
vertical velocities for the North American region for a five-day period. Preliminary results 
suggest the possibility that on the average the kinetic energy of stratospheric motions may 
not result from a conversion of potential energy within the stratosphere by this process. 
The further implication is that stratospheric motions are maintained by the motions in the 


adjacent layers of the atmosphere. 


1. Introduction 


A fundamental question about stratospheric 
motions refers to the manner of the maintenance 
of their kinetic energy. It seems clear that the 
kinetic energy of the motions of the troposphere 
is maintained by the conversion from available 
potential energy as described and verified by 
many investigators, most recently by LORENZ 
(1955) and WHITE and SALTZMAN (1956), and 
that this process is associated with the rising 
of warm and sinking of cold air. On the other 
hand, it has not been established that a similar 
process operates in the stratosphere. This 
preliminary study represents an attempt to 
evaluate from observations the nature of this 
potential to kinetic conversion process in the 
stratosphere. 


2. Procedure 
Following WHITE and SALTZMAN (1956) we 


may write the equations expressing the time 


1 Present address: The Travelers Insurance Company, 
Hartford, Conn. 
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rate of change of the kinetic energy of the hori- 
zontal wind and the total potential and internal 
energy of the entire mass of the atmosphere as 
follows: 


5, | kan - - [oxdm- [ Dam (x) 
M M M 


5 [(@+1)dm= [oxdn+ a) 


M M M 


where dm is the element of mass and the inte- 
gration is carried out over the entire mass of 
the atmosphere, k is the kinetic energy of the 


horizontal wind, 0-7 is the individual time 
rate of change of pressure, a is the specific 
volume, D is the rate of frictional dissi- 


pation of kinetic energy, ® is the geopotential, 
Tis the internal energy, and A is the net rate of 


heat addition. The appearance of the integral 
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fwadm with opposite signs in both equations 
M 


represents the process of potential to kinetic 
energy conversion in the atmosphere which is 
familiarly associated with the vertical exchange 
of warm and cold air. 

The critical problem in evaluating this inte- 
gral lies in the determination of w which is 
closely related to the field of vertical velocity. 
In this investigation © was evaluated by the 
adiabatic relation 


as se 
Ot 
ei DIT cz (3) 
Op Cp 


where T is the temperature, v is the horizontal 
wind vector on an isobaric surface, and C, is 
the specific heat at constant pressure. 

The finite difference form of equation (3) 
was evaluated from constant pressure charts at 
200, 100, 50 and 25 mb over the region of 
North America extending from 30° to 60°N 
and from 70° to 120°W. These charts had been 
carefully analyzed for other purposes by Mr. 
S. Muench of the Atmospheric Circulations 
Laboratory, Air Force Cambridge Research 
Center, and kindly provided us by Mr. W. 
Hering and Mr. S. Muench for this investigation. 
The available data consisted of seven 200, 100, 
and so mb and six 25 mb charts covering the 
period 28 January to 3 February, 1957, a period 
of intense stratospheric warming as descibed by 
Craic and HERING (1959). The wind velocities 
were evaluated geostrophically and the tempera- 
tures hydrostatically. Contour height values 
were abstracted at a grid system of 273 points 
spaced 2.5 degrees apart and the w’s were ob- 
tained at a grid system of 60 points spaced 5 
degrees apart. The time derivatives were ap- 
proximated by finite differences over 24 hour 
intervals. The computation yielded fields of the 
24-hour average values of w and « for the layers 
200—s50 mb and 100—25 mb which were taken 
to be representative of these fields at the 100 and 
so mb levels. 

Using the values of © and « thus derived it 
was possible to sample the integrand for this 
region of the hemisphere for the short period of 
time as indicated, and for the layers of the 
atmosphere at approximately the height of the 
100 and 50 mb levels. Average values of the 
integrand w« over space and time were compu- 


ROBERT M. WHITE AND GEORGE F. NOLAN 


ted and an analysis of the covariance of wa was 
performed in a manner entirely analogous to 
that of WnırtE and SALTZMAN (1956). Following 
their procedure we may write: 


{lo al} - {ol} {el} = {ol el} 
+ {lo + {lol} {al} (4) 


where the single, double and triple primes de- 
note deviations from east-west, north-south, and 
time averages respectively, the brackets indicate 
an average in the east-west direction, the braces 
an average in the north-south direction and the 
bars an average in time. 

The left side of equation (4) is a more repre- 
sentative measure of the integrand average than 
the term {[o«]} alone, since the second term on 
the left side must vanish when the average is 
taken over the entire hemisphere and will give 
rise to spurious non-zero values if not subtracted 
out because of the limited area treated. The first 
term on the right is a measure of the conversion 
process due to overturnings in the north-south 
direction, the second is a measure of the process 
due to overturnings in the east-west direction. 
The last term is associated with temporal pul- 
sations of the space average values of w and «. 


3. Results 


The values of @ obtained are reasonable. 
When transformed by means of the hydrostatic 
equation to equivalent vertical velocities, the 
standard deviation is 0.86 cm sec-! at 100 mb 
and 1.35 cm sec“? at so mb. A typical example 
of the distribution of the vertical velocity in 
relation to the temperature and contour height 
fields at so mb is shown in Fig. 1 for 29 Janu- 
ary, 1957. 

The numerical values of the potential to ki- 
netic energy conversion processes are shown in 


Table tr. 


Table 1. The rate of conversion between poten- 

tial and kinetic energy. A minus sign indicates a 

conversion from potential to kinetic energy. Units 
in ergs gm! sec-!. 


I ln diese 4 
{Lo «]} ele 
— {[eo]} {fag} | Slo)” tals} {Lo oT} {low} {La]}’” 
50 > su) — O.I | 7.5 | 0.1 
100 rab |= 0.5| — 0.3 | — 0.4 | 0.2 
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Bier. 


The first column of Table 1 indicates that at 
so mb the conversion process is opposite to that 
found in the troposphere being from kinetic to 
potential energy in this region of the stratos- 
phere. The sense of this conversion process was 
the same on each of the four individual days 
examined. At 100 mb the magnitude of the 
conversion term is not essentially different from 
zero and the six individual daily values on which 
the mean is based vary in sign. 

A more detailed analysis of this conversion 
process is given in columns 2, 3, and 4 of table 
1. At so mb the principal contribution to the 
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The vertical velocity, temperature, and contour height distribution at so mb for January 29, 1957, 03 
GMT. Vertical velocity units in cm sec}. 


conversion process is associated with over- 
turnings in the east-west direction as indicated 
in column 3. The contributions of the other 
two terms in columns 2 and 4 are negligible. 
At 100 mb the conversion rates associated with 
each of the individual terms are small. 

Since the 100 mb data are representative of 
the layer between 200 and 50 mb, this layer will 
be wholly contained in the stratosphere only 
at high latitudes and will be partially within 
the troposphere at low latitudes. It should be 
possible upon examination of the latitudinal 


variation of the term |w‘«’| to detect whether 


55 
LATITUDE 


“4e_ meon tropopouse 


ao 
- 
- 


Fig. 2. The latitudinal variation of the potential to kinetic 

energy conversion process due to east-west overturnings 

in units of ergs gm! sec~t. Also shown is the latitudinal 

variation of the mean tropopause height as a function of 
pressure. 


there is any systematic difference in the sense of 
the conversion process as one proceeds from 
troposphere to stratosphere. The latitudinal 
variation of this term at both 100 and so mb 
together with a plot of the mean tropopause 
height is shown in Fig. 2. 

It can be seen that the sign of the conversion 
process changes at approximately that latitude 
(43°N) where the mean tropopause crosses the 
100 mb level, being from kinetic to potential 
to the north and from potential to kinetic to the 
south of this latitude. This observational feature 
agrees with the concept of a reversal of the 
potential to kinetic energy conversion process 
from troposphere to stratosphere. 

If these observational findings are further 
substantiated by more extensive investigations 
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ona hemispherical basis and over longer periods 
of time then the implications seem far-reaching. 
If the kinetic energy of the stratospheric mo- 
tions is not maintained by conversion from 
potential energy within the stratosphere by the 
processes associated with the vertical exchange 
of warm and cold air then it must be main- 
tained by boundary interactions with the adja- 
cent layers above and below. Such processes 
involve the vertical transport of existing kinetic 
energy through the top and bottom boundaries 
of the layer as pointed out by STARR (1959) 
and other boundary processes associated with 
variations in the height of the bounding 
pressure surfaces. Should general verification 
of these concepts be obtained, one possible 
implication would be that stratospheric motions 
must to a large extent vary in response to 
tropospheric changes, and that the explanation 
for many observed characteristics of the stratos- 
phere, such as the size and motion of circulation 
systems and the seasonal and latitudinal distri- 
bution of ozone, may very well lie in a better 
understanding of the linkages between the 
stratosphere and troposphere. 


4. Critical Remarks 


As in all such limited investigations care must 
be exercised in generalization before confir- 
mation on the basis of more extensive data. 
This particular study suffers from the following 
deficiencies: 

a. The results are based on a sample of data 
from a small area of the hemisphere and for 
a very short period of time. 

b. The period of time studied was one of ab- 
normal high level temperature changes, and 
the results may not be typical of more nor- 
mal conditions. 

c. The vertical velocity computations on which 
the results fundamentally depend are based 
upon the adiabatic assumption. 
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On Critical Regimes and Horizontal Concentration of 


Momentum in Ocean Currents with Two-Layered System. 


By TAKASHI ICHIYE, Florida State University, Oceanographic Institute 


(Manuscript received August 20, 1959) 


Abstract 


The principles of minimum momentum transport and minimum energy which are commonly 
used in hydraulics of an open channel are applied to the current in a two-layered ocean, in which 
the geostrophic relation holds. For a flow of a finite width over the motionless lower layer 
these principles with the constancy of total volume transport predict the critical values of current 
width and upper layer thickness which will be attained by the flow after losing momentum or 
energy by friction. The two principles lead to critical regimes with almost similar order of 
magnitude. The critical width thus determined is less than 100 km and the upper layer disappears 
at the left edge of the flow in the critical regime, being in good agreement with the Gulf Stream 
and the Kuroshio. The analysis is also applied to a model of vertically continuous density 
stratification. Further the momentum transport of double jets often observed in actual currents 
is shown to be the minimum among multiple jets, next to a single jet corresponding to the 


critical regime. 


I. Introduction 


Many synoptic surveys done recently have 
confirmed the fact that ocean currents such as 
the Gulf Stream and the Kuroshio remain 
narrow and well defined not only along the 
coast but also as far offshore as at least 3,000 km 
from the continent. (FUGLISTER and WORTH- 
INGTON, 1951; ICHIYE, 1956). This is rather 
remarkable when it is recalled that frictional 
forces due to strong shear at the boundaries 
have a tendency to dissipate such narrow 
currents. In order to explain this fact, the late 
C. G. Rossgy developed a novel theory of the 
concentration of momentum (1951) based on 
the momentum principle which has long been 
used in hydraulics of an open channel. As a 
remarkable result of his theory, are obtained 
vertical density and velocity distribution of 
ocean currents which agree well with observed 
ones in the Gulf Stream in spite of the very 
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simple basic assumption. Thus he explained 
the strong current at the upper layer of the 
Gulf Stream as the result of vertical concen- 
tration of momentum transport. However, 
his promising theory had not yet the final 
goal, when his sudden death in 1957 prevented 
appearance of the second part of his paper 
which was expected to discuss the horizontal 
concentration of momentum. 

The author does not pretend to succeed to 
Rossby’s first paper, but does an application of 
the momentum principle and similar energy 
principle, to explain this horizontal concen- 
tration of the current as a critical regime of 
a jet-like current. In Rossby’s original paper 
and Craya’s subsequent paper (CRAYA, 1951), 
the density stratification of ocean currents is 
derived from the principle of minimum mo- 
mentum or energy flux, which gives the verti- 
cal velocity distribution similar to the observed 
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one in the Gulf Stream as a final stage of the 
current. The basic idea of their argument 
which has introduced these principles primarily 
used in hydraulics to oceanography is that the 
ocean current may reach the final stage or 
critical regime specified by these minimum 
principles after it loses momentum or energy 
through dissipative forces. The critical regime 
is determined from the condition of minimum 
momentum or energy with the assumption of 
constant volume transport, specifying a critical 
velocity in an open channel with a uniform 
density distribution. When the density of fluid 
is variable and thus the velocity is not uniform, 
the minimum condition leads to a critical 
regime which specifies a critical density distri- 
bution and corresponding current. 

Both authors considered gravity currents as 
a basic flow, in which a pressure gradient, if 
any, exists only in the direction of the flow. 
However, actual ocean currents are strongly 
geostrophic, that is, a pressure gradient across 
the flow has to balance the Coriolis’ force. 
Further, they are mostly confined in a rela- 
tively shallow layer of a finite width. Thus a 
two-layered system with a motionless lower 
layer well represents a model of the actual 
ocean in spite of its simplicity. When each layer 
has a different uniform density, the velocity of 
the upper layer is proportional to a gradient 
of thickness of the upper layer across the current. 
Then the total volume transport of the current 
with a finite width becomes proportional to 
the difference of the square of thickness at the 
sides. The energy or momentum transport is a 
function of the width of the current, thickness 
of the upper layer and its lateral gradient. 
Therefore these quantities, width and thickness 
instead of density distributions in the previous 
authors’ argument, are considered as variables 
in the variational problem for momentum or 
energy transport. Their critical values specified 
by the minimum condition of energy or 
momentum transport under the constancy of 
total volume transport correspond to a critical 
regime which will be attained by the current 
after its energy or momentum is dessipated by 
frictional forces. 

The critical regime determined from the 
minimal energy principle are generally different 
from those from the minimal momentum 
principle, as discussed by Craya (1951). 
However, critical values of width and thickness 
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by two principles are shown to be of a compar- 
able order of magnitude. Further a critical 
regime in ocean currents “may be only 
hypothetical terminus” which indicates only 
the sense of direction in future change of the 
state, as pointed out by Craya (1951). There- 
fore it is useless to discuss the difference of 
critical regimes determined by the two prin- 
ciples in details. 


2. Momentum transport of a two-layered 
ocean and minimal momentum condition 
in a uniform jet 


A two-layered system is considered, in 
which the current flows only in the upper layer 
with the depth D (Fig. 1). Density of the upper 
and lower layer is assumed to be constant and 
is put to o and Qo, respectively. The vertical 
coordinate z is positive downwards from the 
equilibrium sea surface and ¢ means the devia- 
tion of the surface from the equilibrium state. 
According to Craya (1951) the thrust P of 
the current per unit width is given by 


P=(D +6) go (00 — 0)/200 + (D+) ov? (x) 


for a flow which has a uniform velocity v at 
the upper layer with thickness D+ & and is 
motionless at the lower layer (in Craya’s 
original equation, flux v(D+L)=q is used 
instead of v). This equation can be derived 
directly from the definition of the thrust 
which is the variable part of the integral of the 
sum of pressure and momentum transport 
from the bottom to the surface. The pressure 


of upper and lower layer p and p, are given by 
p=gelz+¢), po=ge(D+ ¢) + gQ9(z- D) (2) 


Since there is no current in the lower layer, the 


Fig. 1. Block Diagrams of a Current in a Two-layered 
Ocean 
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pressure at the bottom py is constant and is 
given by putting z= Hin (2), where H is the 
total depth. 

The integration of pressure from the surface 
to the bottom gives 


D H 
frés+ [pode=*go(D +0 +¢0(D + 
= D 


I 
+) (H- D) += go9(H - D} = pil2g00 + 


+ (D + C)? go (00 - 0)/200 (3) 


in which the first term of the last expression is 
constant. Thus the thrust of the equation (1) is 
the sum of the momentum transport and the 
second term with (D+€)?, as shown in the 
equation (1). 

From the geostrophic relation the velocity 


of the upper layer is 
v = (0plox)|fo = (gif) (02/9x) 


in which x is taken across the current. The 
assumption that the lower layer is motionless 
gives dpjox=0 or IC/dx =(Ae/e)AD/dx, in 
which Ao=o,-0. 


Therefore the equation (1) becomes 
P=(D + 0% o¢4o/200+ 

+ (D +0) o(g4o/fe) (0D/0x) (4) 

Since & is of an order of one meter and D is 


of hundreds meters, we can put D + € & D. 
Further Jo ~ 2x 1073 and also we can take 


g =gAel oy ~ gAg/o. 


The equation (4) becomes 
P=oD[=¢D+(g/fP@D[ax| (6) 


When the current flows only in a belt 
extending from x=o to B as is the case in 
ocean currents, the total thrust is given by 


Me Be (6) 


which tends to decrease always through the 
loss of momentum. On the other hand the 
total volume transport V by the upper current 
remains constant and is given by 
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V= fovDdx Ziele) Di Do) a(7) 


by substituting the geostrophic relation for velo- 
city v, where Dy and D, are the values of D at 
x =o and B, respectively. The equation (7) 
indicates that V is independent of the width 
B. The critical regime for the momentum 
principle can be obtained by making M 
minimum under the condition of constant V. 

When the velocity v is constant within the 
belt, D becomes a linear function of x and is 


given by 

Dep (D; = Dy) x/B 
Then M becomes 

(re 
2 
- 0B [ir D ADS (D, D; + 2 ap") 
(8) 
in which AD= D, - D, and 
D?i= D} D;=2{ Vie 


In Rosssy’s or Craya’s model (1951) the 
minimum condition dP =o leads to the critical 
velocity (g’D)* under the constant transport. 
On the other hand, in the present minimum 
problem the variables are B and Dy (or D,). 
Thus the total thrust M is a function of these 
two variables and its minimum can be deter- 
mined from two conditions. 

Since D, is constant, we can introduce non- 
dimensional variables such as 


£=B/D, and 4,= DD, (f= 0,1) 
where 70 and 7, are related with 
mm 1 (9) 
Then the equation (8) is expressed by 


M= og Dié- 


Ant 2+ nom + 37 (m-no)&?} (10) 


in which 
na = olf? D, 


ise 


The total thrust M is a function of £ and 
no since n, and n, are related by (9). Thus the 
minimum condition for M becomes 


9MIdE=o or E2(H2+ 42+ No M1) = 37 (M1 — No) 
(11) 


and 


2MIdn = 0 or E{n+ni+ 470) = 37(M — No) 
(12) 


These two relations (11) and (12) with the 
condition (9) give the critical values of & and 
ni. If the width or the depth of the upper layer 
is specified, (12) or (11) separately leads to 
critical depth or width with (9). In general 
cases, (11) and (12) should be satisfied simul- 
taneously. Therefore we have 767, =0, which 
gives 79—0 with the condition (9). Thus in 
critical regimes the lower layer comes up to 
the surface at the left edge as often observed 
in the Gulf Stream and the Kuroshio (ICHIYE, 
1956). The critical width B, is given by 


B.=(3g'D,)?f-! or &=B.[D,= (3y)% 


The equation (11) is expressed by & as 


(13) 


EZ = £7 (9 + m) (1 + 27% + NoM) 


This gives &>E£, indicating that &, is the 
maximum of the critical £ under the condition 
that the depth of the current is constant. Simi- 
larly, when & is put into (12), it is concluded 
that the given width € must be smaller than 
&. in order that real critical value of n, or 7, 
may exist. Using numerical constants such as 


do=2 x 107%, feao™*(sec A) and 


D, = 800(m) (14) 
which correspond to the volume transport of 
64 x 10% M?/sec. about equal to that of the 
Gulf Stream, the critical width becomes 70 
km. This width almost agrees with that of the 
swiftest part of the Gulf Stream and the 
Kuroshio as revealed by recent surveys (VON 
ARX, 1952; ICHIYE, 1956). 


3. Variational problem for the minimum 
momentum condition 


In the preceding section a constant velocity 
is assumed in a belt of width B. Generally the 
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velocity is a function of x. Introducing non- 
dimensional variables y = x/D, and 7 =D/Ds, 
the minimal condition for M requires the 
integral 


I= i {n+ 2yn(dnldx)} dy (15) 


to be minimal by taking a suitable function 
n(x) and upper boundary & under the condi- 
tion (9). In (15) & and y are as defined in 
Section 2. 

In order to avoid the bothersome variability 
of €, the funtion 7 and an argument y are inter- 
changed. The equation (15) becomes 


T= fara 2 (0) #3 dn (16) 


in which y = dy/dn. The variation of the 

integral (16) is given by 
m In 
öI= [[Flndxdn + Fy ön (17) 

No 70 

(COURANT and HILBERT, 1953), in which 
F=F (n, x) is the integrand of the equation 
(16), Fy is a partial derivative of F with x’ 


and [F], is an Euler’s differential form. The 
condition that 6I =o leads to 


[Fh=o or Fy =1—2yn(y)* =e (18) 


Further the second term of the right side of 
(17) must vanish and this boundary condition 
gives 

c(önı — Oo) = (mom = 1)0m=0 (19) 
using the condition (9) and the relation (18). 
When the boundary value of 7, (or 7,) is 
specified, 07, =o and any value of c is pos- 
sible. However, when it is not so, c must 
vanish. Thus integration of (18) with n gives 


% = 2(2y} (mt - no) (20) 


The width of the current & is expressed by 


(21) 


The equation (20) indicates that the thickness 
D is represented by a quadratic function of the 
distance x. When the relation (9) is substituted, 
the critical width & becomes a function of ng 


Tellus XII (1960), 2 


E= 2 (2) (mt = Moi) 


OCEAN CURRENTS IN A TWO-LAYERED SYSTEM 


only. It reaches maximum 2 (2+)? for m9 =o, 
giving 
(22) 


This is equal to 113 km for the numerical 
constants of (14) and is larger than that of the 
previous section. 

In more restricted case of fixed n, the 
critical depth distribution is given by 


B, = 2 (29° D} f-2 


X= (oy) f {n/(n? — c)}dn (23) 


The critical width & is obtained by putting 
7 =M in this equation. It should be noted that & 
in this case reaches maximum for c= 10? 
because for any c which gives positive 72—c 
we have 
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When a function 7! (x) and 7! (x) are 
assumed to correspond to different integral 
constant cl and cl in the integral (23), & 
and EU are determined from (23) for cl and 
cl. Further assume that él > El The varia- 
tional problem which requires making the 
integral (15) minimal in an interval (o, él) 
has a solution 7! (x). On the other hand, the 
function 7!" (x) satisfies the Euler’s equation 
only in the partial interval (0, &"). Then the 
integral I given by 7! (x) is smaller than that 
of nl (x) in the range (0, él) (Courant and 
Hilbert, pp. 199—200). Therefore the function 
corresponding to c=n, is the real solution 
of the variational problem, because this gives 
the maximum value of é. 

The integral y (n) of (23) with c=? is 
expressed by the elliptic integral of the first and 
second kind, F (p, k) and E (y, k) (Byrp and 
FRIEDMAN, 1954, p. 82) as 


ey) =] = (en0) {F(, k) - 2E(p, 8) + 
+ 2 tan y(1 — k? sin? p)}} (24) 
in which 


ke == g=sin"1(9 — Nom)! 


Corresponding critical width is given by 
B, = (28 D,)*f-1J (ny) (25) 


in which (2 g’ D,)* f! is equal to 57 km for 
the constants of (14). 
Tellus XII (1960), 2 


3—002421 


153 
— X (2)/(2 198 


1 1.0 my 15 
I 


0,5 


0.5 

N 

? 

1,0 

Ik / / 
Re uud iS 


Fig. 2. Curves of x (n) (27) for different values of 7 

(in full lines) and corresponding & (2y)-# (in broken lines) 

obtained by the momentum principle. (Heavy lines are 

for the function J(7) and corresponding J (7). Fine lines are 
for the function of (20) and & (2y)~# of (21). 


The curves of J (7) are plotted for different 
values of n, in Fig. 2. The curve for n,= 0 
corresponds to the solution (20). The values of 
J (n) are proportional to the critical width 
B. as shown in (25) and the curve for 7,=0 
gives the maxiuium value 2 of J (7) as ex- 
pected. In this figure the curves of x (2 y)? 
and & (2y) ? from (zu) and (21) are also 
indicated. When y)=0, J (7) becomes equal 
to the curve of (20). When mp is finite, J (7) 
is larger than the corresponding function of 
(20) for the same 7, and thus the critical width 
& from (24) is larger than that of (21). The 
curves J (7) are interpreted as the thickness of 
the upper layer which is attained by the current 
after the loss of momentum under the condi- 
tion that the depth of edge of the current 7, 
(orn,) is maintained constant by some mecha- 
nism. On the other hand, all curves from (20) 
except one for 7)=0 do not represent the 
critical regime which gives the minimum 
momentum transport. Thus even if the current 
system eventually has a thickness distribution 
expressed by one of these curves, it will change 
further downstream through the loss of 
momentum and tend to approach the regime 
corresponding to the curve for 7) =0. 


4. Energy principle 


According to Craya (1951) momentum 
principle and energy principle lead to a 
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different critical regime for a stratified current. 
It is necessary to compare the results obtained 
from both principles in the present model. The 
energy flux per unit width e is defined by Craya 


2 


in which p and u are pressure and velocity in 
the general sense. This is given in the two- 
layered system by 


I ‚ I 
e=geAo(D + Ze + = av? ~ g eD + gv? (26) 


in which v is the velocity at the upper layer 
as defined in section 2. 

It is the total energy flux defined like M 
that should become minimum. In a model of 
uniform current within a belt x = o to B, 
computation similar to (5) gives 


B 
2E = 2 fedx = g'o(D, + D:)B + 
0 


+o(g' If} (4 D)?/B (27) 
By introducing non-dimensional parameters, 
E is expressed by 


E=-.og Die {no +m+ylm-n?E?} (28) 


The minimum conditions about two variables 
&and n, lead to 


dEldë=0 or nm+m)&=ylm-no)° (29) 


dElon =0 or (no +) = 27 (m1 — No)? (30) 


Therefore there is no critical regime which 
satisfies both of these conditions. However, the 
relation (29) gives the critical width &, for a 
fixed 7, such as 


& = (m — No) {7 /(m +790) $* 


This attains maximum y? for 7) = 0, giving 
an actual critical width B, = 40 (km) which 
is smaller than that derived from the momen- 
tum principle. Further the equation (30) gives 
the critical depth 7, or 7, under the constant 
value of £. 
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In case of variable v corresponding to 
Section 3 the total energy E is given by 


B 
E= go Dix +-e(g PJ (2D/0x) dx (31) 


0 


Introducing non-dimensional parameters, this 
quantity is proportional to 


L- J {an + y (dnldz)?} dy (32) 


Interchange of 7 and x leads to 


n 
dy dy ÿ. 
t= [Pon +y(#) dn (33) 
No 


The variational problem becomes 6L = 0 
under the condition 7,2—? = 1. The first 
integral of the Euler’s equation is given by 

2n — y (dy|dy)-? =d (34) 
with an integral constant d. 


In a case of unfixed boundary values 7, and 
M, d = o and we have 


(35) 


This equation indicates that the thickness D is 
also à quadratic function of x, but it gives a 
half of the critical width given by (20). 

In a case of fixed mo, Ë = x(n,) becomes 
maximum for d = 2 mp. Thus the criterion 
about the comparison of functions satisfying 
the Euler’s equation in different ranges of 
(0, €) leads to the critical depth distribution 


(36) 


x(n) = (2y)* (m? — no?) 


x (7) = (2y)* (m — no) 


corresponding to (24). This equation also gives 
the maximum critical width &, for No =O 
as in Section 3. The curves of y (2y)-* from 
(35) and (36) and corresponding & (2y)-# are 
shown in Fig. 3 in the same way as in Fir m, 
Comparing these results with those of Sec- 
tion 3, it is found that the critical widths 
derived by the energy principle are smaller 
than those derived by the momentum prin- 
ciple. However, critical regimes themselves 
are highly idealized state for the ocean currents 
as pointed out by Craya (1951). Therefore it 
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— X(2)/ (2 NE 


0 0,5 1.0 


> 


Fig. 3. Curves of y(n) (2y)~# (in full lines) and & (2y)-3 
(in broken lines) obtained by the energy principle. (Heavy 
lines are for the function of (36) and fine lines are for 


(35).) 


is almost impracticable to compare them with 
actual ocean currents. 


5. À system with continuous density 
distribution 

The minimum momentum principles can 
be applied to a two-layered ocean with con- 
tinuous density distribution. This is more 
realistic than the previous models. The most 
appropriate distribution in middle latitudes is 
given by STOMMEL (1954) as 


o =o,(const.); for D>z>o 


00 = Co — de exp {y(D -2z)}: for H>z>D 


(37) 


in which 60=0, —@; and » is a constant equal 
to 3 x 10° (cm1). Since H is assumed to be 


very large, we can put 


exp {»(D - H)}~o 
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The pressure at the lower layer p, is given by 
Polg = @3(D + 0) + @,(z— D) + 
+ de [exp p(D-2)} - ı]v"! = 
= Polg — 0 (H -2) + 6g exp {r(D - 2)}v"! 
(38) 
in which py is the bottom pressure assumed 
constant and is equal to 
Po=@,(D+C) + @(H- D) - dov-! 


The integral of pressure from the surface to 
the bottom is given by 


D H I fi 
Vire + f pode ee deD?/Q, + 5 Poleog an 


(39) 


+ gdov-2(1 —de/204) 


in which we put ¢ < D as in Section 2. 
Owing to the condition that there is no current 
at the bottom or that p, is constant, we have 


0;(9C/Ax) = de (9 D/Ax) 


The geostrophic relation with this equation 
gives the velocity at the upper and lower layer 
v and v, such as 


v= (a 1P)(ODIA), | to 
M = (g'1f)(2DI2x) exp (» (D -2)} 


in which g” = göo / es. By the approximation 
similar to the Section 2 we can put ¢” ~gdo/o, 
in the equation (39). Further, with the approxi- 
mation © + D ~ D, the integration of mo- 
mentum transport from the surface to the 
bottom gives 


D H 
Sov?dz + [ opvadz = 
6 D 


=D (D+) (ax) 


This relation is different from the one in the 
previous model only by the term (2»)-1 which 
is due to the finite velocity at the lower layer. 

Since the first term of the right hand side of 
(39) is variable, the thrust becomes the sum 
of this term and (41). The total thrust expressed 
by non-dimensional parameters is then propor- 
tional to 


€ 
Da) (dnldx}}dx (42) 


in which @ = (2v D.) 1 and y = g”/f°D.. 
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The total volume transport corresponding 
to (7) is given by 


In 


nn: (D, - D,) (Di + Do + 2/v)/f= 


= =o" Dif (43) 


or in non-dimensional form 


Ni — M6 + 4@ (M1 — No) = (44) 


The variation problem which makes 6I = o 
for the integral of (42) under the condition 
(44) is solved by the method similiar to the 
previous cases. In the case of non-fixed 7,, the 
solution of this variation problem is given by 


n 


a(2y')-* Su x} n !dn= 


=2{(n + a) (ne + a) + 
+ a log{(n + a)! - a3} -log {ro + 2)! - a} 
— log{(n + a)? + a} + log {no + a)? + a4}] (45) 


When n, or 7, is fixed, the same discussion 
about the integral constant of the Euler’s 
integral as in Section 3 and 4 leads to the solu- 
tion such as 


n 
x(2y)*=K(n) rs S(n +) (mn? - m) -*dn (46) 
The critical width £, can be obtained by taking 
n = M in this equation. The function K(7) is 
also expressed by elliptic integrals (Byrp and 
FRIEDMAN, 1954, p. 82) as 

(For 7 < &) 


K(n)= 2 (mo + «)!{F (p, k) - E(@, k) + 


+ tan p(1 — k® sin? p)t} (47) 
where 
k? = (a — 9) /(% + no), sin p= 
{7 — mo) ln + no) 
(For 0 > &) | 
K (1) = (no + «) (2Imo)H(R) tk? F(p, k) - 
-E(g,k)+tanp(ı —k® sin? pt} (48) 


where 
R=(n-ax)/2n0 (k’)P=1-k 
sin y = {(n — Ho) /(n + æ)}! 
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Fig. 4. Curves of y(n) (2y’)~# for different 7, (in fulllines) 
and & (2y’)-3 (in broken lines) for the system with con- 
tinuous density distribution. (Heavy lines are for K (nm) 
and K (n,) and fine lines are for the function of (45).) 


Fast {ieee 


The function K(n) and factor K(n,) of the 
critical width & are plotted in Fig. 4. It is seen 
that for the same value of and, K(n) is larger 
than J (7). However, the values of K (n,) are 
smaller than J (nı), except for the curves 
corresponding to 7, less than 0.1, since values 
of 7, determined from (44) are smaller than 
those from (9). As expected from the equation 
(46) the integral K(n) becomes infinite for 
Mo = ©. In this case there is no critical width 
and only the critical depth distribution can be 
obtained by taking suitable integral constant, 
say b, instead of 72 in the integrand of (46) for 
a given width &. This indetermination of 
critical width for a limiting case 7, = 0 is due 
to the finite velocity at the lower layer. When 
Mo ~ 0 the slope of the upper layer dn/dy 
near 7 = n, has to be proportional to 7, in the 
present model, whereas it is proportional to 
Ho! in the model of motionless lower layer. 
This requirement leads to an infinite value of 
x, since 


dy/dn ~ ny or x~ log % 


In Fig. 4 the curves expressed by the equation 
(45) are also shown. The values of K(n) and 
K(n) for the same n, and n are larger than 
those of the integral of (45). This integral also 
becomes infinite for 7, = 0. However, the 
critical depth distribution can be determined 
even for 7, = 0. Since y becomes infinite very 
slowly near 7, = 0 according to its asymptotic 
equation y~log 7, the critical width can be 
determined for 7) very close to o. 


Tellus XII (1960), 2 


OCEAN CURRENTS IN A TWO-LAYERED SYSTEM 


6. Momentum principle applied to a multiple 
current system 


It is assumed that the energy or momentum 
is dissipated from the current in the present 
model. However, there is an evidence that the 
ocean current also receives momentum or 
energy from the surrounding water (IcHIYE, 
1957). Therefore the results of the present 
theory cannot be applied directly to the actual 
condition of ocean currents. 

For instance, it is shown by FUGLISTER (1951, 
1956) that the Gulf Stream and the Kuroshio 
have a multiple structure consisting of a series 
of overlapping currents after it leaves the 
continent. In such a system the theory has to 
be applied to each branch of the multiple 
current. Then it is understood that each branch 
has a critical width and depth distribution 
predicted by the theory. 

On the other hand, it was observed often in 
the Kuroshio that the split of the current into 
two branches was caused by the intense 
dissipation of the momentum (or energy) of 
the current due to storms traveling along the 
current (ICHIYE, 1954, 1955). Further in some 
cases it was confirmed that the total energy of 
the current decreased after the split but its 
volume transport was almost unchanged 
(ICHIYE, 1958). Such examples suggest to us 
that the momentum or energy principle may 
be applied to the possibility of occurrence of 
the multiple structure. 

In the application of the momentum prin- 
ciple, the distribution of the thickness of the 


upper layer is assumed to be given by 
D=D,+AD-x/B+ D, sin (nax/B) (49) 


in which o<x<B. 
When n is an integer, this distribution yields 
the same volume transport for fixed values of 
D, and D owing to the equation (7). Yet the 
velocity distribution has a multiple structure 
for n larger than 2. Then the total momentum 
for different is computed and is compared 
with each other. The argument similar to the 
previous sections leads to the conclusion that 
the most probable mode of the thickness 
distribution is realized by the one which has 
the minimum value of momentum transport. 
The momentum transport for each mode of 
n can be computed by elementary calculation 
similar to (8). When M, means the value 
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corresponding to n = 0 or that given by (8), 
the value of M, for general n can be expressed 


by 


+ p(n) ie Ram — nn) + 
+(AD/B} E OME ii 


— 24r,(- 1) x |} 


(50) 
in which 
A=2¢'(AD)f-2B-*, 1,=D,/AD 
(= DA) 
p(n)=o (for even n) 


p(rn)=2 (for odd n) (51) 
Since only the relative magnitude of M,, is 
important, the values of (M,- M,)A-1 are 
shown for different n in the next table, in which 
the following numerical constants are used 


B=soolkm), D,= 0,1 (km), -AD=21 (km) 
D,=AD or =AD (52) 


DT 2 3 4 5 6 7 8 
I FOO, =) TO 0:2 BAO, 0722100 72E 


II—4.38 2.42 —I.00 1.78 0.00 2.04 1.08 2.70 

The first line corresponds to the case D, = 
AD and the second to D,= — AD. In the 
first case the minimum is attained for n =2 and 
M, for even n is smaller than for odd n. On the 
other hand, in the second case M, is minimum 
and M, for odd n is larger than for even n. The 
thickness for n = 1 in II gives the minimum of 
all M,, because this is very similar to that of 
the critical regime determined in Section 3. 
The modes n=2 in I and n=3 in II yielding 
the velocity distribution consisting of two jets 
give the minimum M, next to this. Therefore 
it is concluded that the double jets are the 
most probable state except a single jet corre- 
sponding to a critical regime. The result 
obtained from the energy principle is almost 
similar. Thus the split of the current into double 
jets can be attributed to the condition that they 
represent the second critical regime which the 
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current reaches after the loss of momentum or 
energy. 


7. Conclusion 


The principle of minimum momentum or 
energy applied to the geostrophic current in a 
two-layered ocean gives the critical regime 
which will be attained by the current as it loses 
momentum or energy but conserves volume 
transport. The critical width corresponding 
to this regime is less than 100 km. This effect 
explains qualitatively the fact that the Gulf 
Stream or the Kuroshio remains narrow and 
well-defined far from the coast. Further the 
critical depth distribution of the current deter- 
mined from the momentum principle predicts 
the vanishment of the depth at the left edge, 
being in good accordance with the observed 
structure of ocean currents. 

The multiple structure of ocean currents 
often observed off the coast is also principally 
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attributed to the second critical regime which 
is represented by double jets, because this state 
has the second minimum energy or momen- 
tum transport. However, the present theory 
indicates only the probable state of the current. 
Thus, for instance, it gives no information 
about the condition for occurrence of double 
jets, which may be treated from a different 
approach about the stability of the current. 
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Abstract 


Laboratory experiments which supplement those presented by STOMMEL, ARONS and FALLER 
(1958) are described. Extension of the analysis to include the effects of bottom friction shows 
that the characteristic spreading of zonal currents, apparent in the original experiments, is due 
to bottom friction combined with the radial depth variation. The experiments confirm the 
theoretical prediction of an intense “‘recirculation’”’ from the western boundary (similar to the 
“wind-spun vortex” of MUNK, 1950) when a localized interior source or sink is present. Appli- 
cation of the basic concepts of the earlier paper to the geometrically more complicated case with 
a partial radial barrier gave good agreement between theory and experiment, but the limiting 
effects of bottom friction were made evident in other examples. 

An eastern boundary current was found to occur in a situation geometrically similar to the 
intrusion of Mediterranean water into the Atlantic Ocean. Simple transformations illustrate 
the application of the equations developed in polar coordinates to the flow of a uniform 
fluid on a sphere, and examples of “recirculation”’ and the spread of zonal flow are worked 


out for a barotropic ocean of uniform depth. 


I. Introduction 


A recent paper by STOMMEL, ARONS, and 
FALLER (1958) hereafter referred to as (A) set 
forth qualitative and quantitative a priori pre- 
dictions of the circulation in a bounded basin in 
which the source of energy for the flow was 
the pressure head set up by the introduction 
and/or extraction of mass at various locations 
in the fluid. Specifically the analysis pertained 
to the circulation of a homogeneous layer of 
fluid in a pie-shaped basin rotating about the 
apex of the sector. Because of the rotational 
constraint the flow could be considered geo- 
strophic and two-dimensional except for certain 
boundary regions in which the viscous and/or 
inertia forces could produce appreciable de- 

1This work was supported by a contract with the Office 


of Naval Research and is Contribution No. 1069 from the 
Woods Hole Oceanographic Institution. 
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partures from geostrophic flow. The predictions 
of the simple theory were borne out by subse- 
quent laboratory experiments. 

This paper presents the results of some ad- 
ditional related experiments and expands upon 
details of the flow patterns of the earlier work. 
As in (A) the experiments are restricted to 
depths of fluid and rotation rates such that the 
flow is primarily geostrophic, but in addition 
to departures from geostrophy near the side 
boundaries the specific effects of the bottom 
frictional boundary layer are considered. It is 
shown that bottom friction is a most important 
consideration for situations where there are no 
unbroken radial barriers, and indeed the effects 
of bottom friction may be observed and are of 
some significance in all of the experiments. 

The analysis described in (A) is extended here 


to somewhat more complicated regimes of ge- 
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ophysical interest, namely the circulation with 
a partial radial barrier and the “recirculation” 
associated with localized interior sources or 
sinks of fluid. 

An unanticipated development was the oc- 
currence of eastern boundary currents in certain 
circumstances. It was asserted in (A) that only 
western boundary currents would occur and 
this assertion seemed to be a justified one on 
the basis of the original experiments. It now 
appears that eastern boundary currents of in- 
ertial character are possible when fluid is intro- 
duced through a slot in the eastern boundary 
in such a manner that viscosity may generate 
negative absolute vorticity. 

Finally the theory for the flat rotating tank is 
extended to a spherical shell of fluid. The 
recirculation phenomenon and the spread of 
zonal flow are computed as examples of what 
might be expected for barotropic flow in an 
ocean basin of constant depth. 


2. Equations for the Interior Circulation 


The polar equations of motion and the conti- 
nuity equation for a uniform fluid with an 
arbitrary distribution of sources (or sinks) are 
(DRYDEN, MURNAGHAN, and BATEMAN, 1956): 


dv, vo I Op 
Hey 2fve= pee 
. (oa + Lior, iv, LI Pv, 2 dve 
4 02 r OP) 1? 72998 7290 (1) 
dve v,vo I op 
yale + 2Qv,= "pr 707 
PVO LL Wen Ve) 10" Ve: > Ov, 
v - = — + z+ (2) 
a r or 1? 1? 990 r90 
dw I Op 
eto 
Kr Ow # Ow in I Ow : I. d?w 
an VO At) a il (3) 
OW SL Of DT 
Dei Tee (4) 


where the r direction is radially outward, © is 
in direction of rotation (counterclockwise) and 
z is vertically upward. P is a source of fluid 
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per unit of volume. Dimensional arguments 
may be used to show that for: ı) small flow 
rates, 2) small horizontal shear, 3) a steady state, 
and 4) large radius the equations reduce to those 
terms which are underlined and which were 
the equations used by Ekman (1905) in his 
studies of wind driven ocean currents. Again 
from dimensional arguments or from the so- 
lutions given below one may see that the depth 
of a viscous boundary layer at a rigid horizontal 
surface should be characterized by the Ekman 
v 
depth De = IT (3) 
aries (z>Dpg) we may omit the vertical de- 
rivatives in 1) and 2) to obtain the geostrophic 
equations 


4 
, so that far from the bound- 


expt Oa NE 
Eu era or (s) 
CP EL: 
st 2Ror 00 (6) 
de _, | 
a 96 Tara 7) 


Therefore for a total depth of fluid h> Dr with 
a free surface the following lower and upper 
boundary conditions apply: 


atz=0: 1, =Ve, =Wog=0 


0 
atz=h: ve, =v 
vu 


„9 ri v oh p oh . 
Wh = 
or r99 at ( ) 


oh R ERS S 
where Fri is the depth variation associated with 
a steady rise or fall of the free surface. The depth 
variation h’ due to the steady geostrophic 


currents need not be considered in the ex- 
nz 


pression for 1, since VV h=0. 

Since the boundary layer circulation de- 
creases exponentially with height, we may 
safely take the upper limit of integration in 
exponential expressions as h=co and so obtain 
as solutions to 1) and 2): 


a4 
= z > 
v,=u-e Pe[ucon —+vsinn— 
( DE DE (9) 
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rs 
DS : FA = 
vo=v-e (usin 7 — -v cos a— |}. (10) 
DE DE 


By substitution of these expressions into 
equation 4) we find: 


Ow -n . z {x du rôrv 
De ee eee es D: 
ee E 29 + =) ur, 


dz Dg 

Vertical integration of (11) from z=o to 
z=h(co) together with the boundary con- 
ditions (8) yields the following equation for 
the geostrophic flow: 


2h v ah oh De du Im MER 
or r 00 A 2x \r90 rr) (12) 


h 
where Q=/ Pdz is the source of fluid per unit 


O 
area. Equations (7) and (12) are sufficient to 


solve for the geostrophic flow given Q, Zr 


oh oh me 
30° u and suitable lateral boundary conditions. 
Introduction of a stream function defined by 


DUO. Mao: 
JO gg MO (7) and (12) and con- 


sidering only radial depth variations yields: 


u 


oh ic ow DE 2 oh 
PON Sana Gi (13) 


Equation (13) is of the same form as that derived 
by STOMMEL (1948) with the following differ- 
ences: 1) the radial depth variation enters in 
place of the latitudinal variation of the Coriolis 
parameter, 2) the effects of bottom friction have 
been derived more completely rather than to 
assume a linear friction law, and 3) no specific 
functional form for the source term has been 
considered here, whereas STOMMEL assumed a 
specific wind stress distribution in order to 
obtain an analytical solution. The above equa- 
tions will now be specialized for specific 
examples of interest. 


3. Transport across an Annular Ring 


As a first test of the adequacy of the Ekman 
theory as developed above a circular channel 
with inner radius 50 cm and outer radius 100 cm 
was set up on the same rotating apparatus as 
had been used in the earlier experiments. Pro- 
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Fig. 1. A schematic diagram of a source-sink experiment 

in an annular ring with the source continuously distrib- 

uted around the outer rim and the sink at the inner rim. 

For relatively slow flow rates the radial mass transport 
occurs entirely in an Ekman layer. 


vision was made for a source of mass con- 
tinuously distributed around the outer rim and 
a sink at the inner wall as indicated in Figure 
1. In this manner a steady mass transport (So) 
across the channel was forced, and since it was 
anticipated that this transport would take place 
entirely in an Ekman layer, the source and sink 
were arranged to feed the bottom boundary 
layer directly. If we assume circular symmetry 
oh oh 
(u=o) and take Q= Amy 
ately obtain from (12) the solution rv = C where 
the constant C must be determined by the 
condition on the total mass transport, 


= O0 we immedi- 


x h 
ro So: (14) 


Substitution of equation (9) into (14) gives 


directly: 


m=—. (15) 


The accuracy of this expression, which 
implies a zonal flow inversely proportional to 
radius, was examined experimentally with a 
source strength S,=8.67 cm?/sec. The follow- 
ing values of the zonal geostrophic flow were 
obtained from successive photographs showing 
the displacements of patches of dye at various 
radii: 


Table 1. 

Radius | v | Sp 
62.5 cm .272 cm/sec 7.35 cm°/sec 
74 +232 7-31 
87 .194 7:45 


£62 
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Fig. 2. Reproduced from Stommel, Arons and Faller (1958) to illustrate the spread 
of zonal flow from the source S, as it moved westward and the convergence of the 
flow as it moved eastward toward the sink S;. 


Values of transport predicted from the zonal 
flow (Sp) averaged 15 % below the imposed 
transport S)=8.67 cmÿ/sec. 

Estimates of the neglected terms of the 
equations of motion show these to be three 
orders of magnitude lower than those of the 
Ekman equations, (2) and (3), which fact is 
indicated directly by the Rossby number, 


R,= _ 2.004. Possible explanations of the dis- 


Qr 

crepancy are: 1) that the bottom boundary 
surface which was a painted metal surface had 
some roughness elements which tended to 
introduce irregular flow and an eddy viscosity 
above the molecular value »=.o1 at 20°C; 2) 
that a free surface stress existed because of some 
properties of the surface film. Other possibilities 
such as a wind stress or a thermal circulation 
were eliminated by a sealed glass cover and 
insulation of the tank. A systematic study of the 
effects of boundary roughness is planned for 
the near future, but it is believed that a surface 
film could easily have caused the observed 
discrepancy between theory and experiment. 
Several cases have been observed where con- 
taminants have caused the surface to act as a 
stationary membrane even in the presence of 
active (1 cm/sec) circulations within the fluid. 
Other experiments have shown that flow at the 
free surface is always somewhat slower than the 
interior flow. In the present example the Ek- 
man transport produced by a membrane would 
be radially inward and would augment the 
transport due to the bottom boundary layer 
with the result that the required values of v 
would be less than with the bottom boundary 
layer alone. 

More extensive experiments of a similar 
nature have recently been performed by Arons 
and Bryan at the Institute of Meteorology in 


Stockholm (personal communication) and by 
GREEN and Arons (Green, 1959) at Amherst 
College. 


4. The Spread of Zonal Flow from a Point 
Source due to Bottom Friction. 


Some of the experiments illustrated in (A) 
and subsequent work have shown a rather 
peculiar spreading of westward flowing zonal 
currents. Conversely, currents directed east- 
ward toward a point sink seem to broaden near 
the western boundary where they originate and 
then converge smoothly toward the sink. 
Figure 2 reproduced from (A) illustrates the 
spread of the westward flowing current from 
the source at the eastern wall near the apex and 
an orderly convergence of the eastward flow 
toward the sink. In Figure 6b the stream from 
the western boundary to the tip of the partial 
barrier definitely indicated radial spread near 
the western boundary and convergence as it 
moved eastward. Other experiments have 
shown this phenomenon more strikingly. 

We consider the flow from a point source in 
the eastern wall toward the western boundary 
where we regard the western boundary current 
as a sink. The circulation in the immediate 
vicinity of a source is exceedingly complicated 
and the analytical solution to be obtained can 
be valid only some distance from the source 
itself. Equations (7) and (12) are specialized by 
letting 
oh = dh oh 


nn = 
dt dO or oe 


Bsus Dis» Dia 
dy 109’ dx Or 


Q 


These conditions represent 1) no interior sources 
or sinks of fluid, 2) a constant radial depth 
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variation, and 3) conversion to rectangular 
coordinates. Then after cross-differentiation to 
eliminate # we obtain: 


dh dv = (= 22v 


+53) =0 + (19) 


2? ; F : P 

The term i will be omitted since it may be 
shown in terms of the resultant solution that 
dv ; bl dv 
a is comparable to DE 
The equation to be considered is therefore: 

dv a 2 dh dv - 
0x? Dr 0x dy © (20) 


ax dy 27 


dy? 


only near the source. 


subject to the integral condition that the trans- 
port be constant at any value of y, 


+00 h 
tif avdzde=Sy: 


=60. 0 


(21) 


In this integration we take h as independent of 
x as a first approximation. 
Equation (20) is similar to the one-dimen- 
2 


. ; Or 
sional heat flow equation ae K 


K is the thermometric conductivity. Corre- 
spondingly, the integral condition on mass 
transport (21) is equivalent to supplying a given 
amount of heat in a pulse at t=0 and x =o an 

thereafter allowing the heat to spread along the 
rod toward x= +00 and x= — ov. By analogy, 
the solution to (20) and (21) is (Carlslaw and 


Jaeger, 1947): 


= owhere 


RTE és) 


7 = han)! (24y} a 
DE . | 
where A= - Fo Inasmuch as À is negative, 
2H — 
Ox 


this solution has meaning only for negative 
values of y although v may be positive or nega- 
tive corresponding to a sink or a source re- 
spectively at x=0, y=o. Equation (22) re- 
presents a Gaussian distribution of v in x for any 
distance y from the source and with a standard 
deviation o, = (2Ay)?. 

Figure 3 illustrates the derived distribution of 
zonal and radial flow, of Ekman flow, and 
the vertical motions for a case with A=—1. 


Tellus XII (1960), 2 


163 


Fig. 3. The theoretical circulation associated with the 
spread of a westward flowing current from a source at 
x =0, y =0 for the value of the parameter À = —1. 


For negative values of x there exists divergence 
in the Ekman layer and a loss of mass to the 
fluid directly above, so that in accordance with 
the depth variation the interior geostrophic 
current must have a component toward negative 
x. The opposite is true for positive values of 
x, and as a result there is a systematic spreading 
of the flow. Along the axis of the stream there 
is clearly a net transport of mass toward positive 
x by the Ekman flow. Qualitatively this mass 
is absorbed on the left hand side of the current 
which is spreading to greater depths and there- 
fore progressively transporting a greater fraction 
of the mass than the right hand side which is 
spreading to shallower depths. However, this 
balance is not quantitatively correct because 
of the approximation h=const. in equation 
(21) which led to a symmetrical solution. In- 
2 


À Vv. N N 
spection of — in comparison with — indicates 
dy? Ox? 


that inclusion of this term would tend to retard 
the divergence of the stream although the effect 
is already small 10 cm from the source for 

= — I(an average value for the experiments). 

Quantitative experimental data to support 
these conclusions are not yet available, but we 
may estimate the spread of the westward flow- 
ing current in Figure 2 and compare this with 
the theory. The equation of the free surface is: 


h=ho(1 +15) (23) 
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Fig. 4. Photographs at 20, 80 and 136 minutes after the introduction of dye showing 
the existence of both eastern and western boundary currents for the case where the 
source water was injected through a slot in the eastern wall. 


(2242 : 
where |= Er and for the experimental values 


2g! 
hg=s cmyal= 1, arrootcm ro nsee and 
MoS so 
y=.01 cm? sec-t it is found that A= - —. 
r 


The radius of the source in Figure 2 is approxi- 
mately r=20 cm which gives A= - 2.5 and 
0x =(sy)#. If we consider that most of the 
dyed source water will be contained within 
the region +2 0, the rimward limit of the dye 
outline is in quite good agreement for various 
distances y from the eastern wall. Dye fronts 
are not adequate for quantitative comparison 
for several reasons including the fact that the 
cross stream Ekman flow causes the dye to be 
systematically displaced toward the left of 
the axis of the flow. 


5. Eastern Boundary Currents 


In (A) it was asserted as one of the essential 
elements which defined the regime to be con- 
sidered that departures from geostrophic flow 
were to be permitted only at western bounda- 
ries where narrow and rapidly flowing currents 
could occur. The function of these currents was 
to preserve continuity of mass in the basin as a 
whole by providing means for radial transport 
but the dynamics of the boundary currents 
were not explicitly considered. In nearly all 
experiments the above assertion was found to 
be a justified one, but in one particular situation 
a pronounced eastern boundary current was 
observed. 

Figure 4 is a sequence of photographs at 20, 
80, and 136 minutes after dye was introduced 
into the source water, there being a source (Sp) 
through a slot in the eastern boundary near the 
rim and a sink (S;) at the apex of the sector. 


As in the previous experiments, the radius was 
a=100 cm, rotation rate Q=1.05 sec}, the 


mean depth h=8 cm, and the ratio of the depth 
at the rim to that at the center h,/h, was 2. 
The volume rate of flow was S,= 110 cm?/min. 
As may be noted from the photographs there 
existed an obvious eastern boundary current 
transporting perhaps 10 % Sy, and repeated 
experiments showed this phenomenon to be 
reproducible. The hypothesis is advanced that 
a small portion of the fluid received negative 
absolute vorticity by means of the viscous forces 
at the edge of the source slot, and that this por-- 
tion of the fluid could take part in an eastern 
boundary current. 

Several theoretical models of boundary cur- 
rent circulations in the ocean have been sum- 
marized by STOMMEL (1958). It is clear from the 
various studies that for linear viscous models 
significant mass transports may occur only at 
western boundaries where it is possible to have 
a balance of changes of relative vorticity (£) 
due to latitudinal (or depth) variation with 
those caused by lateral or vertical viscous stresses. 
In addition, MORGAN (1956) has shown for the 
simplest non-linear models that only western 
boundary currents can be of significance. How- 
ever, inasmuch as the possibility of negative 
absolute vorticity has not been considered and 
since complete non-linear solutions including 
viscosity have not been available for examina- 
tion, there is no reason to rule out all possi- 
bilities of inertio-viscous eastern boundary 
currents. 

The simplest physical arguments suggest that 
if the fluid next to the eastern boundary had 
negative absolute vorticity (¢ + 22<o) it would 
again be possible to balance vorticity changes 
due to the depth variation with those due to 
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Fig. 5. Photographs at 10, 20, and 40 minutes illustrating the absence of an eastern 
boundary current when an interior source was used (compare with Figure 4) and 
the existence of strong “recirculation”? from the western boundary. 


lateral viscous stresses. However, in such an 
event many non-linear terms would become 
important and the simple arguments would be 
invalid. It is suggestive that for the experiment 
of Figure 4 6 % of the source fluid would enter 
with 8 +2 <o given a slot width of 1/4 in., the 
flow rate of 110 cm/sec, and assuming a para- 
bolic velocity profile for the flow through the 
slot. More critical experiments with various slot 
widths and flow rates are suggested as means of 
varying the percentage of fluid entering with 
negative absolute vorticity. 

As a first test of the above hypothesis, the 
experiment illustrated in Figure 5 was con- 
ducted. Conditions were identical to those of 
Figure 4 except that the source water was 
introduced through a tube just away from the 
boundary so that there would be no opportunity 
for the fluid to obtain negative absolute vorticity. 
As may be seen, no eastern boundary current 
occurred. 


6. Recirculation about an Interior Source 


or Sink 


Introduction of the source water in the 
manner shown in Figure 5 eliminated the 
eastern boundary current but caused other 
interesting circulations. It is shown below that 
with an interior point source (distributed over 
a small area by turbulent mixing) an ampli- 
fication of flow toward the western boundary 
is required. The sequence of Figure 5 shows a 
return flow from the western boundary current 
to the source, a circulation which is required 
by continuity in the vicinity of the source and 
the geostrophic constraints in the system with 
variable depth. To analyze this phenomenon 
we now consider the circulation about an 
interior source in the absence of bottom friction. 


Tellus XII (1960), 2 


oh 


From equation (12) setting a Bere o and 
y = o we obtain: 


oh 
MMS —-rQ=0 


(24) 


and by elimination of u with equation (7) we 


find: 


2h 

dv dr 1 2(rQ) 

D CN le 
or or 


Integration of (25) from longitude ©, to the 
eastern wall © = Oz such that the entire source 
is included in this region gives: 


Ph 


or rd 
aha (26) 
or or 


where we have taken ve,=0 (no zonal flow at 
the eastern boundary) and have written Qs = 
0, 
= i QrdO, the zonally integrated source per 
C7 
unit of radius. Using the equation for the free 
surface (23) we obtain: 


ae 2, (Qs 
SET AE OLN ale de 


Equation (27) shows that there should, in 
general, be both positive and negative values 


(27) 
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for vo, and that inflection points ve,=0 occur 
where (=) =o. For a smoothly varying 


source with a single maximum at r, we define 
the radius r (ve,=0) =r’ and write: 


fhodr+ fhvdr= fOsdr=S,. (28) 
0 Ta 0 


m nt 
2) 


m 
T(+) 


The integrals on the left of (28) represents 
transports toward and away from the source. 


We define TIC) 


5 8 the “recirculation” and find: 
0 


EEE“ 


Interpretation of (29) is facilitated by con- 

sidering a specific functional form, say the 
, dr ee Os oF _ (r= 10)" 

Gaussian distribution — =; —e  %* so 

r (2m)? reo 

that r’=r,, the radius of maximum source 

strength, and where o is the radial spread of 

the source. Therefore for a source with small 
radial spread we may write approximately: 


POSTE aa OR 
So ~( 5) {4()- 


a? i To 
A : (270)* 20° 


We see that for small o the recirculation is 
inversely proportional to o. In any physical 
situation a lower limit on the practical size of 
o probably occurs because of mixing associated 
with the source and with the increased circu- 
lation. 

From Figure 5 we estimate r,=85 cm and 
o=5 cm, and using the values a= 100 cm and 


T(+) 


1= 1, we obtain a theoretical recirculation —— 


(30) 


90 
=7.5. This surprisingly large value was checked 
by evaluations from Figure 5. The total volume 
of dyed fluid which had been introduced by 
the time of the first picture was 1100 cm? 
(110 cm*/min for 10 min). The volume of dyed 
water as evaluated from the photograph was 
10,440 cm? or 9.5 Sg At, which indicates that 
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the source water was diluted by mixing at the 
source of 8.5 parts of recirculated clear water. 
This result is in good agreement with the theo- 
retical recirculation considering the arbitrary 
specification of the form and width of the 
source region. Confirmation of these figures 
was obtained from the second photograph of 
Figure 5 which shows a volume of dyed fluid 
10 times the total source input at 20 minutes. 
Comparison of Figures 4 and 5 shows the very 
large differences in the speeds of the zonal flow 
without recirculation (Figure 4) and with the 
additional flow produced by recirculation 
(Figure 5). 

One may well inquire why the phenomenon 
of recirculation did not occur when the source 
water was injected through the slot in the wall. 
Formally the difference lies in the assumption 
following equation (26) where we have taken 
for the case of an interior source veg =O, i.e., no 
geostrophic flow through the eastern wall and 
consequently no radial pressure gradient. For 
flow through the wall and with Q=o we 
would immediately obtain the solution to (25) 
ve, = vez. Qualitatively we observe that for an 
interior source no isobars intersect the eastern 
wall but instead they form a region of high 
pressure extending to the western boundary 
with the source at the extreme eastern end. 
Having in mind this pattern the associated 
geostrophic recirculation and southward flow 
through the source region is obvious. 


7. Partial Barrier Experiments 


At the suggestion of Stommel, some experi- 
ments were performed to test the simple 
hypotheses proposed in (A) under somewhat 
more complicated conditions having geophys- 
ical interest. Figure 6a illustrates the anticipated 
circulation and Figure 6b the observed circu- 
lation in the 60 degree sector described in (A) 
but modified by a partial radial barrier 30 cm 
long extending inward from the midpoint of 
the rim. The barrier thus divided the sector into 
three portions referred to as the N, SW, and 
SE basins. Conditions of the experiment were 
Sg=100 cm/sec, h=8 cm, hy/h,=2, Q=1 
sec! and since no water was removed from 
the tank the sink is represented by the uniform 
rise of the free surface. 

The water budget of the predicted circulation 
is based upon the considerations developed in 
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Fig. 6. The circulation obtained with a source at the apex of the sector and a uni- 

form rise of the free surface in the presence of a partial radial barrier extending 

inward from the center of the rim. a) The theoretically predicted circulation; 

b) a sketch of the observed circulation; c) a schematic diagram for the definitions 
GING SIS, Spelled 


(A) as follows: we consider continuity in either 
the SW or the SE basin as indicated in Figure 
6c where T,, is the transport of the western 
boundary current into the basin, T; is the 
interior geostrophic transport out of the basin, 
and Ty is the upward flux due to the filling of 
the entire sector by the source. Continuity in 
2/42 
SW requires Ty = T1+ Ty. Ty= = a 
is the fraction of S, used to fill the SE or the 
SW basin. T; may be obtained from Ty= 
de 


=f fur dO dz where h is the basic height at 
9 o 


radius r and h’ is the departure from h due to 
the steady geostrophic current. From equation 


1 
(12) setting »=Q==~ =o we obtain u= — 


90 
Bo which permits evaluation of u. Using 
oh 2hr oh 65, en 
ee. and D Te obtain u = — SE 


and ignoring h’ compared to h we find Ti = 
Po 
al 


: S 
for Tr and Ty we obtain Ty= > (x +) 
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(: Hi): From! the above expressions 
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Therefore for I=ı we find Ty= So, and for 
r=70cm we obtain Ty= +.25 S,and Ty= +.75 
S, for either the SW or the SE basin. 

Of particular interest is the origin of the 
source water for the SE basin, for according to 
the simple theory 75 % of this water would be 
recirculated through the SW basin and 25 % 
would come zonally across from the western 
boundary current. This type of circulation was 
confirmed in the experiment by inserting 
crystals of dye in the SW basin to observe the 
flow around the end of the barrier as illustrated. 
Of further interest is the difference in history 
of the water in the western and eastern sides of 
the N basin; that on the west side having come 
directly from the western boundary current 
and 75 % of that on the eastern side having 
passed through both the SW and SE basins 
before flowing northward in the interior. 


8. The Limiting Role of Bottom Friction 


Figure 7a illustrates a partial barrier experi- 
ment similar to that of Figure 6 but with a 
localized sink near the rim as indicated and 
with the western and eastern walls coincident 
half way around the tank. The circulation 
performed in the expected fashion carrying the 
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Fig 7. A series of experiments illustrating the progressive domination of bottom friction 

as radial barriers are removed. a) A partial radial barrier experiment similar to that of 

figure 6; b) the same experiment as in a) but with the complete radial barrier removed; 

c) a spiral circulation pattern required to transport mass from Sp to S;; with no radial 
barriers. 


source water directly to the sink along the 
western boundaries. 

Removal of the radial wall (Figure 7b) left 
an annular ring with partial barriers, a geometry 
somewhat similar to the geography of the 
antarctic circumpolar seas. For this case the 
boundary current from the inner barrier moved 
westward and proceeded around the tank at 
the same time spreading in accordance with 
the theory presented in section IV. After 1 


circumference the dyed source water by-passed 
the outer barrier and the sink as indicated. 
Injection of dye beside the outer barrier clearly 
showed a boundary current flowing away from 
the sink. The interpretation is that the Ekman 
transport toward the rim was oversupplying 
the sink so that an inward flowing boundary 
current was required by continuity. 

An extreme example of the domination of 
the flow by bottom friction in the absence of 
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radial barriers is that illustrated in Figure 7c. Q 


The source and sink were at mid radii and 180 
degrees apart. The complex spiraling pattern is 
evidence of the importance of the bottom 
boundary layer in the absence of radial barriers. 


9. Applications to a Barotropic Ocean 
Circulation 


The example of recirculation about an interior 
source is closely related to the “wind-spun 
vortex” suggested by Munk (1950). In fact the 
theory of STOMMEL (1948) which includes a 
simple form of bottom friction and that of 
Munk (1950) which includes lateral friction 
at the western boundary together contain 
nearly all of the essential elements of the 
theory and experiments presented above. 

A striking feature of the recirculation about 
a source region is the extent of the large scale 
lateral mixing which may involve volumes of 
fluid several times the strength of the source 
itself. Also, the circulation of Figure 6 indicates 
significant mixing processes associated with the 
larger scale recirculations and the radial barriers. 
From this viewpoint it should be expected that 
the circulation of the deep ocean, which is now 
known to have some fairly strong currents 
(SwALLOW, personal communication), should 
be markedly affected by the large scale topo- 
graphy of the ocean floor. Accordingly, the 
distribution of properties may be significantly 
affected by localized mixing processes of large 
magnitude. 

The amplitude of the recirculated zonal cur- 
rents is of considerable interest especially in 
connection with the results recently published 
by ROBINSON and STOMMEL (1959). These au- 
thors described the localized amplification of 
the speed of barotropic currents of the oceans 
(Veronıs and STOMMEL, 1956) by relatively 
small scale irregularities in the ocean bottom. 
In a specific example, an amplification factor 
of 8 was suggested. Neither of the above named 
papers was concerned with nor implicitly in- 
cluded the recirculation phenomenon described 
here, and an additional amplification factor for 
barotropic currents should be considered. As in 
the model, the magnitude of the recirculation 
is primarily dependent upon the scale of the 
source-sink distribution as will be shown. 

The barotropic theory for an ocean of con- 
stant depth D is readily obtained by conversion 
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Fig 8. Conversion from polar to spherical geometry. 


from polar coordinates according to the geom- 
etry of the sphere as in Figure 8. Into equa- 
tions (19) and (20) we introduce the following 
relations for the sphere: 


r= R..cos 9, fre 
Sy 2 0 
rg DE 
We then obtain: 
sob-Qs(y)eRNa 
von East ann) (an 


and f 
T(+)= - fQs(y)dy+ 


0 


‘4 

+ R fd(Qs(y) tan y) (32) 
which may be evaluated for any given latitu- 
dinal distribution of the source. For a source of 
small latitudinal extent and symmetrical with 
a maximum at the mean position y, we may 
approximate equation (32) by letting tan p= 
tan Py to obtain: 


T(+)= = 294 R tan po Qs(y0) (33) 
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where Qs(y,) is the maximum value of Qs(y). 
As was done for the model we consider a nor- 


(v- Yo)” 
mal distribution Qs = | 20 and obtain 
Qs (yo) = as so that the recirculation is: 
7 
T(+)_. 1 Run 
Se a (270)? o (34) 


T+) 
3 
the source region when o<R. To take an 
extreme example we consider an idealized 
source due to precipitation or the wind stress 
convergence associated with a stationary hurri- 
cane having o=20 km at latitude p=25 de- 
grees. The second term of (32) is dominant and 


T(+) 
rm 


example, we consider a large-scale stationary 
system with o= 500 km at 9=60 degrees, and 
TCS) 

S 
formula (33). These figures suggest that the 
estimates of barotropic flow given by Veronis 
and Stommel (1956) and the partition of energy 
between barotropic and baroclinic currents 
(Veronis, 1956) should be reexamined in a 
framework which permits the recirculation 
phenomenon. 

It is of some interest to estimate the spread of 
zonal barotropic flow due to bottom friction 
for the scale of an actual ocean. From section 
IV the formula for the standard deviation of 
the spread is o,=(2Ay)}. For conversion to 
the spherical case we use 


depends inversely upon the spread of 


gives a recirculation = 9! As another 


we obtain =10 from the approximate 


0x dh Dcos y 
sing’ or R sin? 


, Q>Q2 sin y 


Oy 


in the Ekman equations, and »-> », (an eddy 
coefficient of viscosity). We then obtain 


R sin’ o ( “Ve " 
V2 D cos g \2 sin y 
For v,=1 cm?/sec, Q= .729 x 1074/sec-1, p=45 
degrees, D= 5 km, and R=6400 km, we find 
A= -45x10° cm, so that at a distance x= 


1000 km from a point source o,=42 km. We 
may take the total width of the current as 


A= 
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40,= 168 km. This width is latitude dependent, 
but it is reasonable to take 100 km as the mini- 
mum width of a barotropic zonal current except 
where local concentration due to bottom 
topography may be important. m, varies as v.t 
so that only order of magnitude variations of 
v, can be of significance here. 


10. Concluding Remarks 


From time to time the uses and applications 
of model experiments are seriously questioned, 
particularly with respect to inability toimprove 
our knowledge of real geophysical problems. 
The work described here contains a variety of 
experiments which suggest various possibilities 
and limitations for studies of the ocean circu- 
lation. 

In the first place, the experiments have con- 
firmed the qualitative predictions of the general 
ocean current theory suggested by STOMMEL 
(1957), and insofar as dimensional arguments 
can show that the same equations apply to the 
geophysical scale, the experiments represent 
confirmation of that theory. 

Secondly, laboratory studies using a real fluid 
(as opposed to an idealized frictionless fluid) 
are sometimes capable of resolving incomplete 
or inconsistent theoretical analyses such as those 
proposed by HouGx (1897) and by Gorps- 
BOROUGH (1933) each of whom considered 
source-sink circulations on a sphere. As has 
been pointed out by STOMMEL (1957) Hough 
had neither bottom friction nor meridional 
boundaries and could not obtain steady state 
solutions, whereas Goldsborough (having me- 
ridional boundaries) could obtain solutions only 
in certain restricted circumstances and did not 
conceive of boundary currents for the resolu- 
tion of the general problem. The experiments 
described here and in (A) clearly resolve the 
difficulties of the above authors by means of 
experiments which would have been ideally 
suited to examination of their theoretical ana- 
lyses. 

Finally, laboratory experiments may suggest 
circulations of possible oceanographic signifi- 
cance which have heretofore been overlooked. 
Such is the case for the recirculation phenom- 
enon and the spread of barotropic zonal flow 
each of which are contained within the elemen- 
tary theory. 

The eastern boundary current is an example 
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of a circulation which may be of some geophys- 
ical interest but which is not contained within 
the elementary theories and deserves special 
attention. It is tempting to suggest on the basis 
of our experimental results that inertial eastern 
boundary currents may occur where there is a 
sufficiently strong discharge of water from a 
river or other source at an eastern boundary of 
an ocean. A geographical situation which 
resembles the experiment of Figure 4 is the 
saline discharge of approximately 1,680,000 
m/sec (SVERDRUP, 1943) through the Straits of 
Gibraltar and into the Atlantic Ocean. It appears 
to be well established on the basis of observa- 
tions of salinity (NIELSEN, 1912) that at least a 
small portion of the Mediterranean discharge 
moves in a narrow current northward along 
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the coast of Spain, and a diagram of salinity 
would appear similar to the distribution of 
dyed source water in Figure 4. 

Complete velocity profiles are not available, 
but direct observations of currents using 
Swallow floats at the salinity maximum (1,100 
meters depth) off the south coast of Spain have 
shown on one occasion an average speed toward 
the West of 25 cm/sec at the closest observa- 
tional distance 7.5 km from the coast (F. C. 
FUGLISTER, private communication). A shear of 
25 cm/sec/3 km would equal the vertical com- 
ponent of the earth’s vorticity at that latitude, 
therefore it seems probable that some small 
fraction of the Mediterranean discharge may 
have sufficient shear to realize negative absolute 
vorticity. 
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The “Salt-Fountain” and Thermohaline Convection 


By MELVIN E. STERN, Woods Hole Oceanographic Institution 


(Manuscript received January 4, 1960) 


Abstract 


A ‘gravitationally stable” stratification of salinity and temperature, such as is observed 
in the oceans, is actually unstable due to the fact that the molecular diffusivity of heat is 
much greater than the diffusivity of salt. We discuss this stability characteristic and the 
form of the convective motion in the laminar regime. Future studies of this model rela- 
tive to the amplitude of the motion and the subsequent transition to turbulence should 
lead to the formulation of critical observational questions, which will determine whether 
the proposed mechanism is significant in the vertical mixing of the sea. 


STOMMEL, ARONS and BLANCHARD (1956) 
have described an ‘‘oceanographical curiosity”’ 
by noting that if a long vertical tube was 
lowered into the ocean, in such a manner that 
its bottom was exposed to cold fresh water and 
its top to warm saline water, a continuous 
motion could be maintained therein after 
priming the fountain. Their explanation is that 
the ascending (or descending) water in the 
tube would exchange heat but not salinity 
with the ambient ocean and would be acceler- 
ated due to its deficit in salt and density 
relative to fluid at the same level outside the 
tube. The purpose of this note, stemming 
from conversations with Henry Stommel, is 
to point out that in view of the great difference 
between the molecular diffusivity of salt 
(Ks = 1.3x 10-5 cm? sec-! for salinity of 
35 % at 20° C) and temperature RE Ts 
x 107% cm? sec~4) nature provides her own 
convective fountains. If a parcel of small 
radius is given a small vertical displacement it 
will lose its temperature excess much more 
rapidly than its salinity excess and the resulting 
buoyant force may be sufficient to maintain a 


free convection, despite the fact that the mean 
density field increases in the direction of 
gravity. The decrease of mean salinity in this 
direction provides the energy source for the 
convective motions, and if this salt gradient is 
maintained by climatological factors then it 
furnishes a mechanism, whose relative im- 
portance must be determined, for the vertical 
transport of salt and heat in the ocean. 

In order to forestall premature discussion of 
eddy excange coefficients let us first consider an 
isolated system consisting of a horizontal con- 
vection chamber whose top surface is maintain- 
ed at a higher temperature and salinity than the 
bottom surface. An equilibrium state consisting 
of linear variations of temperature and salinity 
may be produced, with the denser liquid below, 
and we shall inquire into the stability of this 
arrangement. This system is subjected to small 
perturbations and it will be shown that there 
is a direct, but superficial, analogy with the 
classical Rayleigh-Benard convection where 
the denser liquid is above. In either problem 
elimination of the horizontal velocity per- 
turbations and the pressure from the four 
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dynamical equations gives the following rela- 
tion between the vertical velocity (w) and the 


buoyancy force (ge/oy): 
(A = vv?) vw = — v2g0/00 (1) 


where » is the kinematic viscosity, v2 is the 
horizontal laplacian, and 1 the complex growth 
rate of the eigen-modes replaces the local time 
derivative. It will be assumed that all molecular 
parameters are constant, in particular the 
coefficients of temperature and salinity expan- 
sion, so that the equation of state is linear in the 
temperature (T) and salinity (S in gms/c.c.) 


perturbations.! Let EIER Br; AS Bs de- 


dz dz 
note the mean temperature and salinity gradi- 
ents in a direction (z) opposite to gravity. Local 
changes in temperature and salinity are brought 
about by advection of these mean fields and 
by diffusion, these relations being expressed 


by (2) and (3) below: 


(A- Krv?)T=Brw (2) 
(A-Ksv?)S= - Bsw (3) 
80/00 =g(aT - oS) (4) 


Equation (4) is the equation of state where 
o(&i cg.s.) is a function of the volumetric 
expansion of water due to a dissolved mass of 
salt. We first investigate the condition of margin- 
al stability and later consider internal waves. 
Thus, when À=0 eqs. (1)—(4) combine to give: 


vv8w=gv2 [SET ES lo (5) 


In the case Bs=0, Br >0 (5) applies to thermal 
convection. For the thermohaline convection 
problem we shall explicitly state that ofr is 
negative and greater in magnitude than of sso 
that the mean density field decreases with z. 


K 
However if oBs/lxBrl > a ~1/100 then the 
Mi 


bracketed term in (5) is again positive and when 
the bounding surfaces are perfect conductors 
of heat and salinity there is a direct corre- 


1 N.P. Fofnoff in an unpublished paper entitled ‘Energy 
Transformations and Stability in the Ocean” discusses the 
effect of the non-linearity in the equation of state on the 
mixing of water masses of different temperature and 
salinity. 
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spondence with Rayleigh’s problem. The analo- 


2 ; Rye ops aß TG 
gy is made by identifying the field] SP + | 
with a new thermal profile te and we there- 

ie 
fore conclude that the condition for marginal 
stability depends on the numerical value of the 
following thermohaline Rayleigh number: 


H4 
Rg epee - Se + | (6) 


where His the vertical height of the cell or the 
chamber. The minimum value of this quantity 
depends on the kinematical boundary con- 
ditions and in general is of order 108. If this 
number is exceeded then the diffusive equilibri- 
um state is unstable. It will be noted in the for- 
mal argument that we need not assume Ps 
and Br to be independent of z and therefore 
the marginal stability analogy ought to hold 
when these basic fields are non-linear, repre- 
senting quasi-steady diffusive states. An example 
of this occurs when warm saline water is placed 
above cold, fresh and dense water. Henry 
Stommel and Alan Faller have performed such 
an experiment and have kindly made available 
the notes of the results. 


Approximately 250 c.c. of warm water (38° C) 
were added to a graduate above 250 c.c. of water at 
(18°C), care being taken to avoid excessive mixing 
at the interface. The salt gradient was then produced 
by the addition of 2.5 c.c. of sea water to the upper 
layer, with a small quantity of menthylene blue dye 
solution added as a tracer. A fairly even distribution 
of the additive was produced by the mixing in the 
upper layer. It was immediately noticed that small 
elements of the dyed saline water began to extend 
downward into the denser region, taking the form of 
thin vertical and continuous filaments that eventually 
reached the bottom. There were four cells per linear 
inch while the speed of descent of the columns was 
about 1 X 107? cm sec~t. These regular structures were 
completely disrupted when a thermometer was inserted 
but re-established themselves within a few minutes. 
The experiment was repeated in a small aquarium 
(10” X 30” X 11” deep). Cold water (15° C) filled the 
bottom half and warm water (49°C) the top. A 
mixture of salt water and dye was introduced to give 
about the same density contrast as in the previous 
case. The qualitative results were the same as above. 


The following points may be noted from 
these observations. The stabilizing thermal con- 
trast was about twenty times greater than the 
destabilizing salt field, measured in terms of 
density alone. The mixed layer separating the 
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two fluids was of the order of three to six 
centimeters and therefore the thermohaline 
Rayleigh number was highly supercritical. The 
horizontal size and shape of the convective 
pattern are different from what one would 
expect by analogy with the thermal problem. 
A qualitative explanation of this latter point 
will be sought by comparing growth rates of 
disturbances of differing sizes at supercritical 
values of Rs+Rr. It will also appear that 
internal gravity waves are damped so that the 
criterion for marginal stability gives the mini- 
mum critical value for self-excited (infinitesi- 
mal) disturbances. In this calculation Bs and Pr 
are constant and “‘free-free’’ boundaries are 
assumed so that the vertical velocity, tempera- 
ture and salinity perturbations are of the form 
sin (maz) exp mi (kx + Ly), where mis an integral 
multiple of the reciprocal height of the chamber. 
Substituting these in (1)—(4) gives the following 
cubic equation for A: 


a +1 2 +r Maus 
\ Ks a Ksa Rearing ye 


À 
=s[¢-D)-zae-n] © 
where 
fé go Bs a _ _ «Br re 
vKs ata? oßs Ze Kr 
FT az = 0? (k? +) = a? - n?m? 
s 


To determine whether neutral gravity waves 
can exist substitute A=i;K,a? into (7). Equat- 
ing real and imaginary parts to zero gives 


(+) er f(r D) (8) 
= Np[f(D-1)+1+r]+r (9) 


By eliminating the frequency (A;) it may readily 
be shown that these equations are incompatible 
if D > 1 and f > 0, from which it is concluded 
that internal waves are damped? and there can 


! Compare the interesting case of a ‘solute’? which is 
distributed so as to stabilize a super-adiabatic lapse rate 
of temperature. For exemple, if f<o, r>D>ı 
(i.e. By <0, By > 0) the salt can stabilize the stationary 
mode and destabilize the oscillatory mode, thereby releas- 
ing the potential energy in the thermal stratification. 
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only be one root whose real part is positive and 
this root must be a real number. 

The approximate value of this root is now 
computed under the conditions: r> D> 1, 
Np>1. In the experiment mentioned above 
Np=7, r=120, D~20. First note that the left 
hand side of (7) is positive and therefore the 


right side implies the upper bound ok 


< aan ~ = Let O (1/D) denote a term which 
is small compared to unity by order 1/D. Then 
we approximate the cubical expression as 


follows: 
(a) (ea) er)” 
)]H+o&D) 


=| r+ 2 (1 +r) REST 
Ww Ks a? Itr 


= [r+ ante + 0(2/D)} | + O(1/D)] 


~r+ nr +r) 
Ks a 


Equating this to the right side of (7) gives 


À. EDIT DM 
Ksa® r+ı+f(D-ı) 1+fD/r 


and when the definitions of D, f, r are substi- 
tuted above one obtains 


À (a3 + mn?) (goBs/vKs) (1 - D/r) 
Ks az (a a3 mr) 35 |gxBr/vKrl| . 
n az (a? + m? 72) feel 
az "(a+ mr) + |gxBr/vKr| 


We already know that at marginal stability the 


size of the cell is me a =O (m). At super- 


critical conditions it follows from (10) that 
disturbances of this size grow at the rate 
AH?/Ks~ Rs|Rr|-!=rD-!. Compare this to 
disturbances which have the much larger hori- 
zontal wave number: a=c|gxß ma where 
c is an arbitrary number of order unity. Since 
ner and Rr>10° eq. (10) simplifies to 
a Durs Rs|Rr|”*,, which has a maxi- 
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mum at c=1. Amongst the class of disturbances 
whose horizontal plan-form consists of close- 

acked squares, the one which grows the fastest 
Fr a side of length 


(11) 


In the first experiment which was cited the 
temperature difference between top and bottom 
was 20° C and the transition region between the 
two fluids was 3—6 cm. Substituting the mean 
gradient into (11) gives a value of L between 
0.2 and 0.3 cm, as compared with an observed 
value of about 0.3 cm. 

It appears then that although cells having a 


horizontal dimension comparable with H are 
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preferred just above the critical point, their 
large size requires that the vertical velocity be 
kept small enough so that the temperature 
excess may be diffused. As the Rayleigh numbers 
are increased, with the salt-temperature ratio 
held constant, the thin columns can reach 
greater velocities and more effectively release 
the potential energy in the salt stratification. 
This explanation was suggested by Willem 
Malkus and I should like to acknowledge the 
fruitful discussions with him and with Bert 
Bolin. 
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Remarks on Correlation Methods in Geophysics 
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(Manuscript received July 29, 1959) 


Abstract 


From consideration of several studies of the sampling properties of the product-moment 
correlation coefficient r it is concluded first that undue concern has been expressed for the 
problem of non-normality in correlation studies in geophysics, and second, that use of geo- 
physically adequate significance tests and confidence limits for r can almost always be achieved 
through use of the simple standard error of r. A growing tendency in climatology and hydrology 
to employ unnecessarily elaborate methods appears to stem from unrealistic emphasis which 
mathematical statisticians frequently place upon theoretical refinements, emphasis that loses 
sight of the limits of accuracy inherent in the very type of data usually subjected to statistical 


analysis. 


I. Introduction 


This paper has two principal objectives: 
First, to question a rather common tendency to 
employ hollow refinements in correlation 
analyses in geophysics; and second, to display 
some empirical results in correlation analysis 
that should be of interest to the geophysicist 
seeking examples that may guide him through 
what is to many a rather bewildering array of 
warnings in the statistical literature on corre- 
lation methods. 

I shall try to show that in most of the correla- 
tion analysis that investigators carry out in geo- 
physical studies especially in the fields of 
meteorology, climatology and hydrology, sta- 
tistical elaborations beyond the simplest and 
most familiar techniques and tests cannot ordi- 
narily be expected to yield refinements of real 
geophysical significance. In addition, I shall 
present arguments for concluding that the pit- 
falls said to be inherent in correlating non- 
normally distributed variates are not likely to 
warrant serious concern on the part of the 
geophysicist. 


Three references concerned with the applica- 
tions of statistical methods in geophysical 
problems (chiefly climatological problems) 
have recently appeared, that of Conrap and 
Potrak (1950), that of Brooks and Car- 
RUTHERS (1953) and that of PANorsky and 
Brier (1958). The first offers only the barest 
details of the sampling theory of the product- 
moment coefficient, yet it is the writer’s view 
that, in most if not all climatological and 
hydrologic correlation analyses, the investa- 
gator is probably well-advised to remain little 
more sophisticated in his statistical treatment 
of his data than are CONRAD and POLLAK in 
their treatment of correlation. Specific reasons 
for this perhaps unorthodox view will next 
be summarized. briefly. 


2. The normality problem 


First, it must be noted that in none of the 
three recently published references cited above 
is it made quite clear that all of the product- 
moment correlation theory ee applies 
rigorously only to normal bivariate popula- 
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tions. On the other hand, the non-statistician 
who consults the literature of mathematical 
statistics is almost invariably disturbed at every 
turn by warnings against applications of nor- 
mal-theory statistics to non-normal popula- 
tions. However, the reader of these latter 
references is seldom given the slightest inkling 
of how large an effect non-normality in his 
data might have on the final conclusions he 
may draw, so it is regrettable that this practical 
question is not confronted in any of the three 
cited references dealing with applications to 
meteorology and climatology. 

“Effect of non-normality” here takes on two 
slightly different meanings associated with two 
statistical questions: First, the investigator 
wishes to know whether a computed coefficient 
constitutes an unbiased estimate of the pop- 
ulation value of the correlation coefficient if 
derived from non-normal data, and second he 
wishes to know whether non-normality can 
seriously interfere with significance tests if 
based on normal-distribution theory. STIDD 
(1953) has suggested, for example, that most 
precipitation data must be regarded as suspect, 
on grounds of non-normality, with respect to 
exactly such questions as these two. 

Concerning the first of these two questions, 
that of bias in r, it is often suggested, though 
typically only in vague terms, that possible 
bias due to non-normality can be avoided by 
use of one or another “normalizing trans- 
formation”. Strpp (1953) urges that in correlat- 
ing precipitation data one use a cube-root 
transformation to suppress what he reports as 
a negative bias in r due to the well-known non- 
normality of most rainfall distributions. 
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The literature contains nothing approaching 
a definitive analysis of the advantages to be 
gained in correlating data normalized by such 
transformation, so it seemed a useful exercise 
to carry out at least a modest number of prod- 
uct-moment correlations on some actual geo- 
physical data of markedly non-normal charac- 
teristics, first using raw data and then using 
each of several transformations thereof. Since 
the data were precipitation data and since the 
investigator may find in the literature, in turn, 
suggestions that such data be transformed by 
square root, or cube root, or logarithmic trans- 
formations (among others), I have tested each 
of those three. The test data were all taken 
from long-record U.S. Weather Bureau sta- 
tions in Arizona. It is well known that pre- 
cipitation frequency-distributions from arid- 
region stations such as those in Arizona are 
notoriously non-normal and usually highly 
skewed. Hence the present analysis should con- 
stitute a fairly extreme test of the true value 
of various normalization transformations. 

The results are displayed in Table 1. Seven 
station-pairs, with record-lengths varying from 
so to 83 years are included in this sample. 
Their records were broken down into “winter” 
and “summer” half-years, “winter” being the 
period from November through April, and 
“summer” the balance of the year. The prod- 
uct-moment correlation coefficients for the 
untransformed (raw) data are indicated by 
“linear”, and the figures in parentheses fol- 
lowing these values are the associated standard 
errors. Table ı represents correlations com- 
puted for fourteen independent bivariate dis- 
tributions, and of the forty-two coefficients 


Table 1. Effect of three normalizing transformations 


nern 
No. of years 
recorder et: 64 76 
RE on cn ao oc 54 (.09) | .38 (.10) 
Square root..." .55 .41 
Summer | Cube root....... 57 .42 
Bose 66 do vue 57 42 
IDM ag con oo bot 78 (.05) | .82 (.04) 
Square root. =... 80 85 
Winter | Cube root....... 72 .78 
ILO 20.00 74 .76 


1 Key to station pairs: A— Tucson, B— Natural Bridge, C — Phoenix, 
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AD A SEC DE | nee 
50 83 64 50 63 
AA 2) TO (GEL) OTELES) FEZNLE8)" | 22.12) 
47 .24 .65 .62 .36 
.46 .20 .65 .63 .42 
.38 .28 .69 .57 .48 
58 (.09) | .65 (.06) | .85 (0.3) | .57 (.10) | .68 (.07) 
.60 .62 85 .57 .65 
52 .63 .83 .60 .66 
.59 OT .81 .62 .64 


D — Flagstaff, E — Yuma. 
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computed from the three types of transformed 
data it will be seen that only six values are even 
as much as one standard error away from the 
linear r. Only one is more than two standard 
errors away from the linear value (Natural 
Bridge-Yuma log data). The results of Table 1 
argue strongly that only negligible gain is to 
be achieved in subjecting even highly non- 
normal data to normalization transformations 
as far as any effects on r are concerned. If such 
minor changes in r-values are produced by 
applying these several normalizing transforma- 
tions to arid-zone precipitation data, it seems 
clear that in analysis involving less skewed 
distributions from more humid regions, no 
significant changes in the value of r could result 
from use of these transformations. STIDD (1953) 
took note only of the fact (confirmed in the 
present studies) that the cube-root transforma- 
tion does symmetrize the frequency-distribu- 
tion noticeably, but he seems only to have 
inferred that this essentially graphical change 
would be accompanied by significant changes 
in values of r computed from transformed data. 
The present results indicate that the latter 
changes are of no consequence.! 

With respect to the second question, cited 
above, that of making significance tests for 
correlation coefficients derived from non- 
normal bivariate populations, the geophysicist 
should realize that theory is at present incapable 
of yielding decisive answers. Statisticians have 
consequently resorted to experimental sampling 
to secure estimates of this effect of non- 
normality; but this work was done so many 
years ago that it seems to have been lost to the 
present attention of most geophysicists. Hence 
I feel that a few examples of such studies 
warrant brief recapitulation. A substantial 
amount of such experimental sampling was, 
for example, done by E. S. Pearson and his 
collaborators early in the 1930's, with the 
objective of testing non-normality effects on a 
number of statistics, including the product- 
moment correlation coefficient. 

In one interesting study PEARSON (1931) set 
up four different parent bivariate populations 
of pure numbers distributed with varying 
degrees of non-normality. Although all of his 


1 For additional data supporting this conclusion, see 
the recently published analysis by Huff and Neil, J. Geo- 
phys. Res. 64, pp. 541—547 (1959). 
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distributions were much smoother, than, say, 
typical precipitation or streamflow frequency- 
distributions, they exhibited skewness and 
kurtosis of a degree that the climatologist 
would regard as rather severe if found, for 
example, in precipitation data. From these 
synthetic populations, Pearson experimentally 
drew 250 successive samples of 10 pairs each 
and then 250 more of 30 pairs each, except for 
one case where he used 395 samples of each of 
these two sizes. In each instance the sampling 
scheme corresponded to sampling from a parent 
bivariate non-normal population of zero 
correlation. From his eight series of samples he 
formed frequency-distributions of r and from 
these he determined the experimental “standard 
error” for comparison with the corresponding 
value computed on the assumption that normal 
theory could be used. For samples of 10 items 
and a “normal theory” standard error of 
0.3333 (the choice of using four digits is 
Pearson’s, not the writer’s), Pearson found 
experimental values of 0.3234, 0.3487, 0.3238, 
and 0.3181. Similarly, for the samples of 20 
items each and a “normal theory” standard 
error of 0.2294, he found experimental values 
of 0.2413, 0.2405, 0.2289, and 0.2308. His 
conclusion should be familiar to geophysicists 
who make use of correlation techniques: 
“Taken as a whole, the extent of agreement 
has proved considerably greater than was 
anticipated before the experiment was com- 
menced, and it seems to indicate that the 
coefficient of correlation is another of the 
criteria based on ratios, whose distribution 
even in very small samples is remarkably in- 
sensitive to changes in the form of the popula- 
tion.” Indeed, one must ask, what geophysicist 
would not view the above agreement between 
Pearson’s normal-theory and experimental 
standard errors as near-perfect? 

In a somewhat similar experimental sampling 
study, BAKER (1930) drew 50 samples of 40 
items (pairs) each taken from a highly skewed 
population. This case, like the one just cited, 
constitutes a test of the case of zero correlation 
in the parent population. Baker drew the 
conclusion that his resultant distribution of 
r-values was “very skew’’, but the geophysicist 
must not be too quickly misled by such 
mathematical-statistical pessimism, since in- 
spection of his results reveals that 29 of the so 
lay within one theoretical standard error of 
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zero, and an additional 15 cases lay within two 
standard errors of zero. This constitutes dis- 
cernible discrepancy but hardly geophysically 
alarming asymmetry. The remaining six of the 
fifty samples fell above (more positive than) 
two standard errors in excess of zero, pre- 
sumably a consequence of the high probability 
of independently drawing two items such that 
both lie below the mean of a positively skewed 
parent population. The conclusion of relevance 
here is that Baker’s study tends clearly to sup- 
port the view that non-normality of even a 
marked degree is not so large that it would 
seriously mislead the geophysicist in most 
problems wherein correlations are used, even 
though Baker, as a statistician, took his results 
as empirical evidence of the seriousness of the 
non-normality problem. 

HALDANE (1949), in a note on the problem of 
non-normality in correlation, summarizes his 
and others’ sampling studies as having led to 
the conclusion that “the normal bivariate sur- 
face may be distorted and mutilated to a 
remarkable degree without seriously influenc- 
ing the frequency-distribution of “r” in the 
case of small parent-population correlation 
(i.e., in the case of importance in making 
significance tests for mere existence of correla- 
tion). 

To see the kind of results that arise in ex- 
perimental sampling from parent bivariate 
non-normal populations with non-zero correla- 
tion, geophysicists should refer, for example, 
to the work of CHESIRE, OLDIs, and PEARSON 
(1932). For the explanation of their sampling 
model, the original paper must be consulted; 
here it is only necessary to note that these 
workers set up two populations with known 
correlation coefficients equal to 0.5000 and 
0.5462 (number of digits chosen by those 
authors) and drew 1,000 samples of 5 pairs 
each, soo samples of 10 pairs each, and 250 
samples of 20 pairs each. One population had 
symmetric triangular frequency-distributions 
for both variates, the other population had one 
triangular symmetric and one asymmetric 
distribution. “Normal theory” values of r and 
of its standard error were computed and com- 
pared with experimental values. Because the im- 
plications of these results are of practical con- 
cern to geophysicists and to the thesis of the 
present paper but are to be found only in a 
paper that is probably not generally familiar, 
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comparative values for just the cases of 10 and 
20 pairs per sample are displayed in Table 2. 


Table 2. Observed and theoretical values of cor- 
relation coefficient r, and standard error, o; for 
500 samples of Io pairs and 250 samples of 20 
pairs drawn by CHESIRE et al. (1932). 


Samples of Samples of 
10 pairs 20 pairs 
14 II | af II 
Obs. .4871| .5289| .4896| .5460 
Mean 7 
ala, .4787| .5242| .4900| .5359 
Obs. 2371-2303 1701161503 
o 
d Th. .2671| .2534| .1780| .1677 


1 I and II designate the two populations samples, 
II being composed of a triangular asymmetric and an 
asymmetric subpopulation (see text). I and II had parent- 
population correlation coefficients of 0.5000 and 0.5462, 
respectively. For deduction of “theoretical” sampling 
values see original paper. 


The authors conclude that there are “real 
differences between the distribution of r in 
samples from these two populations and those 
appropriate for the normal case”, but admit 
that one’s view as to whether these differences 
would be considered serious in practice depends 
upon the degree of accuracy required in partic- 
ular investigations. Putting the author’s apo- 
logies aside as sheer fussiness, I would be forced 
to ask: What geophysicist would see any 
significance in the differences displayed by the 
pairs of observed and theoretical correlation 
coefficients in Table 2 if he encountered these 
in his own research? 


3. Tests for existence of correlation 


If the question of normality has been decided 
in favor of believing that it may be ignored 
there next arises the question of whether a com- 
puted product-moment coefficient may safely 
be taken to indicate that the parent De 
is characterized by non-zero correlation, i.e., 
one next wishes to make a probability state- 
ment relative to the existence of correlation 
(without regard to “true” magnitude). CONRAD 
and PoLLAK (p. 245) quote a rule-of-thumb 
criterion which used to appear in the literature 
of meteorology and climatology fairly fre- 
quently, namely the rule that a coefficient must 
equal or exceed twice its standard error if it 
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is to be taken as evidence of the existence of 
correlation in the parent population ?? We now 
regard it more fashionable to invoke arguments 
based on the t-test. To be sure, it takes only a 
little more time to use a f-test than to argue 
on the basis of a standard error, yet with 
reference to the main thesis of the present 
paper it is relevant to observe the rather un- 
important degree of refinement in all but 
rather small samples (in less than about 20) 
which this introduces into one’s statistical in- 
ferences. In Table 3 are shown, for several 
values of sample size N, values of the product- 
moment coefficient r that will be just signif- 
icant at the 95 per cent confidence level based 
on the t-test, as described, for example, by 
Brooks and CARRUTHERS (1953, p. 220). Also 
shown, for comparison, are corresponding 
values of r, for the same sample sizes, that would 
be regarded as significant at about the same 
level on the basis that r must then equal or 
exceed twice the standard error as given by 


o,=1/(N -2)': 


Table 3. Values of r significant at the 5 per cent 
level as derived for two criteria. 
| 


Sample Size N (x)? (2)? 


Io .63 al 
20 -44 47 
40 31 32 
60 .2 26 
80 22 2 
I0o .20 21 


1 Values of rin Column (rt) are based on the t-test. 

2 Values of rin Column (2) equal twice their respective 
standard errors, i.e., are solution of equation (r) for 
r=20, 


It can be seen from Table 3 that if one used 
only the old rule-of-thumb criterion, he would, 
for allbut very small sample-sizes be employing 
almost exactly the same threshold values that 
he would derive from the t-test for r applied 
at the 95 per cent significance level. For samples 
much smaller than 20 variate-pairs the differ- 


? The philosophical meaning of the phrase “parent 
population” is, of course, obscure when sampling time 
series, but is to be understood here in common-sense 
terms. There is evident need for closer grappling with 
this difficulty that underlies so much of climatology’s 
and hydrology’s use of statistics. 
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ence between the two criteria becomes 
numerically appreciable, but then it seems 
relevant to note that the basically arbitrary step 
of selecting 95 per cent rather than 90 or 99 
per cent as the working significance-level in- 
troduces a difference about as large as the ob- 
served difference in threshold r-values of the 
two columns of Table 3. So light a treatment 
of the difference between, say, a 95 and a 99 
per cent significance-level will, presumably, 
shock the statistician, but except in contexts 
where very important decisions hinge on inter- 
pretations of stastistical results and where the 
basic data are of high accuracy the statistician 
cannot really defend, I believe, his typical point 
of view. I would contend that only in such 
atypical instances as weather modification 
experiments does one really have any reason to 
push statistical sophistication to the limit. And 
even in that example, how many times the in- 
investigator finds himself embarrassed by the 
usually woeful inadequacy of observational 
data upon which to work with the more 
refined statistical tools. 


4. Determination of confidence limits 


The third basis for suggesting that many of 
the niceties of sampling theory of product- 
moment correlation coefficients need not 
greatly trouble the geophysicist who has not 
been able to become thoroughly familiar with 
them concerns the final problem of placing 
confidence limits on a given computed correla- 
tion coefficient. It is theoretically important to 
distinguish such a process from that just dis- 
cussed, for the sampling distribution of cor- 
relation coefficients drawn from a parent 
bivariate population of zero correlation is 
itself normal, but the corresponding distribu- 
tion for a parent population of non-zero cor- 
relation is skewed and (still speaking theoreti- 
cally) this affects significance-testing procedures. 
In almost any reference which the geophysicist 
may consult for aid in interpreting his correla- 
tion analyses, he will be confronted with dis- 
quieting admonishments against using a simple 
standard-error argument in setting his con- 
fidence limits on r. If, being sufficiently dis- 
turbed by these admonishments, he reads 
further, he is led to use Fisher’s hyperbolic 
arctangent transformation, usually called simply 
the 2’-transformation. 
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The technique whereby one transforms r 
into 2’, places selected confidence limits above 
and below this value of 2’ by what amounts 
to a f-test, and then makes the inverse trans- 
formation is well explained in many statistics 
texts, so it requires no comment here. The in- 
terested reader will find particularly clear 
discussions given by OsTLE (1954, pp. 181—185) 
and Mints (1955, pp. 297—309), and a handy 
graphic aid for effecting the r—z’ transforma- 
tion and its inverse is presented by SNEDECOR 
(1946, pp. 1952—53). However, before acquir- 
ing familiarity with the z’-transformation, the 
geophysicist should ask whether use of this 
refinement of sampling theory is justified by 
the usual levels of precision encountered in his 
statistical studies. Table 4 provides information 
that I believe calls for a negative answer. It is 
evident from Table 4 that using plus and minus 
twice the standard error of r (strictly speaking 
one should use 1.96 times the standard error, 
to secure a 95 per cent half width, but such a 
2 per cent refinement is no more materially 
important than the other refinements here 
under discussion) yields estimates of the 95 


Table 4. Comparison of 95 pct confidence half- 
widths of r for two sample-sizes N 


N = 80 
r 


O.10| .46 .42 .46 Ze 2% 22 


0.20] .44 .39 -47 2 .20 22 
0.30 .42 .36 .46 .20 .19 2 
0.40| .39 .31 .45 .19 RT, .20 
0.50| .35 27 .43 "7 
0.60 .30 22 .39 .I4 2 .I6 
070 RE BUG 33 Mes .10 19 
O:60| T7 NEZ .25 .08 .07 .10 
0.90 .09 .06 .I4 .04 .04 .05 


14 = 207. The quantity d’ is the upper half-width of 
the 95 per cent confidence interval as derived from the 
z’-transformation, and d” is the lower half-width of the 
95 per cent confidence interval as derived from the 
z’-transformation. 


per cent confidence half-widths that are, even 
in the cases of poorest agreement, close enough 
to those obtained by using Fisher’s z’ that 
geophysical inferences will seldom be seriously 
distorted by the discrepancy. Since use of the 
95 per cent rather than, say the 90 or the 99 
per cent confidence level is itself basically an 
arbitrary choice made by the investigator, and 
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since substituting either one of these levels for 
the 95 per cent level will alter the half-widths 
by as much as or more than the replacement of 
the z’ method by the two-standard-errors rule, 
little practical justification can be given for use 
of this circuitous transformation in most geo- 
physical studies. The statistician who takes 
issue with so easy an acceptance of this level of 
approximation will have to show that there 
are geophysicists who really do make signif- 
icantly different decisions depending on 
whether their significance-levels come out at 
one per cent or five per cent. I believe they are 
very rare indeed. 

Note that, in general, double the standard 
error gives a half-width intermediate in value 
between the upper and the lower half-width 
so that, in a sense, 20, constitutes a single figure 
rather nicely summarizing the two asym- 
metric limits of refined sampling teory. The 
statistician will, of course, object that it is 
exactly this asymmetry that makes it necessary 
to use the z’-transformation to correct for 
skewness in the r-distribution, and his position 
is qualitatively invulnerable. 

The investigator who has not yet had the 
opportunity to familiarize himself with Fisher’s 
2 -transformation should not hesitate to con- 
tinue employing the simple standard error 
argument unless he feels that he is handling 
geophysical data of such rare precision and 
freedom from observational-sampling varia- 
bility that errors of estimate of a few hun- 
dredths in his confidence limits are going to 
lead him to erroneous inferences. But any such 
unusual data are likely to be of such functional 
simplicity as to obviate resort to statistics in 


the first place. 


5. Historical interpretation and summary 


All that has been said above can be con- 
cisely summarized by saying that investigators 
dealing with climatological, hydrologic or 
other complex geophysical data should take a 
commonsense view of ordinary correlation 
analysis as seen in the context of data-precision 
and should, except in unusual cases, feel quite 
justified in using only the simplest of classical 
sampling theory for the Pearson product- 
moment correlation coefficient. 

Since this thesis may seem retrogressive to 
some, it may be relevant to attempt an explana- 
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tion here of why so many non-statisticians have 
acquired the tendency towards overelaboration 
that so often appears in the statistical discus- 
sions in geophysical papers. It seems to the 
writer, as a result of having recently had the 
curiosity to scrutinize about two dozen statis- 
tics texts and references in relation to this 
point, that a difficulty arises from the fact that 
statistics texts reflect too heavily the viewpoint 
of the theoreticians, and I believe it is fair to 
say that the latter have in recent decades (es- 
pecially since the early work of R. A. Fisher) 
simply pushed their sampling theories in direc- 
tions and to degrees of refinement that are often 
entirely out of proportion to levels of data- 
precision typical of those very areas of research 
which, for reasons of inherent complexity and 
inherent variability of subject phenomena, 
must involve statistical methods. (This his- 
torical speculation must not be taken too 
broadly, for there has certainly been much 
recent theoretical work designed to meet more 
practical investigative needs. Since these more 
recent advances are not well covered in the 
references typically contained in geophysical 
bibliographies, one sees a need for critiques 
on the application of these newer techniques in 
the geophysical sciences.) 

By way of illustration of my above assertion 
that many statisticians have tended to work 
and write in an atmosphere of false precision, 
I would cite the following: In the course of 
preparation of this paper, as I examined in- 
creasing numbers of statistics texts and ref- 
erences, I became aware of a surprisingly 
general tendency for writers of statistics texts to 
compute and cite correlation coefficients to 
more than two significant digits (surely the 
limit that a sense of numerical proportion 
would seem to indicate). Systematic inspection 
of fifteen widely used and quoted statistics 
texts revealed, in fact, that eleven out of 
fifteen writers used three or more digits for r, 
and fully a third gave coefficients to four or 
more digits! Such practices may offer indirect 
explanation of an almost incredible but true 
instance wherein a non-statistician associated 
with a hydrologic project came to the writer 
with a fifty-year basinwide rainfall-runoff 
correlation equal to 0.946 when all years were 
included, and 0.957 when a single very wet 
year was excluded, and asked if this might be 
taken to indicate that a proposed major hydro- 


JAMES E. McDONALD 


~ 


logic experiment on that basin could be tested 
statistically by comparing the correlation 
coefficient for all years up to but not including 
the year of treatment with the coefficient 
found for that data plus the data for the year 
immediately following treatment! That statis- 
ticians ought not too quickly smile at such an 
interpretation of the significance of the digits 
in r is shown not only by the curious results of 
the previously-mentioned tally of fifteen tex- 
tual practices on significant digits but also by 
many statistically recommended but almost 
meaningless practices such as those of the ap- 
plication of “‘attenuation corrections” to r, or 
adjusting r via Sheppard’s correction, or a 
particularly astonishing case which is to be 
found in a text by SMITH and Duncan (1945, 
p. 303) wherein the reader is shown how an 
original correlation coefficient of 0.89576 can, 
by suitable algebraic manipulation, be con- 
verted into the Fisherian ‘‘maximum-likeli- 
hood estimate” of 0.89465 ! Clearly, statisticians 
who would write out r to five digits and then 
theorize on a means of improving it by one 
part in nine hundred are not addressing them- 
selves to the problem of statistical manipula- 
tion of raingauge readings taken in 1908 under 
a rapidly growing shade tree or of streamflow 
data for a month containing a record flood 
which the observer was unable to gauge at 
crest because the access roads were submerged. 

Within the literature of geophysics itself, 
there is no difficulty in finding writers who 
quote r to four places, but as a recent glaring 
instance of unreasonable statistical purism 
within our own field, I would cite viewpoints 
taken by Levert (1958), who criticizes another 
investigator’s failure to test for normality and 
failure to “take care not to use the common 
Fisher’s significance table” in a hydrometeoro- 
logical study that yielded r = 0.99 for a sample 
of 140 pairs! It is very easy, and I speak with 
some experience here, to fall into the error of 
pursuing statistical refinement far beyond one’s 
original goal. But how often does the investi- 
gator or his informed reader have any basis for 
looking beyond the first decimal place in r? 

The geophysicist who has been consciously 
or unconsciously influenced by schools of 
statistical thought so divorced from awareness 
of basic physical and observational uncertain- 
ties may be well-advised to review and simplify 
all of his statistical practices along lines similar 
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to those here followed for the product-moment 
correlation coefficient and its sampling prop- 
erties. If for no other reason than to reduce 
the burden of methodology to be assimilated 
by younger workers in the field of geophysics 
sciences, all superficial complexity of statistical 
methods in geophysics is to be discouraged. 
Where refinements are truly in order, as for 
example, in the evaluation of cloud-modifica- 
tion experiments, they deserve careful atten- 
tion. But only confusion arises from over- 
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sophistication in statistical manipulation of 
data in most geophysical studies. 
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Abstract 


Data on trace substances in air, obtained in mid-Pacific during two summers, are summarized. 
Particulate matter was analyzed for chloride, sulfate and nitrate; tests were performed for ozone, 
nitrogen dioxide, sulfur dioxide and carbon monoxide. The results give information concerning 
the composition of marine air which can assist in evaluating possible natural contributions to 
urban pollution on the West Coast of the United States during the summer months, in addition 


to the geochemical value of the work. 


I. Introduction 


In spite of a substantial effort presently under 
way in determining the composition of polluted 
atmospheres in cities, little effort has been expend- 
ed on the study of natural levels of pollutants 
in influent air, or on the contribution of nature 
to the overall pollution complex of cities. The 
state of our knowledge with regard to gaseous 
pollutants of natural origin has recently been 
reviewed by ALTSHULLER (1958). Studies of the 
composition of marine aerosol have been made, 
among others, by Woopcock (1949, 1950, 
1952, 1953), JUNGE (1953, 1954, 1956, 1957), 
Twomey (1954), LODGE (1955), and by Bygrs, 
SIEVERS, and TurTs (1957). These, however, we- 
re made for special purposes, and do not contain 
the information necessary, in most cases, to 
evaluate the specific contribution of marine 
acrosols to urban pollution. Probably the most 
useful studies in this connection have been made 
1 Assistant Professor of Engineering, San Jose State 
College, San Jose, California. 
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by the group in Sweden (EGNÉR AND ERIKSSON, 
1955, and quarterly reports by a variety of 
authors in Tellus), who have analyzed both 
precipitation and air for a number of substances 
using samples from a dense network of stations 
covering the Scandanavian peninsula and ex- 
tending into contiguous portions of Europe and 
Britain. It may be questioned, however, to what 
extent these stations ever sample unpolluted air. 
The National Air Sampling Network (Public 
Health Service, 1958) now provides particulate 
analyses from a number of nonurban sites, but 
these are subject to the same doubt. 

In 1956, the U. S. Public Health Service and 
the Department of Public Health, State of Cali- 
fornia, undertook a series of joint studies of the 
composition of air at coastal sites in the State 
of California (HOLZWORTH, 1959; HARRISON 
and Lopez, 1959). Even at these sites, consider- 
able evidence was obtained of man-made con- 
tamination. Accordingly, plans were made to 
take measurements from Weather Station 
November, the weather ship maintained by the 
U.S. Weather Bureau and the Coast Guard at 
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Table 1. Techniques used in air quality measurements in mid-ocean. 
T ; 
Substance Tau Reference 
1957 1958 
Oxidant (ozone) 200% Ka) Hagen-Smit, private 
communication 
Phenolphthalin McCabe, 1953 


Nitrogen dioxide 
Sulfur dioxide 


Modified Griess reagent 
Tetrachloromercurate 
method 

Carbon monoxide 


Particulate matter ”High-Volume’ sampler 


wiht glass filter 


(a) Data invalidated by photolysis of reagent. 


a position approximately midway between San 
Francisco and Honolulu (latitude 30°N, longi- 
tude 140°W). 

Two separate series of measurements were 
made. The first set spanned the period from 27 
July to 13 August 1957, during which time 
the station was manned by the USCGC Gres- 
ham.The second set of samples was taken during 
the duty period of the USCGC Taney, from 
31 July to 28 August 1958. Measurements 
were made of oxidant, nitrogen dioxide, sulfur 
dioxide, carbon monoxide and of particulate 
chloride, sulfate, nitrate and organic material. 


2. Experimental 


Some different techniques were used during 
the two seasons of operation. The sampling 
sites were in general on boat deck level; they 
were moved as necessary to avoid contamina- 
tion from stacks, ventilators, etc. Table I out- 
lines the basic techniques used in the two years. 
The analytical methods for gases were checked, 
and the reagents were prepared and calibrated 
by the Air and Industrial Hygiene Laboratory, 
Department of Public Health, State of Cali- 
fornia. The particulate analyses were performed 
in part by that laboratory and in part by the 
laboratories of the National Air Sampling Net- 
work, Air Pollution Engineering Research, U.S 
Public Health Service. The analytical methods 
used at both laboratories were in the main those 
of the National Air Sampling Network (Public 
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Rubber cracking Bradley and Haagen- 
Smit, 1951 

A. W. Brewer, un- 
published work 
Saltzman, 1954 


West and Gaeke, 1956 


Brewer-Mast ozone 
recorder 

Modified Griess reagent 
Tetrachloromercurate 
method 

Indicator tubes (Mine 
Saftey Appliances) 

” High-Volume” sampler 
with both cellulose and 
glass filters 


Shepherd, 1947 


Public Health Service, 
1958 


Health Service, 1958). The Air and Industrial 
Hygiene Laboratory also performed the inter- 
pretation of the rubber cracking determinations. 


3. Results 


Nitrogen dioxide levels were found to be 
below the limit of detection (0.1 pphm) of the 
method used 94 percent of the time. Only a 
single value above 0.5 part per a hundred milli- 
on (pphm) was measured. Only five of the 32 
carbon monoxide measurements (16 percent) 
showed any detectable reaction, and these were 
estimated to correspond to values below 1 part 
per million (ppm). Though determinations of 
carbon monoxide in this concentration range 
are highly inaccurate, it is estimated that con- 
centrations above 0.2 ppm can be detected. 
Detectable amounts of sulfur dioxide (>0.01 
pphm) were found in 42 percent of the samples, 
though 91 percent lay below 0.5 pphm. Seven 
values were above 1 pphm, with a maximum 
value of 4.6 pphm. The occasional high values 
of the first two gases are almost certainly trace- 
able to contamination from the ships. It is less 
certain that the sulfur dioxide levels measured 
are entirely caused by local sources. JUNGE (1957) 
measured concentrations of sulfur dioxide of 
the order of 0.03 pphm in marine air; the bulk 
of the non-zero values lay in this range. 

The first year’s oxidant data were rejected 
when it was shown that they were confounded 
by photolysis of the reagent during the sampling 
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Fig. 1. Cumulative frequency distributions of oxidant in 
mid-Pacific by two methods, August 1958. 


period. Data for the second year are summarized 
in Figure 1. Roughly the expected ratio was 
found between the values given by the Brewer- 
Mast recorder, a potassium iodide instrument, 
and the phenolphthalin method (RENZETTI and 
ROMANOVSKY, 1956). While no attempt has 
been made to prove the logarithmico-normal 
distributions suggested by Figure 1, it is clear 
that the two measures give very similar distri- 
butions of values. It will also be noted that no 
concentrations were obtained in excess of 7 
pphm by phenolphthalin. The absence of any 
higher values, while not conclusive, makes it 
seem likely that the high levels measured on off- 
shore islands by Harrison and LODGE (1959) 
were caused by urban sources, or at least by the 
influence of land. 

Rubber cracking determinations were incon- 
clusive. There was no clear correlation between 
measured crack rates and ozone level as measured 
by either technique. It appears that the uncer- 
tainty in the cracking measurements exceeds 
the total span of measured ozone levels. The 
mean rate of cracking for samples exposed for 
24 hours was 0.245 mm hr.-1. 

The particulate analyses are summarized in 
Figure 2. Separate plots of the two years’ data 
showed them to have similar dispersion, so they 
have been lumped. Nitrates were measured only 
on the 1958 samples. Several points are note- 
worthy. Nitrates were in general an order of 
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magnitude less than chloride, in accordance 
with the findings of JUNGE (1956). Rather more 
surprising is the fact that sulfate levels tend to 
be higher than chloride, in contrast to the 
roughly seven-fold excess of chloride over sul- 
fate in sea water. ERIKSSON (1959) has noted that 
the total cyclic sulfate must be of the same order 
of magnitude as the chloride, and suggests that 
the ’’excess” sulfur may escape from the sea as 
gaseous hydrogen sulfide. While these data give 
no direct evidence in this matter, they should 
substantially exclude urban sources as major 
contributors to the overall sulfate levels meas- 
ured. 

Benzene soluble organic compounds were 
determined on a few samples. The mean level 
was approximately 1.6 ug/m®. A rough calcu- 
lation shows that this loading could be account- 
ed for by assuming a monolayer of organic 
matter on the surface of the sea at the point of 
formation of the aerosol. 


4. Conclusions 


All gases measured, except ozone and possibly 
sulfur dioxide, have in nature normal levels 
below the limit of detection by the present 
instrumentation. Ozone levels in mid-ocean are 
variable; during the period under study, they 
ranged from 1 to 7 pphm by phenolphthalin or 
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Fig. 2. Cumulative frequency distributions of particulate 
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from 0.25 to somewhat over 3 pphm by potassi- 
um iodide. No regular diurnal variation was 
apparent. 

In the particulate state, sulfate levels were 
somewhat higher than chloride, and nitrate was 
an order of magnitude less. Organics were low, 
at a concentration which could be accounted 
for by the pickup of a monolayer on sea-spray 
droplets. 

A reasonably good measure of natural levels 
of pollutants” was obtained for the summer 
season in mid-Pacific. Further work will be 
necessary to determine whether there is any 
marked seasonal variation. 
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Abstract 


Artificial radioactivity in rainwater was measured during the passage of a cold front over 
The Netherlands. The specific radioactivity of the frontal rain increased considerably after the 
frontal veering of the wind. Concurrently the “age” of the fission products decreased from 45 
to 17 days, indicating a supply of atomic debris from recent nuclear explosions in arctic regions. 
The specific radioactivity of the air increased considerably on the day of the frontal passage. 
Trajectories of pre-frontal and post-frontal air indicate that the pre-frontal air comes from a 
source region in the Mediterranean and the post-frontal air from the arctic regions. The diagram 


of the frontal passage is analysed. 


Introduction 


In a recent paper (BLEICHRODT, BLOK, 
DexkKER and Lock, 1959) data showing a rela- 
tion between the amount of artificial radio- 
activity in rain and rain intensity have been 
presented. As a rule the specific radioactivity 
of rain was high during periods of low rain 
intensity and conversely. It is of importance to 
mention that only showers or showery pre- 
cipitation occurring in unstable air masses were 
described. 

The research reported in the present paper 
refers to the changes in specific artificial radio- 
activity of rain and air during the passage of a 
front over The Netherlands. 


The meteorological situation 


On the morning of 31 October 1958, a cold 
front approached the Netherlands coast from 
the west. This front separated two distinct air 
masses, namely modified tropical air in the east 
and relatively fresh arctic or polar air in the 


west. The winds in the tropical air were weak; 
the air mass was rather stable with a high 
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Synoptic surface map of 31 October 1958, 
12.00 GMT. 
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3500 mb 
31 Oct 1958 


Fig. 2. Isochrones of cold front between 00.00 GMT 
and 18.00 GMT (31 October 1958); positions of 
stations. 


humidity in the lower layers (fog over large 
areas). The wind speed was much higher in 
the polar air, which also was rather unstable 
and in which especially over the sea shower 
activity was observed. 

Figure 1 shows the synoptic surface map of 
12.00 GMT. 

Figure 2 gives the successive positions of the 
cold front over The Netherlands and adjacent 
areas between 06.00 GMT and 18.00 GMT. 
It appears that a wave developed in the cold 
front. The top of the wave moved from SW 
to NE over the middle of the country; the 


Fig. 3. Absolute topography of s00 mb level of 31 
October 1958, 12.00 GMT. 


heaviest rains (20 mm or more) were observed 
just NW of the wave top. 

Figure 3 demonstrates the synoptic situa- 
tion at the soo mb level, again for 31 October 
12.00 GMT. It appears that the surface front 
is followed by a deep post-frontal trough in 
the cold air, the successive positions of which 
are also indicated. 

The differences in structure between the 
tropical and polar air can best be appreciated 
from table I which gives the temperatures at 
various levels over De Bilt for 31 October 
00.00 and 12.00 GMT and for 1 November 
00.00 GMT. 

In figure 4 the recordings and observations 
of various meteorological elements during the 
frontal passage, at the stations Ypenburg! and 
Rijswijk! are reproduced. Visibility (a), relative 
humidity (b) and temperature (c) were observ- 
ed at regular intervals at Ypenburg; pressure 


1 The position of the various stations mentioned in 
this paper is given in figure 2. 


Table 1. Upper air temperatures (°C). 


Pressure 
Be 1000 | 850 700 500 400 300 200 100 Tropopause 
Time 
31 Oct. 00.00 GMT 6 2 6 21 32 47 58 60 60 (230 mb) 
31 Oct. 12.00 GMT 8 2 | — 6| —22| — 34 40 52 54 55 (215 mb) 
1 Nov. 00.00 GMT 6 2 12 30 40 47 51 57| — 48 (330 mb) 
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Fig. 4. Diagram of meteorological elements and radio- 

activity observations at Ypenburg and Rijswijk (31 

October 1958). a. visibility; b. relative humidity; c. 

temperature; d. pressure; e. wind speed; f. wind direc- 

tion; g. rain intensity; h. radioactivity; i. “age” of 
fission products. 


(d), wind speed (e) and wind direction (f) 
were recorded at the same station. The rain 
intensity (g) was observed at Rijswijk. Ypen- 
burg and Rijswijk lie at a distance of 1,600 m. 
The front arrived at Rijswijk approximately 3 
minutes earlier than at Ypenburg. This time 
difference may be considered negligible in the 
interpretation of the data. 

The frontal passage is indicated by a veering 
of the wind from southwest to northnorthwest, 
which took place between 11.00 and 11.30 
GMT. During this period high rain intensities 
were observed. Immediately after the veering 
the wind speed increased considerably. 


Radioactivity in rainwater 


Rainwater samples were taken at Rijswijk 
between 8.20 and 15.20 GMT. The rainwater 
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was caught in a galvanized tray of 5 m? which 
was placed with a slight slope on the flat roof 
of the Medical Biological Laboratory. A meas- 
uring glass was placed under an outlet-pipe 
at the lowest point of the tray. Every time 
when about half a liter of rainwater had been 
collected in the glass, the water film in the tray 
was quickly swept down with a screen wiper. 
The total volume of the samples, usually about 
1 liter, and the time during which they were 
obtained were listed. From these data average 
rain intensities for relatively short periods could 
be computed. The results are given in figure 
4 (g), in which the dots indicate the time of the 
middle of the sampling period. 

Before sampling of rainwater was started 
the tray had been thoroughly cleaned with 
tap water in order to remove previously 
deposited fall-out material. 

The artificial radioactivity of the rainwater 
samples was determined as follows. The bulk 
of the beta activity was concentrated by stirring 
the water with 4 g of active coal (Norit C, 
Norit Verkoopmaatschappij, Amsterdam), 
which was then filtered off. The filtrate was 
evaporated to dryness and the residue together 
with the coal containing filter ashed at 500° C. 
The ash was weighed and a fixed amount was 
put in a brass sample cup. 

Seven days after collecting the rainwater, the 
beta radioactivity of the samples was measured 
(standard error < 2.5%,) using an end-window 
Geiger-Miiller tube (Philips 18506) with a 
diameter of 3 cm and a window thickness of 
3 mg/cm”. The beta radioactivity per liter of 
rainwater was then computed from the measur- 
ed activity, the sample weight, the total ash 
weight and the volume of water collected. The 
result is shown in figure 4 (h) (dotted radio- 
activity curve). One radioactivity unit corre- 
sponds with about 5 pc/liter. 

The beta radioactivity of nine ash portions 
was measured in duplicate to investigate the 
possible influence on the changes in specific 
radioactivity of highly radioactive particles 
which can be found shortly after nuclear ex- 
plosions. As the differences between the dupli- 
cate measurements proved to be smaller than 
10%, it seems justified to exclude this in- 
fluence. 

The beta radioactivity of ten samples was 
measured periodically during three months. 
From the decay rate the “age” of the fission 
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Fig. 5. Decay curves of debris collected at 9.10 (A), 
11.10 (0) and 14.40 GMT (@). 


products was determined, using the approxi- 
mate relation: 


A Ar (1) 


where A is the activity at time t and A, the 
activity at £= 1. To this end the values of 
Ai? were plotted versus f and the best 
fitting curve was extrapolated to A-1/12 = 0. 
In this way an indication was obtained of the 
moment when the artificial radioactivity was 
introduced into the atmosphere. 

The decay curves often consist of a combina- 
tion of two straight line segments with a 
slightly bent part in between. In the beginning 
the decay is determined mainly by the rela- 
tively fresh fission products. After a few weeks 
the influence of the older debris becomes 
noticeable. As example three decay curves are 
shown in figure 5. 

The result of the “age” determinations is 
plotted in figure 4 (i). 

With the aid of the obtained “age”-curve 
it is possible to calculate the values of the 
specific radioactivity of the rainwater at the 
date of sampling. The drawn curve of figure 
4 (h) gives the so obtained values of the specific 
radioactivity at sampling time. 

The radioactivity data show the following 
features. It is obvious that the specific radio- 
activity of the rain was generally low before 
the frontal passage with the exception of an 
isolated peak at about 9.10 GMT. There is a 
tendency of the radioactivity to decrease with 
the approaching of the front. During the 


veering of the wind the activity increases about 
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so % and after some variations, in which it 
occasionally drops to the pre-frontal level, a 
marked increase to a value approximately 5 
times the pre-frontal value starts at 12.45 GMT. 
It is furthermore noted that neither in the pre- 
frontal rain nor in the post-frontal rain there 
is a clear indication of an inverse relationship 
between rain intensity and radioactivity. 

The “age’’-curve of figure 4 (i) shows that 
the “age” of the fission products in the pre- 
frontal rain was about 45 days and about 17 
days in the post-frontal rain. The pre-frontal 
as well as the post-frontal “age” values are 
sufficiently uniform and the difference between 
the two is sufficiently large to warrant the 
conclusion that the nature of the artificial radio- 
active material changed completely during the 
frontal passage. 

It should be noted here that the above-men- 
tioned high pre-frontal specific radioactivity 
which was found at9.10 GMT will not be related 
to the high post-frontal values, as exactly for 
this sample an “age’’ determination has been 
made (43 days); this pre-frontal peak must 
therefore be due to unknown processes. The 
presence of a very active particle must be 
considered improbable, as a duplicate measure- 
ment showed the same high value of the ac- 


tivity. 
Airborne radioactivity 


Airborne radioactivity was sampled at 
ground level at the stations Den Helder, De 
Bilt, Rijswijk and Eindhoven (for positions 
cf. figure 2). 

By means of a vacuum pump 30 to 50 m3 
(measured with a gasmeter) of air is drawn 
through a membrane filter (Membranfilter 
Gesellschaft, Göttingen, Germany) in 24 hours. 
The effective surface of the filter is 5.2 cm?. 
The filters are changed daily at 8.00 GMT. 
In order to avoid loss of radioactive particles 
the dust is fixed to the filter paper by treating 
the filter with a solution of formvar in 1,2- 
dichloro-ethane (20 g per liter). 

The beta activity of the filters was measured 
at the Medical Biological Laboratory by means 
of an end-window Geiger-Müller tube (Philips 
18506) three days after sampling to allow the 
natural radioactivity to decay. 

In table 2 the airborne artificial radioactivity 
as measured during the period October 26th— 
November sth 1958 is listed. 
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Table 2. Airborne artificial radioactivity in pc/m?. 


Stations 


Sampling period | Den Rijs- |Eind-|Aver- 


e be 
who lgass 
der as (Z.H.)| ven 

October 26/27 ...| 1.02| 0.47| 0.94] 0.66| 0.77 
» 27/28...| 1.15| I-06} 1.13] 0.58] 0.98 

» 28/29...| 1.15] 0.95} 0.92] 0.91] 0.98 

» 29/30...| 0.74] 0.81] 0.54] 0.52] 0.65 

» 30/37 2.1, 1.291.,7.361.7.481..1.42|31.39 

» 31/Nov.1| 2.98] 2.59} 3.65] 2.79] 3.00 
November 1/2...| 2.57] 2.11; — | 2.16} 2.28 
» 2/2, .| 22.08] "2:40|* — || 71.00] 2.12 

» 3/431) 5.44} 05-09) 223.7] 95:43} *5.35 

» 45 AA Lie 2a 52:07 BE99 297 


1 No data available. 

2 High value probably due to one or more very 
active particles. 

3 Value of Rijswijk (Z. H.) excluded. 


It appears that the radioactivity was rela- 
tively low until the morning of 30 October. 
The samples collected on the morning of the 
31st show an increase of 0.7 pc/m? in the 
average. The “age” of the samples collected on 
31 October was not determined. It is however 
highly probable that the radioactivity in the 
pre-frontal tropical air was the same as that in 
the pre-frontal rain. In that case the rise on 
31 October must be due to processes in the 
tropical air, of which pre-frontal subsidence or 
changes in the turbulence conditions could be 
mentioned. 

A strong increase of 1.6 pc/m? in the average 
is observed on ı November, i.c. after the 
frontal passage. The high activity of the post- 
frontal rain and the high activity of the post- 
frontal polar air go together. The highest 
increase is observed at Rijswijk. The “age” of 
the airborne radioactivity collected during 31 
October—1 November at Rijswijk was 20.5 
days. This agrees fairly well with that of the 
post-frontal rainwater. 

It is well-known that the specific artificial 
radioactivity of the air usually increases with 
height. The question may be asked whether 
the radioactivity at the surface increased so 
strongly solely as a result of a more intense 
turbulent diffusion due to a higher lapse rate, 
higher wind speeds and greater vertical wind- 
shear, or whether the radioactivity on the 
whole was much higher in the polar air than 
in the tropical air. In this connection it is 
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important to note that the specific radioactivity 
of the surface air increased by a factor 2 
approximately and that of the rain by a 
factor 5. 

Now the process of rain formation takes 
place in a thick layer and the specific radio- 
activity of the rain can, therefore, be considered 
more representative for the air mass as a 
whole than the specific radioactivity of surface 
air. Consequently it may be concluded that 
the radioactivity in the tropospheric layers 
above the planetary friction layer in the polar 
air mass was approximately 5 times higher than 
in the tropical air, but that in spite of the 
greater turbulent diffusion the downward 
transport of radioactive matter was not strong 
enough to bring about also a fivefold rise in 
the airborne radioactivity at ground level. This 
seems to indicate that high values of radio- 
activity in the polar air were present in the 
middle of the troposphere and above. 

A new rise in the airborne radioactivity was 
observed during the period 3—4 November. A 
front with a new outbreak of arctic air in 
which the wind veered to a northerly direction 
passed The Netherlands in the afternoon of 
3 November 1958. The high value measured 
at Rijswijk is probably due to one or more 
very active particles. 


Air trajectories 


The explanation of the change in the radio- 
activity regime on 31 October is easily found 
when the trajectories of the pre-frontal and 
post-frontal air are investigated. Under the 
assumption of horizontal motion trajectories 
have been computed for an air mass present 
over The Netherlands on 31 October 00.00 
GMT (representing the tropical air) and an air 
mass which covered The Netherlands on 1 
November 00.00 GMT (representing the polar 
air). Fig. 6, 7, 8 and 9 show these trajectories 
for the 900 mb level (based on gradient winds 
taken from the surface map) and for the 700, 
500, and 300 mb levels respectively. 

It appears that up to s00 mb the pre- 
frontal air originated in the Mediterranean 
area. The wind speed decreased with increasing 
altitude, as usually is the case in a predomi- 
nantly easterly current. The air parcels at 300 
mb originated from the North Atlantic. 

The post-frontal air came from the arctic 
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Fig. 6. Air trajectories at 900 mb level. 
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Fig. 7. Air trajectories at 700 mb level. 


regions in all levels. The wind speeds increased 
rapidly with altitude, so that a relatively strong 
vertical windshear existed. 

According to information published by the 
United States Atomic Energy Commission 
(Harpy, E. P. and KLEIN, S., 1959) nuclear 
test explosions were performed north of the 
arctic circle on Io, 12, 15, 18, 19, 20, 22, 24 
and 25 October. The radioactivity found in 
the polar air is undoubtedly due to one or more 
of these test explosions. Many more data of 
radioactivity measurements would be necessary 
in order to determine at which levels the 
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Fig. 9. Air trajectories at 300 mb level. 


transport of the radioactive material took place. 

It is impossible to state where the atomic 
debris present in the tropical air was injected 
into the atmosphere. 


Further analysis of the data collected during 
the frontal passage 


Figure 4 raises a number of problems related 
to cloud physics and to the structure of fronts. 
For the time being the authors limit themselves 
to point out that apparently $ phases can 
clearly be distinguished in the diagram. 

Phase a, from the beginning to 11.00 GMT. 
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Pre-frontal tropical air with uniform tem- 
perature, uniform humidity and low uniform 
visibility is present over the station. The rain 
intensity varies, there is a light pressure fall and 
slowly veering wind of low speed. The radio- 
activity decreases slowly; the “age” of the 
radioactive material is constant (approximately 
47 days). 

Phase b, from 11.00 to 12.00 GMT. This 
phase is marked by a strong veering of the 
wind at the beginning of the period and of a 
considerable stepwise increase of the wind 
speed at the end of the period. The rain 
intensity increases to a high value in the be- 
ginning of the period and decreases to its 
original value at the end. The visibility increases 
considerably. At the end of the period the 
temperature is slightly higher than at the be- 
ginning and the humidity has decreased. There 
is still a slight pressure fall. The radioactivity 
increases about 50%; the “age” indicates the 
presence of a mixture of old and fresh debris. 

Phase c, from 12.00 to 12.45 GMT. This is a 
phase, which is characterized by an approxi- 
mately constant radioactivity level and a low 
“age” (18 days). Wind direction and speed are 
constant and the rain intensity shows a peak. 
At the beginning of the period the pressure 
starts to rise. The visibility is high. 
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Phase d, from 12.45 to 14.40 GMT. This 
phase is marked by a strong increase of radio- 
activity of low “age”. The wind direction is 
constant. The wind speed is variable but drops 
to a lower value at the end of the period. 
Visibility is high, temperature and humidity 
have slightly decreased and the pressure rises. 

Phase e, from 14.40—end. The radioactivity 
is at a high level. The wind direction and 
speed are constant. The pressure rises further 
while temperature falls and humidity rises 
apparently under influence of the cooling with 
the approaching setting of the sun. 
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Abstract 


Measurements of the radio-active fall-out occurring in weekly rainfall samples at Melbourne 
(Aspendale) are presented for the period June 1958 to February 1959. Radio-active decay 
analyses show the activity to have arisen from nuclear weapons tests carried out near the equator 
during the period April to September 1958. The results indicate that, despite presumably con- 
tinuous dispersion, individual radio-active clouds maintain their discreteness and a limited E—W 
extent even after six or nine months residence in the stratosphere. This property allows a pole- 
ward drift speed of order 1 m/sec and zonal circulation velocities ranging from 6 to 17 m/sec 


to be inferred. 


I. Introduction 


The global fall-out of radio-activity as a 
result of nuclear weapons tests has been the 
subject of much enquiry in recent years. The 
variation of the fall-out pattern as a function 
of both time and latitude has posed a number 
of problems for the meteorologist. Several fall- 
out models have been suggested (Lipsy 1958, 
MACHTA 1958, MARTELL 1959, STEWART et al 
1957) and the difficulty of deciding whichmay 
be the correct one is an indication of the limited 
extent of the present knowledge of the stratos- 
pheric circulation. 

The investigation to be reported here began 
chiefly as a contribution to the International 
Geophysical Year, but the possibility of using 
radio-activity as a tracer of atmospheric motions 
was also considered. Measurements of the radio- 
activity of rain-water were commenced in 
June 1958, and are still continuing. It is felt 
that a detailed analysis of the results from June 
1958 to February 1959 reveal features of the 
high-level circulations which bear directly on 
considerations of the various fall-out models 
being made by other workers. 
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2. Techniques of Measurement 


The most widely used method of measuring 
the radio-activity of rainwater is to evaporate 
the sample into a small planchette and de- 
termine the total f-activity using a thin- 
window counter. This process is extremely 
time consuming in the evaporation stage and 
an ion-exchange technique of the type sug- 
gested by Epvarson (1957) seemed more ap- 
propriate to our problem. 

The ion-exchange resin (Bio-Deminrolit 
FF)! was used in two ways. The normal method, 
called in the following, “regular”, consisted 
of collecting the rain with plastic-coated funnels 
and polythene jars and passing it twice through 
a single column of 10 cc of resin at a flow rate 
of approximately 1 litre/hr. The resin was 
then ashed in a silica crucible and the activity 
determined using an end-window ß-counter 
(Phillips 18514). Tests indicated that 90% of 
the activity was taken up in a single passage of 
the rain-water through the resin. 

The alternative method, called “auto”, is 


1 Supplied by the Permutit Co. of Australia, Sydney, 
Australia. 
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Fig. 1. Radio-activity of weekly samples of rain water. 


particularly suited to sampling at a large 
number of distant field stations in that the rain- 
water passes through the resin at time of col- 
lection. A plastic cartridge containing 10 cc of 
the resin is fitted in an S-tube between the 
funnel and the collecting vessel, the S-tube 
being suitably filled with a measured quantity 
of distilled water to maintain the resin ef- 
ficiency. 

The equipment was calibrated by using 
standard solutions of potassium carbonate. A 
rain-sampling period of one week was used, 
except for the first case (June 1958) in which 
a monthly sample was taken. 

For the few weeks when no rain occurred 
the “regular” collector was washed out with 
one litre of distilled water and this was then 
treated in the usual way. Before each sampling 
period the collecting funnels and jars were 
thoroughly washed with distilled water. 


3. Estimation of Explosion Dates 


An estimate of the date of origin of the 
fission sample can be made by using the decay 
law Le I, Mia of) Hunrentsand’ Barrow 
(1951), in which the activity I; at time t after 
the explosion is related to the activity I, at 
t =1 (In our case t is measured in weeks). 
Using this expression one can write L, = 0.41 1, 
and the unknown (f) can thus be estimated 
by determining the time for the activity of a 
sample to decay to 0.41 times its initial ob- 
served value. 


If the decay of a sample is plotted on a log- 
log scale of activity versus time, linear inter- 
polation can be used to determine f by the 
above method even if the time origin of the 
graph is only approximately correct. In a few 
cases a non-linear plot was obtained, and the 
estimated explosion dates for these are there- 
fore somewhat uncertain. 

It is well known that the exponent in the 
decay law assumes slightly different values for 
different types of nuclear explosions, and the 
dangers of using a constant value of 1.3 are 
well recognised. Apart from this point, which 
will be discussed later, the error in estimating 
t is thought to be about 6 %. 


4. Results 


The results obtained are listed in Table ı 
and presented in Figure 1. 

Two things of interest are immediately obvio- 
us from Table 1. Firstly, the good agreement of 
the “auto” and “regular” techniques in estima- 
ting the explosion date is encouraging. Second- 
ly, it is evident that the estimated explosion date 
varies considerably from week to week. 

Series of United States and United King- 
dom test explosions were conducted in equa- 
torial regions from April to September 1958. 
If the debris from such series of explosions 
were distributed quickly throughout the whole 
atmosphere so that a homogeneous mixture 
was created it would most likely follow a decay 
law appropriate to a mean of the explosion 
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Table 1. Radio-active rain measurements at Melbourne, Australia. 
ee eles ee Re ie Un 


Date off mount mc} Estimated explosion dates 
Collection | RR | curies/litre | Km? 
mm Regular | Auto Series 
a a Ne eg Sl ee i u 
June 1958| 21 23.6 0.49 = | 
7. 7.58 22 3 0.16 — 
14. 7.58 38 122 0.47 os 
21. 7.58 21 164 3.52 | May 25 A 
28. 7.58 07 295 2.20 July 3 B 
4. 8.58 Io 193 1.99 | June 26 B 
II. 8.58 15 92 1.4I | June 3 A 
18. 8.58 32 147 4.70 | June Io June 11 C 
25. 8.58 28 63 2.70 | June 24 June 22 B 
I. 9.58 09 86 0.74 | June 28 June 17 B 
8. 9.58 04 334 1.25 | June 16 — C 
15. 9.58 14 78 1.05 | May 19 — A 
22. 9.58 42 41 1.71 | May 25 May 25 A 
29. 9.58 07 246 1.68 | May 29 == A 
6.10.58 19 56 RO July 2 — B 
13.10.58 23 96 2.15 | July 13 Aug. 4 D 
20.10.58 13 Teer 1.48 | Aug. 14 — E 
27.10.58 32 19.4 0.62 | July 15 July 25 D 
3.11.58 —- — | No measurement 
10.11.58 25 — | No measurement 
17.11.58 80 2.034 PAUSST2 Aug. II E 
24.11.58 20 60 1.18 | June 25 July 27 B 
1.12.59 20 76 1.49 | Aug. 31 — F 
8.12.58 00 No rain 0.44 | Sept. 8 Aug. 25 F 
15.12.58 39 273 0.67 | June 23 — B 
22.12.58 02 325 0.49 | Sept. 4 — F 
u 00 —- — | No measurement 
5. 1.59 No rain 0.59 | Sept. 11 — F 
12. 1.59 OI No rain O18 | July 27 — G 
19. 1.59 00 No rain 0:35. | Une TE — (® 
26. 1.59 u — | No measurement 
2. ++ 06 89 0.58 | July 27 — G 
9. 2.59 04 152 0.57 | July 30 — G 
16. 2.59 19 56 1.05 | Aug. 4 — G 
23. 2.59 00 No rain 0.28 | Sept. 16 — F 


dates and with an exponent approximately 
equal to 1.3. Using an incorrect exponent, the 
estimated explosion dates would progress slowly 
from week to week. This picture would only 
emerge from the present data if the correct 
exponent varied erratically from week to week. 
This in itself implies a degree of inhomogeneity 
inconsistent with the hypothesis of a uniform 
distribution in the atmosphere. 

What appears to be a more plausible explana- 
tion of the variability of the present results 
is that the use of the #13 law is substantially 
correct and that samples from different test 
explosions are obtained in consecutive weeks. 
If this view is adopted for the time being, the 
estimated explosion dates can be used to Sete 
the passage of a particular atomic cloud. This 
has been done in Table 1 in which the re-oc- 
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currence of the same (or nearly the same) date 
is recognised as being due to a sample from 
the same test explosion (Series A, B etc.). The 
mean explosion date for each series can then 
be compared (Table 2) with the dates of 
known nuclear tests obtained from press reports. 

Other United Kingdom tests during this 
period were low yield explosions on August 
22 and September 23, and a megaton explosion 
on April 28. The United States also conducted 
tests on May 11—12 and July 5. No information 
of the yield of any of the U.S.A. tests is avail- 
able, although from micro-barograph records 
here it is believed that large explosions took 
place on June 28 and July 12. 

It is of interest to analyse the arrival time and 
circulation times of each of the atomic clouds 


(Table 3). 
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Table 2. Comparison of estimated and known explosion dates. 


(ee EE eer 


i DRE ea Actual Explosion Dates 
Be Explosion Dates a L' 
Es | 
A May 26 May 25 U.S.A. Test Eniwetok 
B June 28 June 27 U.S.A. Test Eniwetok 
June 27 U.S.A. Test Bikini 
June 28 U.S.A. Test Eniwetok 
July 2 U.S.A. Test Bikini 
C June 12 June 10 U.S.A. Test Bikini 
June 14 U.S.A. Test Bikini 
D July 14 July 11 U.S.A. Test Bikini ae 
E Aug. 13 Aug. 12 U.S.A. Test Johnston Island (Missile) 
F Sept. 6 Sept U. K. Test Christmas Island, High Yield 
Sept. 9 U. K. Test Christmas Island, High Yield 
G July 30 July 26 U.S.A. Test Eniwetok 
Aug. I U.S.A. Test Johnston Island (Missile) 
Table 3. Analysis of arrival and circulation times. 
; Arrival time ist Circulation 2nd Circulation 3rd Circulation 
a (weeks) time (weeks) time (weeks) time (weeks) 
A 7:5 3 5,6,7 ve 
B 3, 4 3, 4 5 7, 10 
C 9 3 (19) = 
D 13 2 == er 
E 9 4 a as 
F 11.5 3 7 
G 19 3, 4, 5 — — 
Mean times 10.6 3.1 6.0 | 8.5 
Corrected 9.0 
Zn a Mean zonal circulation 
N 17 m/sec 10 m/sec 6 m/sec 
I.I m/sec. 


The various values obtained are grouped 
closely about their mean values. A mean drift 
velocity and mean zonal circulations can be 
calculated from the mean arrival time and mean 
circulation times respectively (Table 3). As the 
first arrival of a cloud at this latitude is un- 
likely to occur precisely at the longitude of 
Melbourne, the arrival time will be over- 
estimated on the average by approximately 
half the zonal circulation time. An appropriate 
correction is made in Table 3. 


5. Discussion 


A tentative explanation of the above results 
can be advanced in the following terms. The 
observed fact that fission debris arrived at 


ground level at latitude 38° S for a period of 
several months after a series of test explosions 
clearly implies a stratospheric circulation from 
equator to middle latitudes if one accepts the 
well established residence time of one month 
for radio-active particles within the tropos- 
phere. The radio-active cloud produced by an 
atomic bomb appears to migrate (or disperse) 
polewards from the equator with a limited 
amount of spreading in an east-west direction. 
The latitude dependance of Sr9 deposition 
observed by several authors, Stewart et al 
(1957), MARTELL (1959) would favour the idea 
of migration rather than dispersion. The 
greatest amount of east-west spreading ob- 
served was that of Series A and Series G which 
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were detected for three consecutive weeks at a 
period of four to seven months after the date 
of explosion. This indicates an east-west spread 
of about half the global circumference at this 
latitude. 

Mintz (1954) has presented data on mean 
zonal winds as a function of height, latitude 
and season. The value of 15 m/sec for the first 
observed circulation of the earth by the cloud 
suggests from Mintz’s data a height of approxi- 
mately 100 mb. (16 km). Values of 7.5 m/sec 
and s m/sec for the subsequent circulations 
imply that even greater heights are involved. 
Heights of 400 mb (7 km) to 700 mb (3 km) 
would also give the correct zonal velocities but 
these lie well within the troposphere and can 
presumably be ignored. 

WEXLER (1957) in an examination of the 
Krakatoa dust cloud movement finds evidence 
for a mean meridional drift at 19 km height of 
1.6 m/sec at 10° N decreasing to 0.6 m/sec at 
40° N. Similar values are reported in another 
context by PRIESTLEY and Troup (1954) for 
the region between the latitudes 30° and 35°. 
The value of 1.1 m/sec found in the present 
investigation for a mean meridional drift 
between equatorial regions and 38°S seems 
therefore to be consistent with the above 
notion of stratospheric circulation. 

A number of models of global fall-out have 
been suggested; Lippy (1958) proposes that 
stratospheric debris is distributed uniformly 
over both hemispheres with a stratospheric 
residence time of $—10 years, and attributes the 
observed latitude dependence and increase in 
spring fall-out to rainfall variations and the 
testing of sub-megation weapons in middle 
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latitudes. STEwART et al (1957) find their 
results in agreement with the zonal circula- 
tion model proposed by Brewer (1949) and 
Dosson (1956) in which air rising into the 
stratosphere at the equator descends into the 
troposphere in middle and polar latitudes. A 
recent discussion by MARTELL (1959) attributes 
the increase in spring fallout to Russian tests in 
the immediately preceding autumn together 
with a short stratospheric residence time in 
middle latitudes. 

The present results strongly suggest the 
existence of a poleward drift in the lower 
stratosphere, and indicate that even after a 
period of six to nine months the debris from 
a test explosion is not distributed uniformly 
throughout the stratosphere. During this period 
some of the radio-activity appears at the ground 
in middle latitudes, but there is no indica- 
tion of the mechanism by which the radio- 
activity reaches the troposphere. 

A plausible view would be that the stratos- 
pheric residence time depends on the size of the 
bomb and to some extent the height at which 
it was exploded. 

An examination of fall-out data for other 
parts of the world in the same terms would be 
of obvious value. A sampling period of one 
week appears to be particularly suited to this 
type of study. 
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Abstract 


A systematic variation with season and latitude in the concentration and isotopic abun- 
dance of atmospheric carbon dioxide has been found in the northern hemisphere. In Ant- 
arctica, however, a small but persistent increase in concentration has been found. Possible 


causes for these variations are discussed. 


The content of carbon dioxide in the atmos- 
phere, in contrast to oxygen, nitrogen, and the 
rare gases, has been found to be significantly 
variable (GLUECKAUE, 1951). New data indicate, 
however, that the degree of variability is 
smaller and the variations are more systematic 
than previously believed (STEPANOVA, 1952, 
SLOCUM, 1955, FONSELIUS ET AL., 1956, 
CALLENDAR, 1958, BRAY, 1959). 

Three gas analysers, as described by SmitH 
(1953), equipped with strip chart recorders, 
have been employed to measure the concentra- 
tion of carbon dioxide continuously at stations 
in Antarctica, Hawaii, and California. A fourth 
analyser has been used in the laboratory to 
analyse samples of air collected in glass flasks 
at various places. These analysers provide 
direct comparisons of the partial pressure of 
carbon dioxide in air relative to that in prepared 
mixtures of carbon dioxide in nitrogen gas. 
The concentration of carbon dioxide in these 
reference mixtures is determined manometri- 
cally (KEELING, 1958). 


Contribution from the Scripps Institution of Ocea- 
nography, New Series. This paper represents preliminary 
results of research carried out as part of the International 
Geophysical Year and its extension. Support was provided 
by the National Science Foundation, United States 

Weather Bureau, and the Office of Naval Research. 


The relative accuracy of the data presented 
here is approximately + 0.3 p.p.m. (parts per 
million by volume). The uncertainty in the 
absolute values is considerably larger, however, 
since only a preliminary calibration of the 
reference mixtures has been made. When an 
accurate manometric calibration is completed, 
the absolute accuracy is expected to approach 
+ 0.1 p.p.m. 

Monthly averages of the data from the con- 
tinuously recording stations and from collec- 
tions in flasks are presented numerically in 
Table 1 and plotted in figures 1 and 2. The 
locations of the sampling stations and tracks 
are shown in figure 3. 

Local contamination has been found to occur 
at all three continuous recording stations. At 
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Fig. 1. Variation in concentration of atmospheric carbon 
dioxide in the Northern Hemisphere. 
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Table 1. Monthly Average Concentrations of Atmospheric Carbon Dioxide at Various Stations 


in Parts per Million of Dry Air by Volume. 


Continuous Recording 


Surface Flask Samples on nae 


Stations Aircraft? 
Month Year : 5 Data at} Data at} Limit- 
pee Mauna | La | South ER Ron Lower | Higher | ing Air 
ur Loa Jollat | Pole en Air Air temp. 
temp: | temp. | °C.) 
Sept. 1957 SUITE | ES OC 
Oct. 310.9! 
Nov. 210 0NRS TO 
Dec. 312.02, 313407 311:6 
Jan. 1958 
Feb. 310.8 
Mar. STI.0 | 373,4 302.9 103155 
Apr. 311.4 | 314.4 | 314.3 314.01 927 
May SELF || SET 315.3 
Jun. 311.9 Se 314.91 | —2ı 
Jul. 311.8 | 3229°| 371.0 
Aug. 373.0, 17 372.390. 308.4 
Sept. 312.907 311.6 4) 308.7 373.3 308.31 | 309.6 —18 
Oct. 313.0 310.8 SI OMIS T2 —27 
Nov. 310.6 313.2 
Dec. 32:0 10314.3 
Jan. 1959 37257 |# 3104.09 312-8 Brae —36 
Feb. 313.5 | 314.4 | 313.4 aaa || 0 
Mar. BTA 003754 lignes, 
Apr. 314.7 | 315.0 | 314.3 316.61 | 316.6 —27 
May 315.3 | 315.0 | 314.3 317.3 4 
Jun. 315.2. 315.210 314.4 314.3 NS TA —21 
Jul. 313% 311.5 314-2 
Aug. BER. esos: x 314.2 308.1 310.7 —18 
Sept. SI. 3037. 0.3740 2 OMn Hama —27 
Oct, BLOW | sizer 313 0 3720 —27 
Nov. 311.8 323-7 313.9 3142 732302 — 36 
Dee: 312.5 313.5 315.6 | 314.3 —36 
Jan. 1960 513 4 (135467 317.6 | 315.6 —36 
Feb. BI32700815.4 315.6 —36 
Mar. 314.4 STAAT 


1 Data not shown in Figures I or 
2 Data have been separated into two groups based upon the temperature of the air at the point of sampling. 
Limiting air temperature refers to the temperature separating lower air temperature from higher. It has been chosen so 
that, for each period of sampling, the boundary between warmer and colder air lies at approximately 50° North (see 


figure 3). 


Little America, Antarctica, it was evidently due 
solely to combustion of fuel in the immediate 
vicinity of the station. It could be readily 
spotted from the significant fluctuations in the 
otherwise steady trace of the recorder pen and 
was eliminated from consideration in the initial 
reading of the charts. At Mauna Loa Observa- 
tory, Hawaii, a less prominent variability has 
been found in approximately half of the records. 
This is attributed to release of carbon dioxide 
by nearby volcanic vents; combustion on the 
island associated with agricultural, industrial, 
and domestic activities; and lower concentra- 
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6—002421 


(Hi Kae | 


© SOUTH POLE 
A DOWNWIND CRUISE 


© LITTLE AMERICA 


| 
= a 
es vel 


[eae 
gwuetes4567e9WNI2123456789 WN 1 23 
1957 1958 1959 1960 


MONTH AND YEAR 


Fig. 2. Variation in concentration of atmospheric carbon 
dioxide in the Southern Hemisphere. 
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tion of carbon dioxide in the air transported to 
the station by upslope winds. The values re- 
ported here are averages of data for periods of 
downslope winds or strong lateral winds when 
the concentration remained nearly constant 
for several hours or more. At La Jolla, Cali- 
fornia, the concentration has been found to be 
highly variable. Highest concentrations occur 
during light winds from the north, from the 
direction of Los Angeles; lowest concentrations 
when the wind is from the west or southwest 
and of moderate force or greater. Lowest 
weekly values usually do not differ by more 
than + I p.p.m. during any month, and, within 
a range of 2 p.p.m., agree with other data for 
the northern Pacific ocean. Monthly averages 
of these data, which presumably indicate nearly 
uncontaminated air, are cited here. 

Data for air collected in flasks from aircraft 
flying at 5 to 6 kilometers above the Pacific 
and Arctic oceans and from surface stations at 
the South Pole and on Arctic ice floes show a 
high degree of regularity. 


A clearly defined seasonal trend in concen- 
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Fig. 3. Location of stations and tracks for sampling of 
atmospheric carbon dioxide. Circles denote continuous 
recording stations. Triangles denote flask sampling sta- 
tions. Elevations in meters are given for locations more 
than 100 meters above sca level. The approximate mean 
position of the limiting temperature isotherm, used to 
separate the aircraft data into groups associated with 
higher and lower air temperature, is shown by a dotted 
line. 
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Fig. 4. Variation in concentration and C13/C! ratio of 
atmospheric carbon dioxide, based on Keeling (1961) for 
various locations near the Pacific coast of the United States. 


tration is found at all locations in the northern 
hemisphere. Going from south to north, the 
annual range of concentration becomes greater 
and the month of minimum concentration 
occurs earlier. Separating the samples from 
aircraft into two groups based upon the tem- 
perature of the air at the point of sampling (see 
table 1), the values for the high latitude, or 
polar, air are seen to have a greater seasonal 
range than the comparable values for air of the 
temperate zone. In contrast, data for the south- 
ern hemisphere do not indicate any seasonal 
variation. Data from Downwind Cruise in 
November and December, 1957, suggest that 
the concentrations observed over Antarctica 
prevail at all southern latitudes of the Pacific 
ocean. 

A seasonal variation in the isotopic abun- 
dance of carbon 13 in atmospheric carbon 
dioxide may also exist (KEELING, 1961). In 
figure 4 are plotted the observed seasonal 
variations in isotopic abundance and concen- 
tration, based upon samples collected at various 
stations near the Pacific coast of the United 
States during 1955 and 1956. Changes in relative 
abundance of carbon 13 with concentration 
obey a relationship found for air near plants 
(KEELING, 1958). The isotopic composition of 
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the carbon dioxide associated with the change 
in concentration, according to this relationship, 
is approximately —22 per mil, in good agree- 
ment with the value found for air of forests 
(KEELING, 1958) and for the carbon of plants 
growing on land (CRAIG, 1954, WICKMAN, 
1952). These data, therefore, indicate that the 
seasonal trend in concentration observed in the 
northern hemisphere is the result of the activity 
of land plants. This interpretation receives 
further support from the fact that maximum 
concentrations have been found to occur in 
spring at the outset of the summer growing 
season for plants in the temperate zone; that 
minimum concentrations occur in the fall, 
approximately at the end of the growing season. 
The observed absence of a seasonal trend in the 
southern hemisphere is then to be explained by 
the smaller area of growing plants found in 
the southern hemisphere at temperate and polar 
latitudes. These conclusions should be consid- 
ered tentative, however, since they are based on 
relatively few isotopic measurements. 

Where data extend beyond one year, aver- 
ages for the second year are higher than for the 
first year. At the South Pole, where the longest 
record exists, the concentration has increased 
at the rate of about 1.3 p.p.m. per year. Over 
the northern Pacific ocean the increase appears 
to be between 0.5 and 1.2 p.p.m. per year. 
Since measurements are still in progress, more 
reliable estimates of annual increase should be 
available in the future. At the South Pole the 
observed rate of increase is nearly that to be 
expected from the combustion of fossil fuel 


(1.4 p.p.m.), if no removal from the atmosphere 
takes place (REVELLE and Suess, 1957). From 
this agreement, one might be led to conclude 
that the oceans have been without effect in 
reducing the annual increase in concentration 
resulting from the combustion of fossil fuel. 
Since the seasonal variation in concentration 
observed in the northern hemisphere is several 
times larger than the annual increase, it is as 
reasonable to suppose, however, that a small 
change in the factors producing this seasonal 
variation may also have produced an annual 
change counteracting an oceanic effect. 
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The Influence of Carbon Dioxide Variations on the Atmospheric 
Heat Balance 


By LEWIS D. KAPLAN, Massachusetts Institute of Technology 


(Manuscript received August 20, 1959) 


Abstract 


Net fluxes of radiation in the 15 micron carbon dioxide band at the top and bottom of the 
atmosphere have been calculated for several atmospheric models and with various cloud heights. 
The variation of the fluxes with carbon dioxide amounts is determined, and its effect on tem- 
perature discussed. Plass’ estimate of a temperature drop of 3.8° C due to a halving of the carbon 
dioxide concentration appears to be too high by a factor of two or three. 


I. Introduction 


It has long been the desire of the writer to 
make a careful reinvestigation of the carbon 
dioxide theory of climatic change. It was felt, 
however, that it would be purposeless to 
reopen the subject unless a method of cal- 
culating radiative flux could be found that 
was of sufficient ascertainable precision that 
the results could be accepted as definitive. 
Such a method is now available (Kaptan, 
1959 a). It takes into account the detailed 
structure of the relevant infrared bands and 
its dependence on the thermal structure of the 
atmosphere. 

This paper reports the results of the appli- 
cation of the method to the problem of the 
effect of carbon dioxide variations on the net 
fluxes at the top and bottom of the atmosphere. 
The results are compared with. those resulting 
from the approximate method of Prass 
(1956 a). 


2. Calculation of the radiative flux 


The results reported here are part of a 
larger study of radiative flux in the 15 micron 
CO, band; a detailed account of the calcu- 


lation will be given in the report of the full 
investigation. The method described by the 
writer (KAPLAN, 1959 a) was followed quite 
closely, except that transmission in the vertical 
was converted to hemispheric transmission 
with the use of a transmission-dependent 
scaling factor. 

The 1,000 millibar level was taken as the 
bottom of the atmosphere and the 100 millibar 
level as the top of the atmosphere. The net 
upward flux was calculated at both of these 
levels as a function of thermal structure and 
cloud heights. The bases and tops of clouds 
were chosen in hundred millibar intervals; 
and the lapse-rate of temperature was consid- 
ered constant with height. It was assumed that 
the ground temperature was the same as the 
adjacent air temperature. Four atmospheric 
models were used: the combinations of ground 
temperatures of 0°C and 40°C with lapse-rates 
of approximately 5.6°C and 11.2°C per kilo- 


I I 
meter. Actually the values of : and - were 


used for the exponent x in the equation for 
constant lapse-rate: 


(T/T) —(p/po)* (1) 
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Transmissions were calculated at five-wave- 
number intervals from 565 to 820 cm-!. Thus 
the band considered extends from 562.5 to 
822.5 cm! (12.16 to 17.78 microns) except 
that a ten-wave-number interval at the strong 
central Q-branch was excluded. After the 
transmissions were obtained, the radiation 
integral was evaluated by Simpson’s rule with 
the atmospheric layers chosen thin enough that 
halving them affected the results by only a 
small fraction of a per cent. Some of the results 
are given in spectral detail elsewhere (KAPLAN, 
1959 b). 

The standard amount of CO, in the at- 
mosphere was chosen as 260 cm when reduced 
to S.T.P. This is in line with CALLENDAR’S 
(1958) estimate and only slightly larger than 
KEELING’S (1958). Calculations were made for 
this “standard’”’ amount and for ten per cent 
less and ten per cent more. Calculations were 
also made for half the amount, but only with 


I 
i x corresponding to a lapse-rate of about 


5.6° C per kilometer. 


3. Results 


The results of the calculations are given in 
Tables 1 and 2, in units of ergs cm? sec"? 
(one c.g.s. unit = 0.002064 cal. cm? day-1). 
Table 1 represents the outgoing radiation at 
the top of the atmosphere except for the 


small errors introduced by the truncation of 
the atmosphere at 100 millibars and the 
neglect of the Q-branch region. The errors 
from both of these causes become negligible 
when differences with concentration are taken. 
The errors in Table 2 from both of these 
causes are negligible whether or not differ- 
ences are taken. The pressures given in col- 
umns I refer to the top of an overcast 
in Table 1 and the base of an undercast 
in Table 2. In comparing one column with 
another it should be remembered that the 
temperature height curve for Tgrouna = 40° C 


üb: 
andy ==interseots, that ford Toma = 0 C 
3 


I 
and x = e at about 400 mb. 


The solid lines of Fig. ı show, for the four 
models, the increase in radiation leaving the 
top of the atmosphere due to a decrease of 
10 % in the CO, content as a function of the 
level at which an overcast top occurs. The 
1,000 mb level represents cloudless conditions. 
The dashed lines show the decrease in radiation 
due to an increase of 10 % in CO, concen- 
tration. The solid and dashed lines of Fig. 2 
have the same meanings with regard to the 
net upward radiation at the ground as a 
function of the level at which an overcast 
cloud base occurs. The data for Figs. 1 and 2 
were obtained by subtracting relevant columns 
of Tables 1 and 2 respectively. 


Table 1. Net upward radiation at 100 mb from 15 u CO, band, in c.g.s. units 


| x= 2; lapse-rate & 5.7° C per km 


| 
Lx = = lapse-rate & 11.4° C per km 


Surface 
Temp. (°C) | 2 42 2 4° 

2aral CO: 130 | 234 | 260 | 286 | 130 | 234 | 260 | 286 | 234 | 260 | 286 | 234 | 260 | 286 

(cm STP) 

Cld Top 

(mb) 

Cloudless 57037\54074|53540|53058\91434|86535|85659|84870|47568|46948|46391|74024|72982|72046 
900 54640|51933|51443|51000/88303|83780|82966)82234|43124|42591|42111)68242|67328|66503 
800 51972|49541|49099!48699|84772|80661|79918|79248|38394|37949|37547|01955|61172|60467 
790 48980|46844|46455|46102|80756|77096|76432|75831|33356|32999|32676|55091|54446153861 
600 45578|43759|43427\43124|76109|72946|72369|71846|28014127743|27498|47596|47088146627 
500 41660|40178[39906| 3965917065 5|68029|67549|67113/22386/22196|22024 39418 39047|38709 
400 37062|35929|35723|35535|64085|62048|61675|61335|16574|16454|16346|30569|30324|30101 
300 31534|30762|30622|30495|55932|54527|54269|54036 10786|10721|10663|21204|21065|20940 
200 24600|24214|24143|24079]452506144544144413144294| 5445| 5421] 5400 11800|11746|11696 
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Net upward radiation difference at I0Omb(cgs. units) 


[0] 
200 300 400 500 600 700 800 900 100C 
Pressure at top of overcast (mb) 


Fig. 1. Increase (solid lines) or decrease (dashed lines) in radiation leaving top of 
atmosphere in 15 micron CO, band due to decrease or increase, respectively, of 
CO: concentration by ten per cent. 


If the sky is partly covered with clouds, or Table 3. Assumed distribution of cloud cover, as 
if the cloud tops or bases are not uniform, the viewed from below and aboye 


fluxes or flux differences can be obtained by a | 
3 5 c : Pressure (mb) |900,800|/700|/600|/500|400|300|200 
linear averaging with respect to fractional sky | 
cover. An “average” cloud distribution was | 
used to estimate the effect of clouds on the pires és À SU dot ee ee 
overall flux differences; this distribution is = 
given in Table 3 in units of fortieths of the Fortieths of | | 
entire sky. The values were obtained by minor COHAREPEEN 5 = 13 ES Da A ae eee 
Table 2. Net upward radiation at ground from 15 u CO, band, in c.g.s. units 
I I 
| RP lapse-rate & 5.7° C per km M = lapse-rate & 11.4°C per km 
Surface 
Temp. (°C) 4 “> 4 2 
Total CO 
(cm STP) 130 | 234 | 260 | 286 | 130 | 234 | 260 | 286 | 234 | 260 | 286 | 234 | 260 | 286 
| | | | | | 
Cld. Base 
(mb) 
900 3517| 3260| 3213| 3172| 4872| 4475| 4404| 4341| 6380| 6290| 6209 8787| 8650| 8525 
800 6540} 5994| 5896] 5808| 8985| 8r41| 7991| 7856|11453|11268|11107 15664|15380|15123 
700 9525| 8667| 8515| 8377|13022]11707[11475|11267|16109|15831 15580|22008|21579|2119 
600 12568|11385|11175|10988|17153|15339|15021/14736 20486|20117|19786|28051 a Bu 
500 15750]14219|13950[13708|21501\19150|18741|18375|24610|24155|23747|33865|33157|32523 
400 19152|17244|16908|16609|26197|23260|22752|22298|28448127911 27429|39439|38597|37845 
300 22875|20547|20140|19776|31420|27824|27205|26653|31894|31283 30734|44660/43691|42826 
200 27076|24267|23778)23341|37463|33096|32347|31680134751|34074|33468|49269|48184|47216 
100 32048|28655|28066|27540|44945|39599|38687|37875|36651|35030135283|52679|51503|5045 5 


Cloudless 38232|34064|33343|32699|55687|48832|47668|46634 37154|36420|35761|53790|52581|51504 
_— [1 27127119 st 7199794187 9718093413 74 94130420/35792|53790|52591 515041 
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Fig. 2. Increase (solid lines) or decrease (dashed lines) of net upward flux at ground 
resulting from decrease or increase, respectively, of concentration CO; by ten per cent. 


Table 4. Increase in upward radiation in c.g.s. 
units if assumed total CO, path of 260 cm STP 
is reduced by 50 % or 10% or decrease if path 


is increased by 10%, for assumed cloud distri- 


bution 
| at 100 mb | at ground 
Total CO, path 
(cm STP) 130| 234| 286] 130} 234] 286 
x LEC) 
+ (6) 2751| 421) 381/2960| 438| 391 
40 4618| 704] 635|4801| 700} 622 
5 0 | 148 399 498| 447 
40 765| 689 808| 721 
modifications of the data of TELEGADAS and 


Lonpon (1954). No attempt was made to 
take into account the temperature dependence 
of cloud height; instead, global averages were 
used. The resulting flux differences are given 
in Table 4. 

It is interesting that the total flux change at 
the bottom of the atmosphere is very close to 
that at the top of the atmosphere when the 
cloud distribution is taken into account. This 
implies that any heating or cooling that is 
induced by changes in the CO, concentration 
would be quite uniform with height. The 
dependence of the total flux changes on 
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stability is also remarkably small. It should be 
noted, however, that the numbers will tend 
to zero as the stability tends to zero or infinity. 

The effect of CO, changes on the fluxes 
increases with increasing temperature. The 
magnitude of this variation with temperature 
is such that an induced temperature change 
could be expected to be almost independent 
of the temperature itself. 


4. Discussion 


The striking feature of the results given in 
Tables 1, 2 and 4 is the relative insensitivity of 
the fluxes to the CO, concentration, compared 
with their sensitivity to the thermal structure 
of the atmosphere. Together with the apparent 
constancy of CO, amount in the free air 
(KEELING, 1958), this fact is very promising 
for the problem of determination of thermal 
structure from remote satellite soundings 
(Kaptan, 1959 b). The effect of varying the 
concentration on the fluxes is, in fact, con- 
siderably less than that estimated by Prass 
(1956 a). Plass found a change in the downward 
flux at the surface of the earth of 0.0125 cal 
Cmia*tseor4) OS 72010. 6.SRUNItS; tas ta result of 
halving the CO, amount in a clear atmosphere, 
assuming a surface temperature of 15°C anda 
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lapse-rate of 5.5°C per km. Our result, with a 
lapse-rate of 5.6°C per km, is 8020 c.g.s. units 
for a surface temperature of 40°C and 4,890 
c.g.s. units for a surface temperature of 0° C: 
The value obtained by linear interpolation to 
15° Cis 6,060 c.g.5. units, indicating an overesti- 
mate by Plass of about 44 per cent. Since 
Plass’ method of calculation could be reason- 
ably assumed to be the best hitherto in use, the 
discrepancy is an indication of the danger 
involved in using a method of calculating 
infrared flux that depends heavily on the use 
of empirical homogeneous path data and on 
asymptotic expressions for the transmission. 
Since the time of Plass’ work on this problem, 
better laboratory data on CO, transmission 
have become available; but it is not easy to 
estimate how much of his error is due to poor 
basic data because the curve-fitting involved 
in his method is too subjective to be repro- 
ducible. 

It should be noted that calculations involving 
differencing of the fluxes tend to have much 
larger errors than the errors in the flux calcula- 
tions themselves. Semi-empirical methods are 
therefore inadequate for the problem of deter- 
mining radiative temperature changes in the 
atmosphere. 

For average cloudiness conditions the change 
in downward flux at the surface is given by 
the last three columns in Table 4. For a surface 
temperature of 15°C, halving the CO, content 
results in a decrease in downward flux of 
about 4,110 c.g.s. units, or 0.00589 cal cm? 
min !, If we accept as reasonable Plass’ 
estimate that a 1°C change in temperature 
would result in a change of 0.0033 cal cm? 
min“! in the average outgoing radiation, this 
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would correspond to a temperature decrease 
of only 1.8°C.! This quantity would be further 
decreased if the shielding effect of water vapor, 
which is far from negligible, is considered. 
Although Prass (1956 b) realized that the 
existence of clouds would decrease the effect 
of changing CO, concentration on temperature, 
his clear sky estimate of 3.8° C for a halving of 
the CO, content was used in his discussions 
of the influence of CO, changes on climatic 
change (Prass, 1956 b, 1956 c, 1959). It now 
appears that this number is too high by more 
than a factor of two, and perhaps by as much 
as a factor of three. It would seem, therefore, 
that CO, variations could not play a major 
role in the ice-age cycle unless the CO, 
changes were by an order of magnitude. 

It is also found from Table 4, with the same 
reasoning, that a change of 10% in CO, 
content should result in a change of not more 
than one-fourth of a degree in the surface 
temperature. The magnitude of this change 
does not seem to be large enough to account 
for the secular increase in European tempera- 
tures that was observed in the beginning of 
this century. 
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pCO, in Sea Water and its Effect on the Movement 
of CO, in Nature”: 


By JOHN KANWISHER, Woods Hole Oceanographic Institution 
Woods Hole, Massachusetts 


(Manuscript received August 31, 1959) 


Abstract 


A method is described for measuring pCO, in sea water. A gas phase is analyzed con- 
tinuously by infrared absorption for CO, while it is equilibrated gently with the water 
in a countercurrent column. It has been used to determine the changes in pCO, produced 
by variations of temperature and total CO,. Partial pressure shows large changes for 
small increments in these two independant variables. These properties of sea water are 
useful in estimating the movement of CO, between the atmosphere and oceans. It appears, 
for instance, that most of the fossil fuel CO, released by man has been effective in in- 


creasing the percentage of this gas in air. 


Such diverse problems as biological produc- 
tivity, climatic change, and carbon dating are 
linked to an understanding of the movement 
of carbon dioxide in nature. Pertinent to all of 
these is a knowledge of the concentration of 
this gas in air and its exchange between the 
oceans and the atmosphere. 


Transport within the atmosphere and oceans 
is chiefly by mixing processes. But the move- 
ment of CO, across the sea surface depends 
on diffusion. Since tension or partial pressure is 
the driving force of diffusion, an understanding 
of air-ocean exchange of CO, must include a 
quantitative knowledge of the factors such as 
temperature and total CO, which effect the 
tension of this gas in the sea. 


1 Contribution no. 1103 from the Woods Hole 


Oceanographic Institution, Woods Hole, Mass. 
2 Supported by NSF grant G-4813. 
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When CO, from the atmosphere dissolves 
in sea water it enters the following chain of 
reactions: 


| A=CO, 
CO, + H,O = H,CO, = HCO,- = CO, 


The sums of these forms will be called total 
CO, in this paper. The partial pressure or 
tension exerted by the physically dissolved CO, 
is denoted by pCO,. Harvey (1955) has 
reviewed the considerable work of Buch and 
others that has been done on the chemistry 
of the CO, system in sea water. The amounts 
present as bicarbonate and carbonate ions can 
be found from pH, temperature, and salinity 
with the aid of the apparent dissociation con- 
stants for carbonic acid determined by Buch 
in sea water. Henries law is not valid for the 
relationship between pCO, and total CO, since 
only that component present as a physically 
dissolved gas is responsible for the tension. 
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Buch has computed a series of tables in which 
pCO, and total CO, can be determined from 
pH, temperature, and salinity. The agreement 
between pCO, computed from these tables 
and that measured directly has been rather 
poor if one considers the limited range of 
surface sea water (BUCH 1939, DEACON 1940). 
Part of the cause is due to the inability to 
measure pH with sufficient precision. An error 
of .o2 pH units will change the estimated 
partial pressure by 10 percent. Direct tension 
measurements may be equally inaccurate due 
to the difficulties of measuring a component 
that makes up only 1 part in 3000 of a gas 
mixture. Buch’s tables are, in addition, in- 
convenient for deriving the relationship of 
pCO, and total CO, and also the effect of 
temperature on pCO, when total CO, remains 
unchanged. This is due to the pH changes that 
occur. It was the need for these two relation- 
ships that resulted in the present work. 


There is an additional reason for directly 
measuring the CO, partial pressure properties. 
Buch defined apparent dissociation constants, 
K’ı and K’2 for the carbonic acid system in sea 
water using concentrations of ions rather than 
activities. When a comparison is made with the 
known dissociation constants it is found that 
the carbonate ions exert an activity equal to 
only 2 percent of their measured chemical con- 
centration, (The Oceans, page 205). Some 
recent work of GARRELLS et al (1960) shows 
that this is due to complexing of CO, ions, 
primarily with Mg and Na. Thus the actual 
situation in sea water can be considered: 


Ky Ky 
CO, + H,O = H,CO, = HCO,- = 
Ki Ki 
CO,= + Mg = MgCO, + Nat = NaCO,- 


Other cations have only a minor effect. There 
is a small amount of complexing involving the 
bicarbonate ion. 


In addition, the first dissociation constant, 
K’,, is defined using the total unionized CO, 
for the undissociated carbonic acid. Actually 
only about 1/1000 of this is H,CO 3, so Ky 
is 1,000 times smaller than the true K,. It 
seemed possible that these approximations could 
have an effect on the dynamics of the CO, 
system in sea water. 


Since this work has been finished a paper by 
Bolin and Eriksson has appeared that treats 
theoretically much of what I have done ex- 
perimentally. The agreement with my mea- 
sured values is rather good which indicates that 
the possible effects mentioned above are not 
drastic. Their work will be referred to further. 


Method: A method has been developed for 
the continuous recording of CO, partial pres- 
sure in sea water. The apparatus is shown in 
Figure 1. A Fernbach flask containing about 
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Fig. 1. Diagram of apparatus for pCO, determination. 


3 liters of sea water is surmounted by a gas 
equilibrating column. Water is pumped from 
the flask to the top of the column where it 
flows down by gravity over a packing of Berl 
saddles and re-enters the flask. The air stream 
enters at the bottom of the column and leaves 
at the top. Both countercurrent flow and a 
large surface facilitate tension equilibration 
between the air and the gases in solution. The 
entire apparatus is placed in a temperature bath. 

A closed 500 ml. volume of air continuously 
circulates through the column, a drying tube, 
infra red CO, analyzer, and pump. Since the 
ratio of sea water to air volume is 6:1, very 
little CO, must be supplied by the water to 
bring the air into equilibration at the normally 
occurring values of pCO,. The three liters of 
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pCO, IN SEA WATER 


sea water contain about 135 ml. of total CO,. 
The amount in the air at 300 ppm is only 
.I percent of this. A correction has been com- 
puted at the higher tensions. 


The analyzer is a dual path type with a CO, 
filled pneumatic detector, (model 200, Mine 
Safety Appliances Co.) The beam chopper is 
quartz which does not absorb appreciably in 
the infrared until 4u. Thus only the longer 
wave length absorption bands are effective. 
This results in a 2,000 to 1 differentiation against 
water vapor. In spite of this the gas stream 
must be dried. If it were not, that amount of 
water vapor in 20°C saturated air would 
appear as a 20 ppm CO, signal. The infrared 
principle offers a 10 to 100 times increase in 
sensitivity over previous methods in determi- 
ning CO, in gaseous mixtures. 


It is necessary to frequently calibrate the ana- 
lyzer with gas mixtures of known CO, con- 
tent. These are made by volumetrically diluting 
100 percent CO, with nitrogen in a 20 liter 
stainless steel spirometer. The CO, is injected 
from a syringe. If pure gases are used the ac- 
curacy depends primarily on measurement of 
volume, with some errors introduced by tem- 
perature differences at which the volumes are 
determined. Given values could be repeated to 
within 1 percent. 

The amount of total CO, in the sea water 
has been varied in different ways. Distilled 
water was saturated with 1 atmosphere of CO, 
gas at a known temperature to produce a solu- 
tion of known concentration. This was checked 
by measuring its electrical conductivity. Small 
amounts of the solution were titrated into the 
3 liters of sea water from a syringe, while 
pCO, was continuously recorded. Introduction 
of equivalent amounts of CO, free water 
caused no pCO, shift. 

Biological material, fish and sea weeds, were 
also used to introduce known amounts of CO». 
Their CO, output rate was initially determined 
by enclosing them in a container and using gas 
analysis. The material was then introduced into 
the sea water for varying lengths of time to 
provide the desired CO, increase. Photo- 
synthesis of sea weeds, measured by oxygen 
increase, was used to remove known amounts 
SHlO, 

The increase in CO, from biological activity 
and the normal gas exchange of a stored sample 
Tellus XII (1960), 2 
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require that the water for pCO, determinations 
be immediately at hand if the value is to repre- 
sent that which occurred in the ocean. For this 
reason, I have been able to make measurements 
only on water from Woods Hole Harbor. This 
has a salinity of about 31.5 %. The buffering 
capacity of sea water for CO, will increase in 
proportion to salinity so the data given can be 
extrapolated. 


Results 


Figure 2 shows the large changes in pCO, 
that result when the temperature of a single 
water sample is varied. The total CO, remains 
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Fig. 2. Effect of temperature on pCO, in a local water 
sample. The cross points were taken with increasing 
temperature and the circles with decreasing to show 
that the effect is reversible. 


unchanged during a run since the system is 
sealed and the biological activity in the water 
sample is negligable. The CO, partial pres- 
sure at the sea temperature at the time of 
this run was 285 ppm. The air at this time had 
a concentration of 315 ppm. The sea was 
distinctly undersaturated with reference to the 
atmosphere. 
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The tension increases about 4 1/,%/° C. 
Harvey (1955 pg. 6) states that pCO,, in- 
creases 1 %/°C and this is what one finds from 
Buchs tables at constant pH. However, with 
constant total CO,, pH decreases about. .o1 
units for each °C increase. It is this pH change 
and its resulting effect on the various equilibria 
of the carbonate system that accounts for the 
larger observed sensitivity of pCO, with 
temperature. It is unnatural to speak of tem- 
perature variation with constant pH since this 
can only happen if the normal temperature in- 
duced pH change is buffered with the addition 
of acid or base. This latter is in effect a change 
in alkalinity which has been shown not to 
occur to any extent in sea water. Revelle and 
Suess arrived at a 7 percent change per °C in 
some undescribed manner. BOLIN and Erics- 
son unfortunately did not include the effect 
of temperature in their theortical analysis of 
the CO, system. 


Figure 3 is a plot of the temperature varia- 
tion of pCO, with varying amounts of total 
CO,. The parallel straight lines on a logarithmic 
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Fig. 3. Temperature variation of pCO, with different 
amounts of total CO,. 


plot show that the effect of temperature is 
constant in percent per degree regardless of 
the initial tension. The meaning of these curves 
is clear. If one is to measure pCO, meaning- 
fully to 1 percent then the temperature must be 
regulated to .2°C and this must be related 
with equal accuracy to the sea temperature. 


The variation of pCO, with changes in total 
CO, is shown in Figure 4. It is readily apparent 
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Fig. 4. Variation of pCO, with changes in total CO,. 
The 10° C values were taken on the same water sample 
after the 20. 


that sea water undergoes large changes in pCO, 
from the variations in CO, that are regularly 
produced by biological and physical processes. 
This increase in tension must be due to the 
increase in the physically dissolved component 
of the carbonate system. It can be computed 
(using the solubility of CO,) that about a 
1/13 of any CO, added to sea water stays as a 
physically dissolved gas. Fresh waters low in 
alkalinity should show at least 10 times as great 
a tension change since most of the CO, stays 
in solution. My determinations on distilled 
water have confirmed this. 


Discussion: 


The measured CO, partial pressure proper- 
ties of sea water can be applied to some con- 
siderations of the movement of CO, in nature. 
The curves in Figures 3 and 4 can be used in 
estimating, for example, the amount of CO, 
that must enter or leave a water mass to keep 
pCO, constant when the temperature is varied. 
The important point to consider is the large 
sensitivity of tension to changes in both tem- 
perature and total CO,. 


Different levels in the oceans show large 
variations in CO, partial pressure. Intermediate 
ocean depths are typically low in O, due to 
biological activity and are thus high in CO,. 
These layers. may be nearly anaerobic in the 
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tropical Pacific and Atlantic and the partial 
pressure will therefore be as high as 2,000 to 
3,500 ppm. I measured a pCO, of 3,000 ppm 
in the nearly anaerobic bottom layer of a pol- 
luted harbor at Woods Hole. The presence of 
such high values is probably of little significance 
to the air-ocean mixing of CO, since waters of 
this type rarely reach the surface. The high 
pCO, layer is separated from the mixed surface 
layer by a temperature produced density 
difference which inhibits vertical mixing. As 
small amounts do mix across this thermocline 
they are diluted into the mixed layer with 
relatively little effect. Since the higher CO, is 
accompanied by nutrient salts it may be largely 
fixed by photosynthesis before the elevated 
tension can be relieved by CO, diffusing into 
the air. This general argument may be valid 
for all cold bottom water formed at high 
latitudes which moves at depth to the tropics. 

The surface water does show more limited 
variations in pCO, which are unrelated to 
mixing or to temperature change. COOPER 
(1937) found that these changes in the English 
Channel accompany the seasonal variations in 
photosynthesis and respiration of marine orga- 
nisms. I have found such changes in the coastal 
water at Woods Hole. The partial pressure is 
below atmospheric during a plankton bloom 
when photosynthesis is fixing CO», and is 
higher than the air when the organic material 
is decaying. One might expect the largest 
excursions in Artic waters during the long 
hours of daylight. Buch has found a pCO, of 
150 ppm in the summer near Spitzbergen. 

It appears that biologically produced pCO, 
gradients at the sea surface are larger than those 
due to seasonal temperature changes. From the 
data presented here a 20° seasonal temperature 
shift requires about 1.5 ml/l difference in CO, 
to keep pCO, constant. The seasonal photo- 
synthetically fixed CO, may be several times 
this. The two effects partly nullify each other 
since the water has more CO, made available 
from spring warming at the very time that 
photosynthesis is maximum. Unfortunately 
this reasoning depends largely on scanty data 
from coastal areas rather than the open sea. 

In theory a knowledge of pCO, in space and 
time along with pertinent data on the physical 
state of the sea surface could be used to calculate 
net CO, fluxes. In practice this would be very 
difficult. However, surface pCO, does give 
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qualitative information about where and in 
what direction the gas is diffusing. Isotope data 
provides an independent means of arriving at 
an average rate over a long period. Craic 
(1957) calculated that a given CO, molecule 
has an average residence time of 7 years in 
the atmosphere before diffusing into the ocean. 

There has been considerable discussion of the 
effects and fate of the CO, that has been releas- 
ed by man’s burning of fossil fuels. This has 
equaled more than 14 percent of the amount 
already in the atmosphere, (REVELL and Suzss 
1957). Bouin and ERIKSSON (1959) have pre- 
sented a comprehensive discussion of the 
problem and have shown that most of the in- 
crease must still be in the atmosphere. The 
problem is essentially one of estimating how 
much of the potential increase has been buf- 
fered by the seas capacity to absorb CO,. The 
pCO, properties presented here provide a basis 
for estimating this. BOLIN and Eriksson used 
a theoretically derived figure. 

Any higher atmospheric pCO, must ap- 
proach tension equilibrium with the seas and 
in the process some net CO, will pass into the 
water bringing its partial pressure up and 
lowering that of the air. Since the oceans con- 
tain some 60 times as much CO, as the air it 
would not be surprising if they had the 
ability to buffer out a large part of any atmos- 
pheric increase. We need consider here though 
only that portion of the oceans in touch with 
the sea surface, the mixed layer above the 
thermocline. Wooster (CRAIG 1957) has esti- 
mated from oceanographic data that the 
average depth of the thermocline is 75 meters, 
therefore less than 2 percent of the ocean water 
is in diffusion equilibrium with the atmos- 
phere at one time. This contains only slightly 
more total CO, than is in the atmosphere. 
Referring to Figure 4, a TO percent increase in 
partial pressure is produced by a .6 percent 
increase in the total CO,. Thus one finds that 
less than 10 percent of any atmosperic increase 
would be removed by the amount of water in 
the mixed layer at one time. Since the same 
water does not remain in the mixed layer due 
to continual exchange with bottom water, 
some estimate of the rate of turnover of the 
mixed layer is needed. From oceanographic 
current data it seems unlikely that the mixed 
layer could have been renewed more than 2 
to 5 times in the last so years. It is clear even 
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from these rough figures that any CO, released 
from burning fossil fuels has been largely 
effective in increasing the concentration of this 
gas in the atmosphere. Over longer periods of 
time, hundreds to thousands of years, all of the 
deep water will equilibrate at the surface and 
something like 90 percent of any atmospheric 
increase will be absorbed. However, man is 
burning fuel at an increasing rate so if one 
considers only the solubility effect described 
here, the CO, should continue to increase. 

Bolin and Eriksson have presented a mathe- 
matically formalized version of this argument 
using the same two-layer model of the sea. 
They show that the amount of CO, increase is 
quite insensitive to a wide range of turnover 
rates for the deep ocean. This is a consequence 
of the small amount of water in the mixed 
layer aided by the poor buffering capacity for 
CO,. Even if one assumes a short turnover 
time, which means the surface layer has been 
renewed many times, the total amount of 
water passing through it can absorb only a 
small part of any atmospheric CO, increase. 
CALLENDAR (1958) concluded from a considera- 
tion of actual measurements that there has 
been a 10 percent atmospheric increase. The 
above argument support this view. 

REVELL and Suess (1957) arrived at the much 
lower estimate of 3 to 6 percent for the 
atmospheric CO, increase due to fossil fuels. 
This difference is largely due to their using a 
model of the sea which assumes it is a well 
mixed reservoir. This is difficult to rationalize 
with our present oceanographic knowledge. 

Fossil fuel CO, contains no CM so its in- 
introduction decreases the specific activity of 
the atmospheric CO,. The effect has been 
measured and found to be less than a few per- 
cent. This might at first seem incompatable 
with an atmospheric CO, increase of several 
times that. This can be explained by considering 
both the dynamic exchange and the tension 
equilibration taking place in the surface layer. 
It has already been mentioned that the water 
in the surface layer at one time can absorb 
less than 1/,, of any atmospheric increase. But 
if diffusion back and forth across the sea sur- 
face is rapid enough, the surface layer CO, 
would completely mix with that in the air. The 
Suess effect decrease would be about 50 percent 
since these are nearly equal reservoirs. Craic’s 
(1957) estimate of 7 years as the average turn- 
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over rate of air CO, indicates this is nearly the 
case. His figure also demonstrates that tension 
equilibration with the atmosphere is nearly 
complete for water passing through the mixed 
surface layer. Isotope data only gives informa- 
tion on the size of the effective reservoir, inde- 
pendant of hydrographic considerations, with 
which the atmospheric CO, has been mixed. 

The long term fate of the present or any 
other CO, increase is complicated. The solu- 
bility buffering described here will reduce the 
present increase to less than 3 percent. The 
slightly greater total CO, in the sea water will 
increase the solubility for CaCO . Sea water 
is in contact with this in bottom sediments and 
in animals such as corals. Bolin and Eriksson 
show that this increases by 3 times the capacity 
of sea water for CO,. Thus the eventual in- 
crease left in the atmosphere will be less than 
1/4) of the original. 

Further compensation of the greater amount 
of CO, may result from an increase in photo- 
synthesis both on the land and in the sea. This 
would result in some of the newly released CO, 
being fixed as organic material. The increased 
Ca ion might concievably accelerate the forma- 
tion of animal and plant CaCO,. In the event 
of no further CO, disturbance these inter- 
dependant factors would eventually produce 
what we can consider the long time average 
value of atmospheric CO ;. The tendency 
towards equilibrium of carbonates with CO, 
(Urey, 1952) must be only one of many factors 
determining this value. 

The regulatory mechanism just described 
would have a time constant of at least many 
centuries so oscillations of such periods might 
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Fig. 5. Diagram abstracted from Redfield to show that 
large amounts of CO, have moved through the ocean- 
atmosphere system in geological time. The quantities 
are expressed on the basis of atmospheric CO, equals tr. 
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be expected. Figure 5 shows the amounts of 
carbon in different forms in nature expressed 
in units of atmospheric CO, (REDFIELD 1958). 
It is clear that in geological times large amounts 
of CO, have passed through the atmosphere- 
ocean pool. Since this flux has involved such 
cataclismic events as volcanic activity there 
must have been variations in the amount of 
CO, in transit. 

Prass (1959) and others feel that a rise in 
the earths surface temperature will result from 
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more CO, in the air because of the increased 
absorption of thermal radiation which would 
otherwise be lost to space. It is tempting to 
attribute such things as climatic variations in 
historical and geological times to such a factor. 
The complicated atmospheric dynamics in- 
volved force one to view a cause and effect 
relationship of CO, and climate with caution. 
Fortunately man is providing his own experi- 
ment to test this hypothesis, (REvELLE and 
SUESS). 
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Abstract 


From the CO -data of the years 1955—59 for Scandinavia the annual cycle was determined. 
The picture is changing from year to year showing time displacements of the maxima and 
minima, which can, however, be placed “in phase” by shift of the time axis. One gets a mean 
value of the annual cycle of the years 1955—59 with two maxima and two minima. The total 
annual mean value obtained for Scandinavia is for 1955 326 ppm, 1956 321 ppm, 1957 323 


ppm, 1958 315 ppm, 1959 331 ppm. 


Introduction 


In most investigations on atmospheric CO;- 
content the interest has been concentrated on 
seasonal variation, where the usual summer 
minimum is assumed to be associated with the 
assimilation of vegetation. 

There are very few data on the complete 
annual variation of CO,. Furthermore they 
are difficult to interpret due to the fact that 
during the vegetation period considerable 
fluctuations have been observed between day- 
and night values (Huser and POMMER, 1954). 
This phenomenon can be observed most 
markedly in the first few meters of the at- 
mosphere. As most investigations are confined 
to this layer, they can be used for studies of 
the CO,-content in the “free” atmosphere 
only to a limited extent. This is particularly 
evident, if one considers that the difference 
between the day- and night values can be as 
large as 100 ppm while at the same time the 


1 Presented at the “Conference on atmospheric chem- 
istry at the International Meteorological Institute in 
Stockholm, 20—22 May 1959”. 


2 Part of the research reported on this document has 


been sponsored by the U. S. Weather Bureau under 
contract No. Cwb—9759. 


Fig. 1. The net of stations. 
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contribution effect of the vegetation on the 
whole atmospheric CO,-distribution amounts 
only to about 2 % (KALLE, 1943). 

Since 1955 there exists in Scandinavia a 
network of about 19 stations,! where air 
samples are collected three times a month 
(ERIKSSON, 1954) (FONSELIUS, KOROLEFF and 
Bucu, 1955) (Fig. 1). The data obtained are 
published in Tellus regularly. A condensed re- 
port of the results obtained during the first 
year was published by Fonsgttus, KOROLEFF 
and WÄRME (1956). These data show consider- 
able seasonal variations but good agreement be- 
tween the stations, except for the Finnish ones. 
In the present paper an investigation of the 
data from 19$5—59 is made and critically 
reviewed. 


Accuracy of measurements 


Already at the sampling errors can be in- 
troduced through air respired by the man 
doing the sampling. This danger is greatest in 
calm weather. Also in calm weather respiration 
from the vegetation during the previous night 
can spread into the lower layer of the air 
whereby an entirely wrong picture of the 
CO,-content of the free atmosphere is ob- 
tained. Samples, which are taken during calm 
weather have therefore received particular 
attention in the evaluation. 

Leaky gas burettes may be influenced by 
outside air. At the sample collection the burette 
is, therefore, provided with overpressure. If 
this overpressure is noticed when analyzing 
the sample there has been no contamination 
during transport and storage. In the chemical 
analysis an accuracy of about + 3 ppm can be 
expected. Since May 1959 a gas analyser 
(IRGA) has been used instead of the chemical 


analysis and has an accuracy of measurement of 


= I ppm. 


Evaluation of the data 


The data obtained by analyses are normally 
plotted chronologically as Fig. 2 shows, making 
it possible to check the reliability with refer- 
ence to the errors discussed earlier. Samples 
taken at calm weather and from burettes 
without overpressure are denoted by open 
circles. 


1 Inthe Scandinavian network the Dutch station VI 
(Vlieland) is also included. 


Tellus XII (1960), 2 


7—002421 


The data scatter considerably, which, in 
many cases, may be due to errors; the extreme 
values are, indeed, not always identical with 
the “‘reliable” ones. For the values in the figure 
it is therefore necessary to find a curve rep- 
resenting the annual cycle. This has been 
done for the stations Fl and Br in Figs. 2 and 
3. A seasonal variation with maxima and min- 
ima at different dates from year to year can 
be observed. From this smoothed curve the 
monthly mean values are obtained and shown 
in Figs. 4 a and b—now for all stations. For a 
period of one year one obtains generally two 
maxima and two minima,which by displacement 
of the time axis can be put “in phase” (Figs. 5 
and 6). In some cases the periodicity is inter- 
rupted by small irregularities; sometimes a 
maximum or minimum is very weak. This 
is particularly the case in the autumn, but the 
periodicity is, in general, resumed, thus 
making the reconstruction of the annual cycle 
possible. The values of 1959 make a general 
exception. Here an independent maximum 
between the existing spring- and autumn 
maxima has appeared. Again the deviating 
tendency of the Finnish data is striking as 
already shown by FonsEttvs et al. In the further 
discussion the Finnish data are therefore treated 
separately from those of the remaining Scan- 
dinavian stations. 

In Figs. 5 and 6 the time adjusted curves 
from Figs. 4 a and b for different stations are 
shown for each of the years together with 
their average. With the exception of the values 
from 1959 rather uniform curves with two 
maxima and minima are obtained. These 
maxima and minima, consecutively marked 
I—2—3—4, are represented in Fig. 7, in which 
the time of the maxima and minima for every 
station is entered. The annual mean curves 
from Figs. 5 and 6 are collected in Figs. 8 and 9. 

In order to compare the annual mean curves 
with each other a representation as in Fig. 10 
was chosen. First the deviations of the monthly 
mean values from the yearly means have been 
computed and the average of these deviations 
for each month. This again has been done 
separately for Finland and for the remaining 
Scandinavian stations, and cover the years 
of 1955—59. The striking fact about this 
figure is that the Finnish data now can be 
fitted into the remaining Scandinavian ones, 
if the time axis is displaced about 3 months. 
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FLAHULT 


Figs. 2 and 3. The distribution of 

CO;-data for several years at the 

station Fl (Flahult) Fig. 2, and Br 
(Bredkälen), Fig. 3. 


As it would be very interesting to connect 
these annual variations with climatological 
events, a program for this purpose is at present 
under preparation but so far consistent results 
are not yet available. The monthly mean 
values in Figs. 4 a and b have been used to 
compute the annual mean values for every 
station found in Fig. 11. By graphical inte- 
gration one obtains an annual mean value for 
the whole area. A comparison of these and 
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1959 


averaged monthly means obtained by graphi- 
cal integration are shown in Table r. 


Table 1. Annual mean values of CO,-content in 
ppm 

1955 1956 1957 1958 1959 

Scandinavia without Finland! 319 320 323 308 326 

Jette) ain eae SR ES ae Re 337 326 325 335 340 

The whole Scandinavia?.... 326 321 323 315 331 


1 From monthly mean values. 
? From the graphical integration. 
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Fig. 4 a. The monthly mean values for all stations. 
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Fig. 4 b. The monthly mean values for all stations. 
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Fig. 5. The annual cycle of all stations without the Finnish, a—b_. marks the dis- 
placement of time axis for getting the curves “in phase”. 


Discussion of the results 


as Comparison with other investigations 


There are, as already Stated, very few in- 
vestigations ranging over a whole year. In 
STEINHAUSER’s (1958) data one can find an 
annual cycle for 19 57 showing a winter maxi- 
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mum in December—January, and a summer 
minimum in May—June (Bien 12). 
(1959) Presents a collection of a great number 
of European investigations in which he 
transforms the monthly mean values into 
Percentage of the maximum month. He 
receives a winter maximum in January and a 
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Fig. 6. The monthly mean values for the Finnish stations. 


summer minimum in August. EFFENBERGER 
(1951) used partly the same data as Bray; 
he determined the deviations of the monthly 
mean values and obtained an annual cycle 
with maxima in January to March and besides 
that in August to October, as well as minima 
in June to July and in December (Fig. 10). 
Here, indeed, a good coincidence with the 
Scandinavian investigations is found. 

Huser and POMMER (1954) have used 
continuous measurements over a period of one 
year. The annual cycle, however, becomes 


unrecognizable because of the big amplitudes 
between day- and night values. 


b. Time sampling with respect to the daily cycle 

of the CO;-content 

In the Scandinavian investigations the sam- 
pling time is chosen to be 1300 hr local time, 
but Huber’s and Pommer’s investigation 
points to the necessity of taking the daily 
cycle into consideration. The daily cycle 
shows, as is well known, a considerable CO,- 
excess near the ground during the night. This 
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Fig. 7. The date of the maxima and minima marked with 1—2—3—4 
as shown in Figs. 5 and 6. 


excess is consumed by the “CO,-hunger” 

of the assimilating vegetation in the course of 
the forenoon. In case the consumption of CO, 
cannot be covered sufficiently by the night 
excess, CO, will be taken from the atmosphere. 
According to measurements performed by 
Huser (1951) over a corn field a well-develop- 
ed daily cycle still can be observed at a height 
of 22 m. In the data by CHapman et al. 
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Figs. Sand 9. The annual mean values collected from the 
representation in Figs. 5 and 6 by putting curves “in 
phase”. The time for displacement of time axis is shown 


at the lower picture. 
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Fig. 10. Comparison of deviations of monthly mean 

values in relation to respective annual mean content for 

the Scandinavian stations and for Finland separately. 

The time axis for the Finnish data is displaced about 3 

months. The ‘‘Europe” data are taken from EFFEN- 

BERGER 1951, Fig. 1, and correspond to selected data 
from about 50 years. 
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Fig. 11. The distribution of annual mean values for every station. The annual mean value for the whole Scandi- 
navian network is found by graphical integration, 


(1954) the influence of assimilation at a height 
of 152 m—also over a corn field—is still 
noticeable. 

A result of vertical sampling during flight, 
shown in Fig. 13, indicates the influence from 
ground up to a height of 800 m, thus to the 
presumptive height of the atmospheric bound- 
ary layer in which the turbulent exchange is 
still effective. The daily cycle of CO, is, of 
course, influenced by the extent of the tur- 
bulent exchange of the lower atmosphere. In 


the Scandinavian project this effect has been 
taken into consideration by chosing 1300 LT 
for the collection of air samples, as at that 
time the exchange of the ground air is best 
developed. Moreover, in the choice of stations 
the proximity of excessive vegetation and 
industry has been avoided. 

Some of the results from an investigation in 
the vicinity of Stockholm are here of interest. 
It is seen from Fig. 14 that at the two collection 
points, a local factor is clearly present in spite 
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Fig. 12. The annual cycle for 1957 at Wien, from STEINHAUSER 1958, Table 1. 
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of the small distance. Here “golfbanan”’ signifies 
measurements at a golf course, about 10 km 
north of Stockholm, “bryggan” those at a 
place only 1 km east of the golf course at the 
shore of a narrow bay of the Baltic with an 
easterly wind. In order to avoid the influence 
of vegetation as much as possible, a series of 
experiments was also done on a small rocky 
island situated in the archipelago on the east 
coast of Sweden. Again an easterly wind was 
blowing. Fig. 14 shows the results. Here, the 
nightly respiration from vegetation is missing, 
but the day values increase somewhat. 

Later measurements at the golf course on 
July 28, September 1 and December 2, 1959 
show a displacement of the day minimum to 
a later hour and decreasing amplitude. 


Fig. 13. Vertical CO:-distribution at a flight 28.10.1959, 
1400 LT. 


Tellus XII (1960), 2 


daily samples. 
samples taken 3 times every month. 


c. Consideration of single samples for ascertaining 
monthly mean values 


The small choice of samples—about 3 a 
month—brings up the question how far 
these represent the monthly mean value. To 
investigate this a series of measurements on 
daily samples taken at the golf course may give 
some information. In Fig. 15 the daily data 
are entered. Out of these current values the 
following ones: 


328, 332, 323 ppm for the 11th, rath, 13th 


Golfbanan 1959 
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Fig. 14. The daily cycle of CO;-content at the golf 

course (golfbanan) on several occasions and comparison 

with a place 1 km away at a Baltic creek (bryggan), 
and at a small island (Dalarö). 
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Fig. 15. Daily CO;-data, 1300 LT from a golf course 10 km north of Stockholm. 


of August and 316, 314, 312 ppm for the 
15th, 16th, 17th of September 


corresponding to the appointed days of sam- 
pling in the Scandinavian network give the 
monthly mean values for August 328 and for 
September 314 ppm. The correct values are 
for August 317 and for September 310 ppm. 
Even more easy to compare are the single 
samples mentioned by STEINHAUSER (1958). 
In Fig. 12 the monthly mean values of the 
samples taken three times every month— 
the first, the tenth and the twentieth—are 
compared with the monthly mean values of 


the daily samples (as in Tab. 1 of Steinhauser’s 
paper). Although these monthly mean data 
deviate from each other up till ro ppm, the 
annual cycle can still be distinguished and does 
well agree with that using all data. 

Besides the present routine program in 
Scandinavia preparations are being made for 
continuous measurements at two or more 
stations in Sweden. Further investigations on 
the vertical distribution of CO,-concentration 
are also planned. 

The writer is indebted to Dr. E. Eriksson 
and Mr. S. Fonselius for many discussions and 
useful suggestions. 
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Abstract 


A method of making CO; standard gas mixtures varying between 1 and 400 ppm in concen- 
tration with an accuracy of + I ppm is described. An accurately known volume of pure CO: 
gas is mixed with CO: free dry air in an evacuated vessel. A method to take samples from 
this vessel and to analyse the samples in an IRGA gas analyser is described. The use of the apparatus 


for the Scandinavian CO: network is described. 


1. Description of the apparatus (Fig. 1). 


The apparatus consists of a round bottomed, 
thick walled Pyrex glass flask (A) of about 20 
liter volume, which can be evacuated safely. 
The volume of the flask is estimated by weighing 


Q 


with water. The flask is fitted with a rubber 
stopper (B), covered on the inside with a silicone 
lubricant. (Glass joint better.) Through the 
stopper are inserted four capiallry tubes (Cj, 
C,, C3 and C,). The first capillary tube (C:), 
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which has an inner diameter of 0,1 mm and 
10 mm? volume, and reaches to the middle of 
the flask (A), is connected to the gas burette (D) 
through the two way stopcock (E,). The burette 
(D) is graded in 100 parts of a ml and can be 
connected to the mercury bulb (G,) or to the 
gas cylinder (I) through the three way stopcock 
(F,). The CO, gas cylinders manufactured by 
Svenska Kolsyrefabriken, Stockholm, contain 
99.5 % pure CO, and 0.5 % argon. The Hg 
bulb (G,) is connected to the open niveau tube 
(H,). The second capillary tube (C,) is con- 
nected to the gas pipette (J) (volume 250 ml) 
through the two way stopcock (E,) (similar to 
E,). By turning (E,), (J) can be connected to 
the gas analyser (K) through the tube (L). The 
gas pipette (J) is connected to the three way 
stopcock (F2) (similar to F,). The stopcocks of 
(J) are called (V) and (W). (J) can be connected 
to the mercury bulb (G,) (similar to G,), which 
is fitted with a level tube (H,), or to the free air 
through (F,). The tube (L) is fitted with a three 
way stopcock (M,) wich makes it possible to 
let in air in different parts of the system. 
Between (L) and the analyser (K), an air dryer 
(N) is connected, containing Dehydrite (Mag- 
nesium Perchlorate) which does not absorb CO, 
The tube (Y) of the gas analyser (K) is con- 
nected with the glass tube (P), which is fitted 
with a Hg manometer (R) and a stopcock (Q,). 
(P) is through the three way stopcock (M,) 
connected to the vacuum pump (O). (M3) can 
connect (O) with (K), with (A) through the 
capillary tube (C,) which is fitted with the 
stopcock (Q,) or with both. This makes it 
possible to evacuate the tube (Y) in the gas 
analyser (K) or the vessel (A) or the whole 
system. Through the fourth capillary tube (C,) 
the vessel (A) can be connected with the As- 
carite tube (S) fitted with the Calcium Chloride 
tube (T) or with the Hg manometer (U), using 
the two way stopcock (Es) (similar to E,). 


2. Mixing of Standard gas 


The Hg bulb (G,) is lowered so that the Hg 
level is on the height of (F,) with the connection 
between (D) and (G,) open. No air is allowed 
to stay below the stopcock (F,). The rubber 
tubing which connects the gas cylinder (I) with 
(D) is removed from (F,) and (F,) is turned to 
connect (D) with the free air. A glass vessel is 
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kept under the tip of (F,) and by help of a 
rubber bulb some air is blown through (D) in 
order to remove Hg in the stopcocks. The 
rubber bulb is connected to the free end of the 
stopcock (E,). Then the bulb is removed, (Ej) 
is still kept open to the free air, (I) is connected 
to (F,) and CO, gas is blown through (D) 
during several minutes. Then the gas stream is 
stopped and (E,) is closed. (F,) is turned so that 
(D) is in connection with (G,) and (G:) is raised 
until the Hg level in (H,) reaches the zero mark 
on the burette (D). Then (E,) is opened some 
seconds in order to let the CO, in (D) reach 
room pressure. The levels of (H,) and (D) are 
adjusted to exactly the same height and the 
exact reading of the level in (D) is taken. The 
burette is kept in room temperature for one 
hour and the eel is adjusted to room pressure, 
and the reading is checked. During this time 
the vessel (A) is evacuated to vacuum (less than 
t mm Hg) through the capillary tube (C;) with 
the stopcock (Qs) open and the stopcock (Mg) 
connected both to (C;) and (P). The stopcock 
(Q,) is closed during this operation. The vacu- 
um can be controlled on the manometer (R). 
The stopcocks (E,), (E:) and (E;) have of 
course to be kept closed. After about 15 minutes 
the flask (A) is evacuated and the stopcock (Q;) 
is closed and the pump stopped. Then CO, 
from (D) is introduced in (A) by careful opening 
of (E,) so that (D) is connected to (C,). When 
the Hg level in (D) rises, the bulb (G,) is raised 
simultaneously in order to control the volume 
in (D). When the level in (D) reaches the 
desired value or better when it is a little below 
it, the stopcock (E,) is closed and the level in 
(D) is adjusted to room pressure by adjusting 
(G,) and comparing to the level in (H,). 

By careful operations of (E,) the introduced 
amount of CO, can be adjusted to exactly the 
desired volume. After this operation the stop- 
cock (Eg) is opened in order to let in CO, free 
air through the absorption tube (S). A drying 
tube is connected to (S) in order to protect the 
Ascarite from moisture. The flask (A) is filled 
with air in about half an hour. The pressure is 
then checked with the manometer (U) by 
turning (Es) to connect (A) to (U). The temper- 
ature in (A) is then checked on the thermometer 
(X). The flask (A) now contains a known 
amount of CO, in a known volume at known 
temperature and room pressure. The pressure 
is read off on an ordinary barometer. 
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3. Analysis of standardsamples with the IRGA 
The Hg bulb (G,) is raised so high that with 


(F2) in connection to (J), Hg enters into the gas 
pipette (J). The stopcock (W) in the pipette (J) is 
closed and the stopcock (V) is opened. (E,) is 
turned to connect (J) with the tube (L), (M,) 
is turned to connect (L) with (N), (Q,) is 
opened, (M) is kept as before and the system 
is evacuated with the vacuum pump (O). The 
vacuum is checked on (R). Then (Eg) is turned 
to connection with (A). (J) is now filled with 
standard sample. (W) is opened and (G,) is 
lowered until the Hg level reaches (W). The 
stopcock (Q,) is then closed and the pump is 
stopped. (E,) is turned to connect with (L) and 
by raising (G,) the sample is brought into the 
tube of the IRGA analyser (K). The level in (J) 
is adjusted to 760 mm pressure by comparing 
to the level in (H,) using the difference between 
room pressure (from barometer) and 760 on 
a mm scale behind the tubes. The temperature 
in the IRGA is checked (in our apparatus it is 
constant at 32°C in a thermostated room). The 
reading on the voltmeter is noted. A second 
sample is taken from (A) as described above. 
The first sample gives often a somewhat high 
reading due to CO, in the capillary tube (C,) 
which has not mixed with the air. It is therefore 
important to have a very small capillary diam- 
eter in (C,) and to take several samples from 
A). 
making several standard mixtures with 
different concentrations of CO,, a standard 
curve can be drawn. Fig. 2 shows a standard 
curve made with our apparatus. It is a linear, 
when a suitable range is selected. Ordinary air 
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samples can then be analyzed in the IRGA using 
the same kinds of gas pipettes as (J). In the 
sampling procedure a small overpressure is 
given to the samples as described by Fonsezrus 
and KOROLEFF (1955). The system is evacuated 
to the stopcock (V), the overpressure is released 
through (F,) by opening (W) and then the Hg 
bulb is raised after turning (F,). (V) is pee 


and the sample is introduced as described above. 


4. The IRGA gas analyser 


The IRGA gas analyser model SB 1 with 
fluorite windows manufactured by Sir Howard 
Grubb Parsons & Co LTD, Walkergate, New- 
castle upon Tyne 6, England, has been found 
to be very suitable for our purposes. A Philips 
L. F. Amplifier-voltmeter GM 6017/02 has been 
used as indicator. The gas analyser will not be 
described in details here, a good instruction is 
delivered from the manufacturer. However 
some instructions have to be included. The 
construction of the analyser is schematically 
shown in the figure. (Y) is the sample tube 
which can be evacuated to vacuum and (Z) is 
the standard tube containing CO, free air. The 
analyser is used on the highest sensitivity and 
the zero reading on the voltmeter is adjusted 
with a standard gas containing 275 ppm CO. 
The analyser is then calibrated using different 
standard mixtures as described above. Every 
day the analyser is checked with a standard gas 
and adjusted to give the right reading by using 
the balance shutter. Several checkings with the 
standard gas are made during a day’s operation 
of the apparatus. 

Our equipment is now used for the routine 
analyses from the Scandinavian network of 
CO, sampling stations described by FONSELIUS 
and KOROLEFF (loc cit). It may be pointed that 
the IRGA gives results compared to dry air, 
while the earlier technique (FONsELIUSs and 
KORLEFF loc cit) gives the results compared to 
air saturated with water vapor at about 20° C. 
A correction to dry air is easily made when the 
temperature in situ and in the laboratory during 
the analysis is known. The accuracy of the 
infrared technique is somewhat better (around 
+ 1 ppm) than in the chemical method (+ 3 
ppm) and around 40—so analysis can be run 
in one working day. The use of the equipment 
requires however more training and experience. 
In order to test the accuracy of the method 
two test series of 8 samples were analyzed. 
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The test was made in the following way. An 
air pump was connected to a T-tube to which 
were attached two series of 4 gas pipettes con- 
nected to each other by small pieces of rubber 
tubing. An air stream was sucked through the 
gas pipettes during 1 hour on the roof of the 
Institute building. The pump was stopped and 
the stopcocks of the gas pipettes were closed. 
The samples were analyzed for their CO, con- 
centration. As shown in table ı the difference 
between the samples in each series was about 


+ I ppm. 


Table 1. 
Series I Series 2 
a | b a | b 
318 ppm 317 ppm | 315 ppm | 317 ppm 
318 » S TOME) DES Mu 3150 
318 » 319 » ZI 316 » 
Brom) STONE) ZU) 317) 
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Abstract 


The location of the cosmic ray equator off the west coast of Africa has been established from 
measurements obtained with a shipboard neutron monitor which made ten crossings between 
November, 1956, and October, 1958. From the individual determinations of the point of 
minimum cosmic ray intensity, the mean geographical position was 6.7+0.8° N at 14° W. 
Direct measurements of zero magnetic inclination were made at the same time in the immediate 
vicinity of this point. 

Although solar and geophysical activity varied considerably during this period, the location 
of the cosmic ray equator appears to have remained fixed within approximately 1°. However, 
the possibility of a small progressive change of about the magnitude of the uncertainty cannot 
be precluded. 

Comparison of the results with plots of measurements of the magnetic elements along the 
ship’s route reveals a significant difference between the location of the points of maximum 
horizontal intensity and minimum cosmic ray intensity. 


I. Introduction 


Ambiguities generally inherent in the map- 
ping of isocosms on the basis of data previously 
obtained by various investigators do not apply 
to determinations of the location of the cosmic 
ray equator. In particular, temporal variations 
in the cosmic ray intensity introduce a con- 
siderable uncertainty in the comparison of 
measurements obtained at different times, and 
the problem of correctly normalizing data 
obtained in a survey encompassing large di- 
stances is quite complex. Furthermore, it 1s 
well established that changes in the primary 
spectrum are produced by an externally-acting 
modulation mechanism. Thus, for example, it 
would be incorrect to assume that the locus of 
the so-called “knee” represents an isocosm at 
high latitudes (POMERANTZ, POTNIS and AGAR- 
WAL, 1959). 

It appeared reasonable to expect that these 
factors do not seriously influence determina- 


Tellus XII (1960), 2 


tions of cosmic ray intensity minima. The 
errors in the location of the cosmic ray equator 
introduced by the normalization procedure are 
probably smaller than statistical and systematic 
uncertainties. Therefore, in order to obtain 
additional evidence relating to the possibility 
that the trajectories of cosmic ray particles are 
significantly affected by distortions of the 
earth’s magnetic field by ionized interplanetary 
gas, repeated determinations were made of 
the location of the cosmic ray equator during 
a period characterized both by a high level 
and extreme variability of solar activity. 

The greatest discrepancy between the posi- 
tion of the cosmic ray minimum and the loca- 
tion of the equator based on the dipole repre- 
sentation occurs off the west coast of Africa. 
It was previously reported on the basis of pre- 
liminary results (POMERANTZ, SANDSTROM and 
ROSE, 1958) that at longitude 14° W the cosmic 
ray minimum occurs at approximately 15° N 
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Fig. 1. Routes of M/S Lommaren and M/S Stratus. Nov. 
1956—Jan. 1960 


of the conventional geomagnetic equator, and 
is inconsistent with other equators based upon 
the dipole approximation, but lies very close 
to the isoclinic equator representing a surface 
measurement of the dip angle of the earth’s 
magnetic field equal to 0°. Hence, this region 
appeared to afford optimum conditions for 
monitoring the cosmic ray intensity minimum 
over a protracted period of time in the hope of 
distinguishing between the two alternative ex- 
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planations of the discrepancy: (1) local magne- 
tic field anomalies and (2) effects of conducting 
interplanetary gas. 


2. Experimental Procedure 


The construction of the pile and counters in 
the shipboard neutron monitor designed for 
the present program was identical in every 
respect with that utilized earlier during cruises 
of HMCS Labrador and USS Atka (Rose and 
KATZMAN, 1956). The pile and all associated 
equipment were installed in a thermostatically- 
controlled air-conditioned hut which was car- 
ried aboard M/S. Lommaren from 15 Oct. 
1956, until 20 Jan., 1958, and was subsequently 
transferred to M/S Stratus, where it remained 
between 9 March, 1958, and 11 Feb., 1959. 
Fig. 1 shows the routes of the two ships. The 
station operated continuously during four 
cruises of Lommaren from Sweden to South 
Africa and return, and during two cruises of 
Stratus from Sweden to Australia and return. 
The southbound voyages of Stratus followed 
the same route as Lommaren. Thus, during the 
period November, 1956, to October, 1958, the 
neutron monitor traversed the equator off 
West Africa ten times, eight aboard Lommaren 
and two aboard Stratus. 


Oct.11, 1958 


GEOGR. LAT. 
0° 10° 20°S 


Fig. 2. Ship board monitor, M/S Stratus, equator crossing, 6 hour averages. 
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3. Results 


A. Determinations of Cosmic Ray Equator 


Fig. 2 is a plot of a typical set of data. Each 
of the points comprises a six-hour interval. 
The raw bi-hourly counting rates were cor- 
rected for variations in barometric pressure 
with the usual exponential factor and single 
absorption coefficient (— 0.96% /mm Hg) (Po- 
MERANTZ, POTNIS and AGARWAL, 1959). The 
data were entered on punched-cards, and data 
reduction including continuous internal con- 
sistency checks was performed by machine. 
All of the data were normalized for world-wide 
variations on the basis of the Uppsala neutron 
monitor records. 

For each of the crossings, a machine deter- 


n 
mination of the function y= >> a, x" best 
oO 


representing the points between 30°N and 
20° S was performed by the method of least 
squares. In every case the analysis was carried 
out for n = 2 and 3. In a few cases computa- 
tions were extended up to n = 5. The effect 
of the choice of # upon the position of the 
minimum was not appreciable. 

Similar treatment of data for crossing nor- 
malized on the basis of the Huancayo neutron 
wionitor records revealed that, within the 
resolution of the experiment, the position of 
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Table 1. Summary of Cosmic Ray Equator 
Determinations at Longitude 14° W. 


Location of | Location of 
Cosmic Ray | Dip Equator, 


Date Equator, 1955 Chart 
Geographic Geographic 
Latitude Latitude 
November 1956.. 7.5°N 8.3°N 
January, 1957... 8.4° N 
MATCH rose 6.4° N 
May, 1057-40 6.3° N 
‘July, 1057 ma LEN 
August, 1957.... 6.8° N 
November, 1957 . OEAN 
December, 1957 . 5.5°N 
APTI, TO58...... 6.2°N 
October, 1958 ... 6.9° N 8.6° N 
Mean Value ..... 6.7°+o.8° N 8.4° N 


U.S. Navy Hydrographic Office Observations of 


Zero Dip. 
; 5 Latitude 
Date Longitude] Latitude at 14° W 
March 9, 1956....| 12°%31’W|. 7° 28’ N|’ 7.0° N 
February 8, 1957..| 12° 35’W| 7° 32°:N| 7.0° N 


the minimum was not dependent upon the 
selection of the reference station, even when 
the most adverse conditions (greatest varia- 


bility) prevailed. 
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Fig. 3. The curves represent daily averages of the nucleonic intensity at Uppsala. The short horizontal lines indicate 
the periods of ship passages through Equatorial Region. 
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In six cases, the position of the minimum was 
the same for # = 2 and n = 3, and shifted 
less than 1° in the remainder. Table 1 gives a 
summary of the location of the cosmic ray 
equator thus determined for each of the cros- 
sings, with n = 2. The individual uncertainties 
are approximately 1°. The indicated standard 
error in the mean is computed from the ob- 
served scatter of the points under the assump- 
tion that the position of the cosmic ray equator 
is constant. The location of the magnetic dip 
equator according to the 1955 U.S. Navy 
Hydrographic Office isoclinic chart, taking into 
account the predicted secular variation, is 
also listed. 

The differences between the positions of the 
cosmic ray equator and the dip equator in- 
dicated on the chart appear to exceed the un- 
certainty in the former. The Hydrographer of 
the U. S. Navy has kindly provided data con- 
cerning two actual measurements of zero dip 
in the immediate vicinity of the cosmic ray 
intensity minimum. These observations are also 
summarized in Table 1. The small differences 
in longitude were taken into account by a slight 
extrapolation. The mean value of the cosmic 
ray intensity minimum coincides with the ob- 
served position of the directly-measured dip 
equator. 


B. Effects of Solar Activity 


It has been pointed out that the purpose of 
this program was to determine whether the 
position of the cosmic ray equator is dependent 
upon the level of solar activity. The observa- 
tions summarized in Table 1 extend over a 
period characterized by a profusion of specta- 
cular geophysical effects of solar origin, and 
great variability in the level of solar activity. 
Of particular relevance to the present discus- 
sion is the fact that the cosmic ray intensity 
fluctuated appreciably. Fig. 3 shows the daily 
average intensity of the nucleonic component 
during the periods of the ship’s passage through 
the equatorial region (i.e. geographical latitude 
30° N—20° S). It is evident that the present 
determinations of the position of the cosmic 
ray minimum represent widely different con- 
ditions in the interplanetary medium. 

Indeed, the nucleonic intensity displayed 
large short-term fluctuations, as well as con- 
siderable changes in base “level”. If inter- 


planetary gas clouds were responsible for both 
cosmic ray intensity fluctuations and distor- 
tions of the earth’s magnetic field at large 
distances, a variation in the position of the 
cosmic ray equator, correlated with the pri- 
mary flux, might occur. 

As seen in Fig. 4, the results are not incon- 
sistent with this possibility. The solid line re- 
presents the mean value based upon the as- 
sumption that the equator did not move, 
whereas, the dash-dot curve is the best-fit 
straight line determined by the method of 
least squares without this restriction. The latter 
has a slope of 1.5° per 12 % decrease in the 
nucleonic intensity at sea level at a high 
latitude. 


INTENSITY 
#600 4700 4800 4900 5000 5100 5200 


Fig. 4. Cosmic ray equator at 14° West ard nucleonic 
intensity, Nov. 1956— Oct. 1958. 


Consideration of all factors involved, how- 
ever, leads to the conclusion that the position 
of the equator remained constant within ap- 
proximately 1°, despite significant changes in 
the primary flux. In any event, it is clear that 
comparisons between single determinations of 
the position of the cosmic ray minimum by 
shipboard monitors are subject to considerable 
uncertainty. 


C. Comparison with Magnetic O! servations 


It is of interest to compare the location of 
the cosmic ray equator with magnetic observa- 
tions in the equatorial region along the route 
followed by the shipboard neutron monitor. 
Fig. 5 is a plot of the magnetic elements vs 
geographic latitude. The data from which 
these curves were prepared were kindly pro- 
vided by the Hydrographer, U.S. Navy. 

There is a significant difference between the 
location of the points of maximum horizontal 
intensity and minimum cosmic ray intensity. 
Although the tracks of the magnetic survey 
did not cross the ship’s route at the dip equator, 
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Fig. 5. Magnetic elements along the ships route. 
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the previously-discussed close correspondence 
of the zero vertical intensity and zero dip point 
with the cosmic ray equator is suggested. 
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Abstract 


Instead of the time-distance graph, the P-S diagram is effectively employed to estimate the 
value of Poisson’s ratio in the earth’s crust. This diagram is advantageous, especially to detect 
its local anomaly or distribution. The obtained results are as follows: 

1) Applying this method to eighteen earthquakes occurring in south-western Japan, the 
domain with anomalously large Poisson’s ratio was found in the Kyushu district. This hints 
to the possibility of an existence of a local magmatic reservoir. 

2) The azimuthal distribution of the crustal Poisson’s ratio in case of earthquake occurrence was 
investigated and it was found to be considerably different in the push- and pull-zone of initial 
seismic motion. This may be related to the state of stress accumulation causing the earthquake 


occurrence, 


I. Introduction 


In general, Poisson’s ratio in the earth’s crust 
is considered to be nearly constant. Its value is 
slightly less than 0.25 in the relatively shallow 
depths of the crust and is about 0.25 in the 
lower layers (GUTENBERG, 1955). Recently it 
was ascertained that the earth’s crust in Japan 
has a distinctly local character (NisHIMURA and 
KisHIMOTO, 1956). In the present article, this 
character is clarified in more detail from the 
investigation of Poisson’s ratio. And in the 
latter part, the distribution of the crustal 
Poisson’s ratio in case of earthquake occurrence 
is discussed. 


2. P-S diagram method 


Theoretically, Poisson’s ratio is related to the 
velocity ratio of the P-wave to the S-wave ac- 
cording to the following equation. 


N I 
Ve a SET 


where Vp and Vs are the velocities of longi- 
tudinal and transverse wave respectively, and 
o is Poisson’s ratio. 


In order to estimate the value of Vp/Vs, 
especially its local anomaly the P-S diagram 
method is effectively employed. In this diagram 
the arrival time of P-wave (Tp) is taken as 
abscissa and that of S-wave (Ts) as ordinate. 
The observed value at one station corresponds 
to one point. The value of Vp/Vs is given by 
the following formula 


Ph tS te 
Vs pee 


where Ty is the origin time at the hypocenter, 
being common for the P- and S-waves. Com- 
pared to the usual time-distance curve analysis, 
this method has some merits mentioned below. 

1) In order to construct the time-distance 
graph, the position of the hypocenter and the 
origin time at the hypocenter must be known 
and to estimate the value of velocity, some as- 
sumption about the subterranean structure is 
inevitably required, whereas in the P-S dia- 
gram analysis, only the origin time at the hypo- 
center is necessary and the value of Vp/Vs can 
be directly obtained without any assumption. 

2) Furthermore, in order to detect only the 
local anomaly or distribution of Vp/Vs, the 
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Fig. 1. The location of observation stations and epicenters. @ represents the observa- 
tion station and X represents the epicenter. The dashed line shows the horizontal 
contour of the region with anomalously large Poisson’s ratio. Three earthquakes 
A), B), and C) are mentioned in the section 3. 


origin time is not always necessary, and they 
are immediately obtained from the distribu- 
tion of points in this diagram. However, it is 
remembered that the value obtained by the 
P-S diagram is the mean value which is 
estimated along the seismic wave path from the 
hypocenter to the station. 


3. Local character of Poisson’s ratio 


In this section, the crustal Poisson’s ratio in 
south-western Japan is investigated, especially 
with regard to its local character. For this 
purpose, eighteen earthquakes are selected, 
twelve having rather deep foci (40—130 km 
depth) and the remaining six having shallow 
foci (10—30 km depth). The location of epi- 
centers and observation stations are shown 
in Fig. 1. In the following, three examples are 
mentioned. For details, reference is made to 
the paper by KAMITSURI (1959). 


Earthquake A) 1953 Jan. 234 115 47™ 59 
(A =33.3°N 9=132.2°E, H=40 km) 


The points are separated into two groups. 
The upper group points are those of the sta- 
tions in and near the Kyushu district and the 
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Fig. 2. The P-S diagram of the earthquake on Jan. 23, 
1953. Vp/Vs (1.85) corresponds to o (0.29), and 1.70 
corresponds to 0.24. 
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lower group points are those of the stations in 
the Shikoku and Chugoku districts. (Fig. 2) 


Earthquake B) 1955 Mar. 284 18% 12™ 20° 
(A=29.4°N 9=130.1°E, H=60 km) 


A straight line through the points of stations 
in the Shikoku and Chugoku districts is drawn. 
The remaining points which are those of the 
stations in the Kyushu district deviate up- 
wards from this straight line. (Fig. 3) 


Earthquake C) 1957 Nov. 14% 22h 29™ 09° 
(A=34.0°N 9=134.2°E, H=40 km) 


All points are those of the stations in the 
Shikoku, Chugoku, and Kinki districts and 
they are completely arrayed on a straight line 
within the limit of accuracy. (Fig. 4) 

In all cases, the points of stations in and near 
the Kyushu district deviate upwards, and their 
deviations amount to about 1.5 ~ 3.0 sec in the 
arrival times of S-waves. The shallow earth- 
quakes, the foci of which are situated in the 
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Fig. 3. The P-S diagram of the earthquake on Mar. 28, 
1955. Vp/Vs (1.75) corresponds to & (0.26). 
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Fig. 4. The P-S diagram of the earthquake on Nov. 14, 
1957. Vp/Vs (1.74) corresponds to o (0.25). 


upper layers of the crust do not show such an 
effect as seen in the above-mentioned cases. 
From these, it is concluded that the crustal 
Poisson’s ratio in the Kyushu district (0.28 ~ 
0.30) is anomalously larger than that in the 
neighbouring districts (0.22 ~ 0.25) and 
that the anomalous domain is situated at the 
depths from about 20 km to 40 km. The 
horizontal contour of this domain is roughly 
shown in Fig. 1 with the dashed line. 

Several years ago, one of the writers 
(NISHIMURA, 1950) reported that from the long 
period tiltmetric observation at Makimine 
(A = 32°38°N, @ =131°24’E), the earth 
tidal amplitude of the M,-component shows a 
large fluctuation beyond the observation error, 
and that this is attributed to the variation 
of properties of subterranean material as- 
sociated with the seismic activity in the sur- 
rounding area. Furthermore, the Kyushu 
district is one of the most volcanically active 
regions in Japan. By more extensive study from 
these various standpoints, it will be expected 
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Fig. 5. The push-pull distribution of initial seismic 
motion of the Daishoji-Oki earthquake. 


to clarify the existence of the domain with 
anomalously large Poisson’s ratio in more 
detail. 


4. Azimuthal distribution of Poisson’s ratio 
in case of earthquake occurrence 


It is generally accepted that earthquakes are 
caused by stress accumulation in some limited 
region of the earth’s interior and subsequent 
release of such energy according to the occur- 
rence of the fracture-like event. The detailed 
process is, however, far from being completely 
understood. 

As described previously, the P-S diagram 
method is very advantageous to detect any 
local distribution of Poisson’s ratio. In the 
present section, this merit is profitably used 
to investigate the azimuthal distribution of 
Poisson’s ratio in case of earthquake occur- 
rence. 

For this purpose, three shallow earthquakes, 
the foci of which are situated in the earth’s 
crust are examined, one being accompanied by 
aftershocks, and the other two belonging to 
the earthquake swarm occurring in the Shikoku 
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district. The report by KAMITSUkI (1959) may 
be consulted for details. 


1) The Daishoji-Oki earthquake 
1952 Mar, 07¢ 16h 32™ (1 = 36° 45 N. 
@ = 136° 20° E, H = 20 km, M = 7) 


The push-pull distribution of initial motion 
is seen in Fig. $ which shows the so-called 
quadrant-type pattern. The P-S diagram is 
shown in Fig. 6. 

In this figure, the points are arrayed on the 
two branched lines. The points on the upper 
line are those of the stations in the pull-zone, 
while those on the lower line are of the stations 
in the push-zone. The gradients of these two 
lines are 1.70 and 1.90 respectively and the 
former is a little larger and the latter is much 
larger than the normal value of 1.67 
(YosHIYAMA, 1957). Reflecting upon the con- 
struction of the P-S diagram, the gradients of 
these two lines represent the apparent values of 
Vp/Vs. If the origin time point is determined, 


40 


Fig. 6. The P-S diagram of the Daishoji-Oki, earthquake, 
(Solid circle: push, open circle: pull). 
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the true value of Vp/Vs can be obtained. As- 
suming the origin time at the hypocenter being 
common for the P- and S-waves, the corre- 
sponding point slides on the 45°-gradient line 
through the origin of this diagram. From the 
time-distance curve analysis, the origin time 
is presumed to be plotted on the segment 
between two branched lines. Accordingly, the 
value of Vp/Vs increases from about 1.60 to 
1.68 in the push-zone while it decreases from 
1.80 to 1.68 in the pull-zone. 


2) The Tokushima earthquakes 
1955 Apr. 104 115 o9 128 (A = 34.1° N 
œ@ = 134.2° E, H = 10 km, M = 4.6) 


1955 May 184 138 39™ 245 (A = 33.8° N 
@ = 134.3° E, H =0~20 km, M = 4.9) 


In the Tokushima prefecture, Shikoku 
district, small earthquakes have frequently oc- 
curred. Their magnitude are reported as from 
about 4 to 6. Two examples are shown in Figs. 


Fig. 7. The P-S diagram of the earthquake on Apr. 
tod rrh 09m 128, 
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ig. 8. The P-S diagram of the earthquake on May 
18d 13h 39m 245. 


7. and 8. An effect similar to the case of the 
Daishoji-Oki earthquake is found. 

From the above-mentioned examples, it is 
found that the crustal Poisson’s ratio in case of 
earthquake occurrence is considerably different 
in the push- and pull-zone of initial motion. 
As no local anomaly of Poisson’s ratio is found 
in the region concerned in the present investiga- 
tion, the above-mentioned effect must be attrib- 
uted to the variation of the crustal Poisson’s 
ratio which is related to the stress accumulation 
causing the earthquake. 


5. Summary 


In order to investigate the crustal Poisson’s 
ratio, especially its local anomaly or distribu- 
tion, the P-S diagram method is effectively 
employed. The advantage of this method is 
that the anomaly of Poisson’s ratio can be 
directly detected without any assumption about 
the subterranean structure. 
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As a first treatment, the local character of 
the crustal Poisson’s ratio was investigated. In 
the Kyushu district of south-western Japan a 
certain limited domain with anomalously large 
Poisson’s ratio was found, and from the anal- 
ysis of earthquakes having various focal depths 
it is estimated to be situated at the depths 
from about 20 km to 4o km. 

Recently it is stated that low-velocity layers 
with large Poisson’s ratio exist in the earth’s 
crust and upper mantle. From the results of the 
present investigation, it is ascertained that there 
exists a layer with anomalously large Poisson’s 
ratio in particular regions such as the Kyushu 
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district. But the correlation between such 
anomaly of Poisson’s ratio and the velocity 
variation with depth should be investigated 
more in the future. 

Secondly, the azimuthal distribution of Pois- 
son’s ratio in case of earthquake occurrence 
was investigated and it was found that the 
crustal Poisson’s ratio is considerably different 
in the push- and pull-zone of initial seismic 
motion. Especically this effect is distinctly ob- 
served in the domain close by the epicenter. 
This may be related to the state of stress accumu- 
lation causing the earthquake. 
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EDITOR S NOTE 


Observational data, which since 1957 have been published in Tellus under the headings “Current data 
on the chemical composition of air and precipitation” and “CO, values in Scandinavia” will from this 
issue be printed as a separate report, appearing once a year. Orders for this report should be mailed to: 
The Chemical Laboratory, International Meteorological Institute, University of Stockholm, Lindhagens- 


gatan 124, Stockholm K, Sweden. 
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